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ABSTRACT 
— ^ 

The NASA Skylab Sensor Performance Evaluation task 
SPE— S193— 004 is concerned with estimating the precision and 

I 

accuracy with which the S-193 Scatterometer measured the back- 
scattering cross section of ground ; scenes . These estimates 
were derived from data collected during Skylab missions. For 
this study, homogeneous ground sites were selected and S-193 
Scatterometer backscattering cross section data analyzed. 

The precision was expressed as the standard deviation of the 
scatterometer-acquired backscattering cross section. In 
special cases, inference of the precision of measurement was 
made by considering the total range from the maximum to min- 
imum of the backscatter measurements within a data segment, 
rather than the standard deviation* For Skylab missions 2 
and 3 a precision better than 1.5'dB is indicated. 

The indication of the measurement accuracy was derived 
from various comparisons. A theoretical scattering formula, 
most suitable to the surface model, was selected. Ground 
parameters were used to evaluate the theoretical values of 
backscattered cross sections of homogeneous sites. Aircraft- 
acquired backscattering cross sections were analyzed to verify 
and supplement the theoretical values. Through this tedious 

i 

procedure, the most appropriate set of backscattering cross 
sections was generated for certain sites. As a final step, 
the differences between the actual ■ measured values and those 
developed using aircraft-acquired data together with mathe- 
matical scattering models were computed. These differences 
were indicative of the accuracy of measurement. 


IV 



This procedure indicates an accuraty of better than 3 dB 
for the Skylab 2 and 3 missions. The estimates of precision 
and accuracy given in this report are for backscattering cross 
sections from -28 to 18 dB. Outside this, range the precision 
and accuracy decrease significantly. 
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'1,0 INTRODUCTION 

The goal of this investigation is to' provide a statis- 
tically based estimate of the precision and accuracy with 
which the S-193 Scatterometer measured the scattering cross 
sections of the ground scenes. 

The term precision is used to imply repeatability of 
data from sample-to-sample with no regard to the bias between 
the true value and the measured value of the scattering cross 
section. Thus, precision is of significance to investigators 
who are interested in differences between scattering cross 
sections for various ground scenes. The term accuracy, on 
the other hand, is used to imply a measure of the bias errors, 
plus the repeatability. The accuracy estimate is of impor- 
tance to investigators who utilize the absolute value of the 
scattering cross sections for correlations with a phenomenon 
of interest . 

The classical method for determining precision/accuracy 
would be to subject the system to a known environment and 
compare its output to a known standard. Additionally, an 
error analysis would be performed to place an upper bound on 
the measurement error of the system. However, accurate 
standard instruments are not available for S-193 Scatterometer 
data comparison. Original test data is also sometimes inade- 
quate for placing the necessary bounds on the system parameters 
required for the classical analysis. This is due , to the 
difficulty of knowing (and being able to simulate) the exact 
thermal environment that the S-193 Scatterometer experienced 
during the Skylab data-gathering missions. 
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Therefore, the estimates contained herein of the 
precision/accuracy of the scatterometer system, for the 
actual data-gathering periods, were, based on the sensor 
analysis and comparisons of the actual S-193-acquired data 
with values of the backscattering cross sections obtained 
by aircraft sensors and the cross sections computed by 
analytical methods (using actual ground data) . 

The sensor specification control documents were studied 

r.. . 

to arrive at the latest configuration of the scatterometer 
system. The system performance has been summarized in 
appendix A. The realized antenna performance is also noted 
in this appendix. The relationship between input power to 
the antenna and the system output recorded on the tape was 
obtained by simulating the elements of the scatterometer 
system on a computer. Several parameters were changed to 
study their influence on the output. These results have 
been discussed in appendix A. The sensor mathematical model 
and the preflight systein parameters gave the baseline preci- 
sion value for the scatterometer measurements. 

Since the estimates of precision/accuracy are based on 
the analysis of processed data, it is important to have an 
understanding of the data processing algorithms and possible 
sources of error due to processing. In appendix B the review 
of Skylab and aircraft scatterometer data processing is 
presented. The processing of the S-193 Scatterometer- 
acquired data has been reflected in the interpretation of 
the precision and accuracy estimates derived from processed 
backscattering cross section data. 
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In determining precision and accuracy, the crucial step 
is to locate homogeneous sites. This insures that the vari- 
ance caused by the variability in surface roughness, biomass 
cover, intervening medium properties, and surface dielectric 
properties is minimum. The procedure for selecting homoge- 
neous sites involved examination c\f the roughness and dielec- 
tric parameters as well as the intervening medium parameters. 
The medium and surface parameters were either measured by 
ground based sensors and/or airborne systems. The aircraft- 
acquired photographic, microwave, and laser profiler data 
was used to categorize the homogeneity of the scene viewed 
by the S-193 Scatterometer . The Skylab-acquired S-190 
photographs, S-193 Radiometer, and S-194 Radiometer data 
supplemented the data derived from ground and aircraft 
measurements. The precision was then expressed as the vari- 
ance or standard deviation of the scatterometer backscatter. 

In special cases, inference of the precision of measurement 
was made by considering the total range from maximum to 
minimum of the backscatter measurements within a data segment, 
rather than the standard deviation. 

The determination of the measurement accuracy was 
accomplished by using various comparisons, including a theo- 
retical scattering formula most suitable to the surface 
model being selected. Ground parameters were used to 
evaluate the theoretical values of backscattered cross 
sections of homogeneous sites. Aircraft-acquired back- 
scattering cross sections were analyzed to verify and supple- 
ment the theoretical values. Through this tedious procedure, 
the most appropriate set of backscattering cross sections 
was generated for each site. As a final step, the 
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differences between the actual measured values and those 
developed using aircraft-acquired data, together with 
mathematical scattering models, were computed. These 
differences were indicative of the accuracy of measurement. 

It should be emphasized that any accuracy determination will 
be pessimistic because ‘it will not be possible to specify 
the upper bound in the errors involved in arriving at the 
mathematical scattering cross section values. This problem 
is compounded by the fact that data from the sources con- 
sidered are gathered fr.om various vantage points, e.g., on 
the ground, from aircraft, and from spacecraft. 

It should be noted that data for determining precision/ 
accuracy was gathered over the same or quite similar ground 
sites during Skylab 2, 3, and 4 missions. The data analysis 
leads to the variation of the precision/accuracy of the sensor 
during the period of Skylab missions. This provides a useful 
input to the design of future spaceborne active microwave 
systems: Exact determination of the precision/accuracy is 

not possible. The values obtained with the procedure given 
in this report will enable the sorting of erroneous data. 

In fact, a number of S-193 cross section measurement values 
were found invalid with this procedure. Furthermore, the 
investigations also revealed the reason for this invalid 
data. 

The study presented in this report is only for the 
Skylab 2 and 3 missions, A precision of better than 1.5 dB 
and an accuracy of better than 3 dB is shown for the scatter- 
ometer backscattering cross sections. 
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The review o£ the Skylab 4 data disclosed the presence 
of two problems: 

(1) The scatterometer received signal decreased by 
approximately 24 dB throughout the mission. 

(2) The roll scan operation began to malfunction 
on January 9, 1974. 

As a result of these anomalies, special processing 
procedures were implemented in the NASA production data 
processing program. It is anticipated that the precision/ 
accuracy of the scatterometer will, have degraded during 
Skylab 4 mission. However, the extent of this degradation 
cannot presently be. defined. . 
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2.0 GROUND DATA USED FOR COMPARISON WITH 
S-193 SCATTEROMETER MEASUREMENTS 


The estimation of the precision/accuracy of the 
scatterometer measurement requires a detailed knowledge 
of the ground scene. Proper consideration was given to 
the acquisition of the ground data during Skylab missions. 
These ground data requirements were specified in the Task 
Implementation Plan , (ref erence 1). The Sensor Performance 
Evaluation (SPE) sites were chosen so that the values of 
backscatter ing cross section obtained by measurement using 
ground based/aireraft systems and/or theoretically evaluated 
using pertinent ground parameters, could be compared with 
confidence with the S-193 Scatterometer measurements. In 
particular, the following criteria were used for the selec- 
tion of SPE sites (references 1 and 2) : 

• The site should be relatively uniform in such 
parameters as surface roughness and dielectric 
constant for at least one resolution cell. 

• The roughness range and dielectric constant of the 
surface should be known so that the backscattering 
cross section (oq) could be computed using theo- 
retical techniques. 

• The targets should be selected so that their a^’s 
cover the dynamic range of S-193 Scatterometer. 

For this, the high, medium, low, and no reflectivity 
sites were chosen. 

• The scenes over which data was previously collected 

were preferred over those for which no data 

has been reported in the literature. 


2-1 



• Sites for which simultaneous data was gathered using 
S-190A, S-193 Radiometer, and S-194 Radiometer 
sensors were preferred to ones where these sensors 
were not operated. 

• To reduce the effect of errors involved in computing 
exact attenuation because of heavy clouds and rain, 
no rain and cloud cover under 50 percent were speci- 
fied as unique (mandatory) test conditions. 

. ' * ‘ 

The SPE site criteria were satisfied by choosing ocean 
targets (high, medium, and low seas) and relatively smooth 
uniform ground targets (Great Salt Lake Desert and White 
Sands, New Mexico). Deep space was chosen to give the no 
reflectivity condition. 

Detailed ground data requirements were stated in the 
S-193 Quick-Look II Plan (reference 3). Subsequently, these 
requirements were revised to reflect available resources for 
the SPE ground truth effort. In this section, a description 
of the ground data which will be utilized in this report will 
be presented. Actual ground data will be given when S-193 
Scatterometer data compa'risons are illustrated. The term 
’’ground data" is used in this report for any data (other 
than that acquired by S-193 Scatterometer) which will be 
used in determining S-193 Scatterometer backscattering cross 
section precision/accuracy. In general, three types of 
ground data have been used to classify the scene sensed by 
the S-193 Scatterometer. These are discussed in sections 2.1, 
2.2, and 2.3. . 

There are several experimental results of the backscat- 
tering cross sections from ocean and ground scenes reported 
in the literature. Theoretical models have also been given 
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for scattering from these surfaces. These tlisoretical and 
experimental results have also been used to generate data for 
comparison with Skylab- acquired S-193 Scatterometer data. 

2.1 SKYLAB -ACQUIRED DATA 

The data from the following sensors was used to gain 
an understanding of the site: 

• S-190A photographs 

• S-193 Radiometer/Altimeter 

• S-194 Radiometer 

The photographs helped to verify cloud conditions, 
rain, and general ground features such as vegetative cover 
and open water bodies. The photographs will not be repro- 
duced in this report, but the results of the review of these 
photos will be utilized. 

The S-193 and S-194 Radiometer data has been used to 
verify the homogeneity of the sensed area. Additionally, 
these measurements were employed to calculate the reflec- 
tivity of various ground scenes. 

The S-193 Altimeter can yield the backscattering cross 
section, for the following pitch angles: 0°, 0.4°, 1.3°, 2.65°^, 

7.56°, and 15,6°. Where available, this data has been used 
for comparison with the S-193 backscattering cross sections 
for similar sites. 
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2.2 AIRCRAFT -ACQUIRED DATA 


The aircraft sensors specified for the collection of 
,the data vrere: 

• Metric camera 

• Laser profiler 

• Litton Industries LTN-51 Navigation Computer 

• 13.3 GHz Scatterometer (or 13.9 GHz Radiometer/ 
Scatterometer) ■ 

• Multifrequency Microwave Radiometer (MFMR) 

The data from metric camera, 13.3 GHz Scatterometer, 13.9 GHz 
Radiometer/Scatterometer , and MFMR were used similarly to 
that of the Skylab-acquired support data (see section 2.1). 
The laser profiler data was used to categorize the roughness 
scale of the SPE site. Roughness power densities and surface 
correlation functions were computed from the laser profiler 
data. These parameters are vital to the selection and eval- 
uation of the theoretical models of radar backscatter from, 
ground scenes. The Litton Industries LTN-51 system can be 
used to determine the surface wind velocity. 


In concluding this section, it should be noted that an 
anomaly (very high backscatter) was found in some portion, 
of the aircraft-acquired 13.3 GHz Scatterometer and 13.9 GHz 
Radiometer/Scatterometer data. However, questionable data 
was not used for comparisons in this report. 
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2.3 DATA ACQUISITION WITH GROUND-BASED SYSTEMS 

Direct measurements o£ environmental and surface 
conditions were also made within the limitations imposed 
by the necessary resources available to acquire and process 
the data. Fortunately, extensive weather data was compiled 
by the Space Meteorology Group (U.S. Weather Bureau). This 

I 

group is located at the Johnson Space Center and provided 
much of the weather information needed to plan Earth 
Resources Experiment Package (EREP)' passes during Skylab 
missions. The weather maps provided by the Space Meteorology 
Group present gross information about the ocean surface wind- 
speeds, cloud cover, and extent of rain on a global scale. 

The requirements for the ground-based data acquisition 
were as follows: 

Ocean sites 

— Radiosonde data 

— Water temperature and salinity 

/ 

— Sea state (significant wave height) 

— Sea surface wind velocity 

Land sites 

— Radiosonde data 

— General surface characteristics, 

including dielectric properties (from 
available geological information) . 

National Oceanographic and Atmospheric Administration 
(NOAA) has taken wind and wave measurements using aircraft 
and surface instrumentation. These measurements were not 
taken for Skylab 3 (SL-3) mission. All NOAA data was reviewed 
and only applicable data was used. 
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Wind and wave data was also Tecorded at Buccaneer Tower 
(by NASA/ JSC, Flight, Operations Division) off Galveston 
Island. This data was uj.ilized in this report. 

For the ground sites, previous measurements made by 
U.S. Geological Survey were reviewed. Ground conditions 
were also ascertained using photographs and laser profiler 
data. In certain cases the weather at the time of the 
Skylab overpass was verified by calling the weather bureau 
nearest the SPE site. 



3.0 THEORETICAL MODELS FOR RADAR CROSS SECTIONS 
FOR TERRESTRIAL SURFACES 

Numerous approaches have been advanced to explain 
scattering from rough surfaces . The mathematical models 
which have shown considerable prjomise in explaining the 
experimental observations have several constraints. In 
general, these constraints are (reference 2): 

• Most models assume a rough boundary between air and 
a homogeneous dielectric. Thus, the subsurface 
nonhomogeneities are not included in the models . 

• The radar wavelength is assumed much smaller than 
the dimensions of the radar resolution cell. 

• The radar system is assumed to be sufficiently far 
from the resolution cell. The incident wave at the 
surface of the site is assumed plane. Furthermore, 
the reflected wave at the radar is also assumed 
plane . 

• The rough surface is considered to have isotropic 

statistical properties over the radar resolution 
cell. ' 

.• The rough surface is considered stationary random 
process . 

• A uniform dielectric constant is assumed over the 
resolution cell. 

• Most models assume only two or three statistical 
parameters to describe the rough surface (standard 
deviation, mean slope, correlation distance , etc .) , 
In practice, these parameters are rather troublesome 
to obtain for terrestrial surfaces. 
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A. summary o£ the available methods for calculating 
for rough surfaces is given by Barrick (reference 4). The 
mathematical models for the computations fall in three 
categories : 

• Semiempirical models 

• Geometrical models 

• Statistical models 

3.1 SEMIEMPIRICAL MODELS 

Semiempirical models offer the simplest results and 
require little analytical derivations. These models involve 
one or more arbitrary constants which are determined from 
an agreement between the model and the measured results. 

These constants must be chosen for each class of rough 
surfaces. The most common semiempirical models use Lambert’s 
law : 

a = C, cos^ 6 (1) 

0 1 

and the generalized Lambert’s law: 

Qq = C 2 (cos0)^^, (2) 

where and C 2 are constants obtained by best fitting 

the equations to the measured data, as is the constant n . 

0 is the angle of incidence (see figure B- 3 , appendix B) . 
Both of these models apply to some terrestrial surfaces that 
scatter very diffusely (reference 4) . These laws do not 
apply to the sites chosen for Sensor Performance Evaluation. 
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3.2 GEOMETRICAL 'MODELS 

Geometrical models assume a surface composed of simple 
shapes arranged randomly on a planar area. From the knowl- 
edge of the scattering from simple shapes and proper boundary 
conditions, the field scattered from these surfaces can be 
calculated. These m^odqls partially take, into account the 
multiple scattering, i.e., mutual interaction among the 
simple shapes. The scattered field is usually easily cal- 
culable. Among the various shapes considered in the litera- 
ture are bosses on a conducting plane and infinite cylinders 
arranged with a random spacing upon a plane sheet (refer- 
ences 4 and 5) . Beckmann (reference 5) has calculated the 
scattering cross section of a rough surface by approximating 
the surface with infinite half -planes tilted at variable 
angles with respect to the horizontal plane. These planes 
are all arranged in one direction. Purely geometrical models 
will not be considered in this report because one can seri- 
ously question how accurately such a model could be applied 
to rough terrestrial surfaces with composite roughness. 

3.3 STATISTICAL MODELS 

The model most applicable to the surfaces selected for 
the S-193 Sensor Performance Evaluation is the one where the 
roughness is characterized by suitable statistical parameters 
Statistical models for a rough surface treat the height of 
the surface from the mean planar surface as a random variable 
The most commonly used height distribution function is the 
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Gaussian distribution function 


PC?') given by (reference 4): 


pco = 


hV27 


exp - 


2h‘ 


(3) 


where ? is the height from the mean planar surface and h 
is the standard deviation of this, height. The roughness in 
the horizontal direction can be described by introducing a 
surface correlation function. The correlation equation 
requires that the height of the rough surface above every 
point (x and y) in the mean surface plane be multiplied by 
the height above a point x distance away, and the 
product be integrated over x and y and divided by the area 
in the mean surface plane defined by the integration limits. 
The limits are permitted to become infinite (reference 6) . 

Mathematically, the correlation function can be expressed 
as (reference 6) . 


B(x) = lim 





?(x,y) ? (x+Ax,y+Ay)dT 


(4) 


-T 


T = [(Ax)^ + CAy)h^^^ 


where 
of the surface. 


and T denotes the extent 


3-4 


EEIPP'-OI^UCIBILITY 01*' Tlib 

ffrp 'X'SAL PAGE IS POOii 



I 


The surface correlation function has the following 
properties (reference 6) : 

• B(0) = h^ 

• As , the statistical dependence of the height 

at (x, y) and the height at (x+Ax, y+Ay) will 
decrease, i.e., lim “ mean surface height = 0 

= mean square surface slope = 

T = 0 

The surface can be categorized in terms of the root mean 
square (rms) height h (for the probability density func- 
tion of equation (3), the rms height is equal to the standard 
deviation of the surface height 5 from the mean planar 
surface 5=0 )• The surface for which h<<X(X is the wave- 
length of the radar system) is called ''slightly rough." 

"Very rough" surfaces have h>>X . 

Slightly rough surface scattering can be treated 
mathematical] V by us'ing perturbation techniques, whereas 
for the very rough surface, the analysis involves an 
asymptotic method. Presently, no mathematical methods have 
been developed for the surface where h is of the order of 
a wavelength. 

The choice of a correlation function for terrestrial 
surfaces has been a source of considerable controversy in 
the published literature on the subject of radar scattering 
(reference 2). Experimental determination of B(t) is 
cumbersome since it requires the knowledge of the surface 
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profile. Laser profiler data was requested as part of the 
ground data for the SPE sites. The laser data can be 
processed to yield the correlation function and the surface 
height distribution function. This is what was done. 

For the theoretical investigatioils correlation functions 
are usually chosen on two bases — first to allow complicated 

integrations to be carried out in the * analysis and second to 

- . 

yield a best fit between theoretical* and experimental values 
of the versus 0 curve. One widely used correlation 

function with correlation length Z is the Gaussian. 

B(t) = h^ exp (5) 

Gaussian-correlated surfaces have continuous slopes at all 
points with a total mean square slope of 



The exponential correlation function 

B(t) = h^expC- I t 1/il) (6) 

has also been used in many theoretical models (reference 5) . 

Surfaces with exponential correlation function are 
jagged and have vertical facets (jumps) (reference 4) . For 
such surfaces, the surface slopes and all higher order surface 
derivatives can be undefined or infinite at many surface 
points . 
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Previous investigations by Krishen (reference 7) have 
shown that neither daussian nor exponential correlation 
functions gives an accurate fit to the backscatter versus 9 
data from ocean surfaces. A better fit was observed with 
Gaussian correlation .function for the gravity waves and 
exponential correlation for the ca’pil’lary waves. For the 
land sites, the choice of a correlation function will 
entirely depend on the type of terrain. 

The correlation distance I defines the region of 

* < 

statistical dependence of the surface height. The surface 
heights of two points separated by a distance greater than 
i are essentially statistically independent of each other. 
Points within the distance I are statistically dependent. 
Thus, the correlation distance is intuitively the average 
horizontal extent of ripples or irregularities in the rough 
surface . . 

3.4 APPLICABLE THEORY TO SCATTERING FROM 
SPE LAND AND OCEAN SITES 

Three theories which have received attention and show 
promise of efficient interpretation of experimental data 
are the Kirchhoff method, the small perturbation theory, and 
the composite scattering theory. A discussion of these 
methods will be presented in this section. Most of this 
section has been taken directly from Kaufman's report 
(reference 2) . 
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3.4.1 


Radar Scattering Theories 


Kirchhoff method 


The field scattered by a rough surface Sp Cx' ,y' ,z') 

^1 

is formulated according to Huygen’s principle and is given 
by the Stratton-Chu integral: 


tp (X’ ,y ,z’) = i 5 . dt 


1 f 

w J 


[- jwy^Cn X H)\jj 


+ (n X X + (n • S) V\(j]ds (7) 

In equation (7), is the portion of the interface 

illuminated by. the electric field S and magnetic field 
H , and is contour of , The coordinate system is 

illustrated in fijgure 3-1* Cp and Pq are respectively the 
free space permittivity and permeability, and 


I e 


- jjc r 


with r = [(x* - x)^ + Cy' - y)^ + (z’ ' S)^] 


1/2 


( 8 ) 
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(B) 


Figure 3-1. — Scattering from the rough surface. [C(x,y,5) is 
any point on the rough interface ind PiCx’ ,y' ,z') is the 
observation point. (B) is the far field approximation 

of (A).] 


n is the unit vector normal to the surface reflecting the 
incident wave , a) is the radian frequency of the radar, 

k is the propagation constant, and ’ ^ is the rough scattering 
surface. 

At each point C of the illuminated rough surface the 
fields E and H are the- sum of an hpcident plane wave field 
( 2 ^ and and a reflected field ( 2 ^ and ^2) so that 

^ + ^2 

( 9 ) 

H = Hj * 

The fields on the surface are evaluated by tangent 
plane approximation. The assumption is that at the point 
C , the surface is locally flat. This approximation 
restricts the use of this method to surfaces where the 
radius of curvature of surface roughness is large compared 
with the radar wavelength. Under this assumption the 
incident and reflected fields can be related by using 
Fresnel reflection coefficients. 

Once the fields at the rough surface have been 
established, these fields may be entered into the scattering 
integral and the scattering integral adapted for the 
calculation of the mean power , l/ 2 { [<E2E2>] } 

scattered from the surface to the radar system at a distance 
Rg . The brackets <> indicate that an averaging process 
involving the statistical parameters of the surface must be 

A 

carried out, while E2 and E2 are respectively the 
complex magnitude (involving both magnitude and phase) and 
the conjugate complex magnitude of the scattered field 
received at the radar. 
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In carrying out the averaging of the .scattered power. 

It will be assumed that the surface height variation about 
the mean surface plane has the Gaussian distribution of 
equation (3) while t'he 'surface correlation function has the 
form of equation (5). A rather complicated integral is the 
result of the averaging process. Certain approximations 
required to facilitate its evaluation are: 

• £<<T, i.e., the correlation length of the surface 
roughness is much less than the dimension of the 
illuminated area. 

• Kaufman’s procedure (reference 8) requires also 
that all rough surface slopes not exceed 0.3 in 
magnitude and that the absolute- .'alue of the surface 
dielectric constant be much gi cater than 1.4 (at 
least 1 . ) • 

Using the relation (reference 9) , 


a 

0 


= 47TR' 


Average power scattered back to radar 
Power incident at rough surface 


( 10 ) 


three normalized radar cross sections, c , a and a „ 

’ oHH oVV oC 

‘^oHV *^oVH “ calculated depending upon the polariza- 

tion of the transmitted and incident waves . 


Barrick gives the result (reference 4) : 
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and 


'^oC " °oHV ■ ^oVH ~ ’ 

where S^, are the root mean square (rms) slopes in 
X and y directions, (x direction being taken to lie along 
the surface in the plane of incidence), R(o) is the 
Fresnel reflection coefficient evaluated at 9=0 , i.e., 
at the normal angle of incidence, R(o) is given by: 

R(o) = -■ 

1 

where e is the complex dielectric constant of the rough 

surface. For a surface with isotropic slope distribution 

2 2 2 2 

” ^y ~ ^ » where S is the total mean square 

slope given by = 4h^/il^ . 

Kaufman (reference 7) provides a more complicated but 
possibly more exact result since he does distinguish 

between ^qYV ‘^oHH ' ’ 


‘^oHH “ "T ^h ® ^ ^1 ® 


-tan 6 


%VV = jl V 8 8 ^ exp C13) 


'^oG " ^oHV " '^oHV 



«fproducibility of TxiM 
P/VGF is poob 
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where and are the Fresnel reflection coefficients: 


D _ e cos 9 - (e-sin^9)^^^ 
^ e cos 6 + (E-sin^e)^^^ 


and 


R = cos 6 - (£-sin^9)^^^ 
^ cos 9 + (E-sin^ey^^^ 


The factors and K 2 are functions of 6 , and are 

given by: 


= 1 


COS0 

Ce - sin'e)!/^ 


tan 0 


K. = 


£ (1 - e) sin 0 

(e - sin^0) (e cos 0~+ ['e - sin^0] 


(14) 


A few comments about the above equations are in order. 

For the case of a perfectly conducting surface, the Fresnel 

coefficients reduce to R, = 1 and R = -1 so that 

h V 

Kaufman and Barrick’s results are equivalent for such a 
surface. Essentially both sets of results can be interpreted 
as the radar cross sections available from a rough surface 
as the incident wavelength X approaches zero. Also, the 
cross-polarized cross section is identically zero in the 
zero wavelength limit. Kaufman (reference 8) has shown, 
however, that the cross-polarized normalized radar cross 
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section behaves in the following manner as X->0 


oC 


2 2 I h 

8-rr^h sin20 ^ 



( 15 ) 


Experimental data gathered over gently undulating 
natural surfaces shows good agreement with Kirchhoff’s 
scattering ijiodel up to the incidence angles of 25° for VV and 
HH polarizations, from equation (15) does not show good 

agreement with the experimental data. This model has there- 
fore been. used in evaluating S-193 Scatterometer precision/ 
accuracy at the lower incidence angles. 


Small perturbation method 

For surfaces where the rms height is smaller than the 
incident radar wavelength, a small perturbation technique 
can be applied. Two methods have been proposed in the 
literature. In the Bass and Bocharov method (reference 10), 
the homogeneous problem with boundary conditions on the 
rough surface is converted to a nonhomogeneous scattering 
surface with boundary conditions on a plane. For more 
detailed discussion and applicability of this method, the 
reader is referred to Krishen (reference 11) who applied 
this method to calculate scattering from a rough layer. 

Once the field on the surface is evaluated with this method, 
the field everywhere can be calculated using the Stratton-Chu 
integral . 

Rice (reference 12) gave an extension of Rayleigh's 
scalar solution for solving vector problems. The method is 
based on mode representation of the scattered field. The 
scattering rough surface 5(x,y) and the scattered field 
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are expanded into Fourier series. The coefficients of the 
scattered field are ^evaluated u§ing boundary conditions on 
the surface, The following assumptions are involved in this 
method: ... , 

• k C(x,y) < 1.0; i.e., roughness height is small 
compared to .wavelength X , 5(x,y) is the height of 
the rough surface above the ‘mean, and k the 
magnitude propagation vector. 

g r g £ 

• ^ > -dr < 1.0, i.e., the slopes are relatively small. 

o X dy 

i \ . 

• ^3x^ > = < the roughness is 

isotropic. < > indicates average over an ensemble 
of surfaces. This is not essential to the solution 
but has been employed in developing cross section 
expressions. 

• T >> X , L, i.e., dimension of the illuminated area 

is large compared with the correlation distance L of 
the surface and wavelength X of incident radiation. 

In the small perturbation method, multiple scattering 
and shadowing are not neglected. For the Bass and Bocharov 
method, the third assumption above is not necessary. 

Since the surface C(oc,y) is a random variable, the 
Fourier coefficients associated with it are also random 
variables and may be averaged statistically. Such an 
averaging process yields an average value of for 



small-scale roughness (reference 4): 


and 


^oVV " 4irk^la^^|^ cos^ 9 W(t)'-- 
^oVH " ^oHV " ° 


( 16 ) 

(17) 


where k = 2 it/X , W(t) is the roughness spectral density 

of the surface with t defined as t = 2k sin 6 and with 

a,, and a given by: 
hh vv * 


a 


hh 


[cos 6 + [ 


e - 1 

T~WW 

e - sin d j 


OL 


VV 




For rough terrestrial surfaces, the roughness spectral 

density appearing in equations (16) and (17) can be related 

to the correlation function of the surface through a Fourier 

transform. If the correlation function is Gaussian, i.e., 

2 2 2 

B(t) = Z exp (-T /L ) [see equation (4)], then the radar 
cross sections become: 


.2 „2 .2 I . |2 /-___ 4 , 


2 .2 . 2 , 


OoHh " ^-k^' L‘' (cos 9) exp (-k L sin 9) (18) 
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^oVV ~ ^ i ^ (cos^6) exp (-k^ L^ sin^Q) (19) 

where E and L are the rms height and correlation 
distance o£ the small-scale rough surface. 

Valenzuela (reference 13) used Rice’s theory to obtain 
depolarization from slightly rough surfaces. The depolariza- 
tion is a second order effect. Figures 3-2 and 3-3 have been 
taken from Valenzuela ' s paper and illustrate the depolariza- 
tion predicted using small perturbation theory. A Neumann 
spectrum for the fully developed sea has been used to arrive 
at the results given in figure 3-3. 

The small perturbation method is useful in the low 
frequency limit and can therefore be employed to a class of 
slightly rough surfaces when large radar wavelengths are 
used. A comparison of the theoretical and experimental 
results over a slightly rough surface has been given by 
Wright (reference 14) . As pointed out in his paper, the 
measured compare very well with those predicted by this 

theory. Experimental data over ocean surfaces for angles 
of incidence greater than 20° also shows reasonable agreement 
with the theoretical results predicted using small perturba- 
tion theory. 

Composite scattering theory 

Nearly all natural rough surfaces possess a composite 
structure where small-scale roughness appears superimposed 
on large-scale roughness. Mathematically, the treatment of 
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e ANGLE OF INCIDENCE (DEGREES) 


Figure 3-2. — Polarization dependehce' for backscattering from a 
slightly rough surface for various dielectric constants. 

(Valenzuela, 1967) 



e ANGLE OF INCIDENCE 


Figure 3-3. -Depolarized backscattering cross section per unit area 
for a sea with directional wave spectrum for various windspeeds. 

(Valenzuela, 1967) 
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such a surface is quite formidable, but to a first-order 
approximation one can merely add the cross section calculated 
for the large-scale roughness alone to that caused by small- 
scale roughness. Of course, as 0, the angle of incidence, 
approaches zero, the return from the large-scale structure 
is dominant, especially for surfaces' with gentle slopes. 

Considerably more detail regarding the composite surface 
theory is given by Barrick (reference 4) . In conclusion, 
returns given by equations (18) and (19) should be added to 
the returns given by the large-scale roughness for angles 
of incidence over 20®. 

3.4.2 Radar Return from the Sea - 
A Qualitative Approach 

The calculation of radar return from the sea presents 
a challenge somewhat different but just as thorny as that 
for terrestrial targets. Because of the large number and 
variability of the parameters that produce a particular 
sea state, it is especially difficult to describe that 
state accurately. The windspeed, duration, fetch, and 
direction at the water surface, the ocean currents, con- 
taminants such as oil, the effects of distant storms that 
propagate disturbances of the sea with low loss over vast 
distances, bottom variations, and local weather can all have 
an effect on radar return and are difficult to assess in 
practice . 
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Skolnik (reference 15) gives some of the more useful 

• 't 

terms for describing the sea surface structure: 

Wind wave — a wave resulting from the action of the 
wind on a water surface. While the wind is acting on it, 
it is a sea\ thereafter, it is a swell. 

Gravity wave — a wave whose propagation velocity is 
controlled primarily by gravity. Water waves more than 
5 cm in length are considered gravity waves. 

Capillary wave — a wave whose propagation velocity is 
controlled primarily by the surface tension of the liquid 
in which the wave ds traveling. Water waves less than 
2.5 cm long are capillary waves. 

Fetch — (1) an area of the sea surface over which seas 
are generated by a wind having a constant direction and 
speed; (2) the length of the fetch area, measured in the 
direction of the wind, in which the seas are generated. 

Duration — the length of time the wind blows in 
essentially the same direction over the fetch. 

Swell — ocean waves that have traveled out of their 
fetch. A swell characteristically appears more regular 
for a longer period and has flatter crests than waves 
within their fetch. 

Sea — waves generated or sustained by winds within 
their fetch; opposed to swell. 
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Wave spectrum — the distribution of wave heights (or 
square of the wave .height) with respect to frequency of 
the wave. 

V ' ‘ 

Sea state — the numerical or written description of 
ocean roughness, often referred to as numerical code and 
expressed in terms of the significant wave height. 

Significant wave height — the average height of the 
one-third highest waves of a given wave group. (Height 
is the vertical distance between a crest and a trough.) 

Fully developed sea — the maximum height to which ocean 
waves can be generated by a given wind force blowing over 
sufficient fetch, regardless of duration, as a result of all 
possible wave components being present with their maximum 
amount of spectral energy. 

Sea waves are generated by the wind and differ markedly 
from a swell in physical appearance and in their affect on 
radar return. Individual sea waves are more peaked than 
pure sine waves and are skewed in the direction of propaga- 
tion. They are irregular, chaotic, short-crested (length 
along the crest is of the same order of magnitude as the 
wavelength), mountainous, and unpredictable except in a 
statistical sense. Sea waves contain many small waves 
superimposed on the larger waves, and their spectrum covers 
a wide range of frequencies and directions. 

Swells are more regular than sea waves, longer crested, 
more rounded tops, and more predictable. Their spectrum 
covers a narrow range of frequencies and directions, with 
periods falling between 30 seconds to 5 minutes. 


3-21 



Both wind-generated sea waves and swells can be 
incl’j.ded in the category of gravity waves. Typically, 
the period of gravity waves varies from about 1 to 30 
seconds. Gravity waves are also dispersive, i.e., waves 
of longer wavelength propagate faster than waves of 
shorter wayelength. Ultragravity wave' (also known as 
high frequency gravity waves) have periods of about .1 
second to 1 second. 

« .1 

» 

Capillary waves have periods less than approximately 
0.1 second. Like sea waves, they are generated by the 
wind, but surface tension rather than gravity is the force 
controlling their characteristics. Capillary waves are 
fairly sensitive to the wind. In constrast, if the wind- 
generating gravity waves stop, they continue to run and 
become swells. When capillary waves interact with the 
longer gravity waves, the capillary waves appear to be 
concentrated, at times, on the forward face of the gravity 
wave just before the sharp crest. Capillary waves are 
significant in radar return at the higher microwave fre- 
quencies (X band or greater). 

Wave height is not fixed in relation to the wavelength 
but depends on the wind generating it. Any wave becomes 
unstable and breaks if the angle formed by the crest exceeds 
120°; wave height also can be no greater than one -seventh 
the wavelength. Once the wind is blowing, it takes a finite 
time for a se' to develop. The term fully developed sea 
describes the condition which exists when the ocean waves 
have reached their maximum height generated by a given wind 
force over a given fetch. 
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A typical plot d£ versus 9 , the angle of 

incidence, would appear as in figure 3-4. In the quasi- 
specular region near vertical incidence, the radar echo 
is fairly large with measured values at 9=0 lying 
between 0 and +1-8 dB. As in the case of the "very 
rough” terrestrial targets, this enhanced return near 
vertical is apparently due to specular-like return from 
facet-like areas on the sea surface that are oriented in 
the direction of the radar. 


Above some transition angles (in the order of 20° 
incidence angle) there is little likelihood of significant 
return from the facets making up the sea surface.. Most of 
the return now appears to be due . to sea surface small-scale 
structure such as spray, foam, and capillary waves. Such 
return is relatively isotropig and accounts for the plateau 
region of the o versus 9 plot. At still higher angles 
of incidence falls sharply, but this region is of no 

concern in the S-193 Scatterometer measurements and will 
not be discussed in this report. 


At near normal incidence, measurements (reference 16) 
indicate that decreases with increasing wind on the 

sea. Little difference is seen between ^oVV 

at near normal incidence. In the plateau region 
increases as the wind rises. At low windspeeds 
the plateau region is considerably less than ^ qYy » hut 
as the wind increases, increases faster than ^oVV ’ 

so that with rough- sea conditions there is relatively little 
difference between the two returns . 
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Figure 3-4. — versus 6 for the sea (general observed dependence) 
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3.4.3 Oceanic S-193 Targets - Applicable Theories 

Since the sea state is so sensitive to wind, it is to 
be expected that no single model can predict o for the 
sea. Thus, this section will discuss three models, each 
pertaining to a particular range of sea state, that can be 
used to theoretically determine . These are (1) low 

sea (rms wave height less than X/4tt ) , (2) medium sea 
(h >X/47t and WCt), the roughness spectral density of the 
surface, dependent strongly upon the wind magnitude), and 
(3) rough sea (h >> X/4tt with W(t) weakly dependent upon 
wind) . 

For low seas the radar return contains a coherent and 
a noncoherent return. The noncoherent return is quite small 
and diffused and can be calculated by equations (18) and 
(19) if the sea rms height is known. The coherent component 
has a strong return at 0=0, with side lobes elsewhere, 
and is given by the following equation. (The results of 
Barrick (4) and Kerr (17) suggest this expression.) 

2 

cos^6 exp (-4k^ cos^G) (20) 

where it is assumed that the illuminated area of the sea is 
a circle of area A = ira- , ' is a Bessel function of the 

first kind . . 

For the moderate sea, the composite rough surface model 
will apply; the resultant cross section is the sum of the 
cross sections due to large (h >>X ) and small (E<< X ) 
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scale roughness. Comparison of the composite model theory 
and NASA/ Johnson Space Center 0.4 GHz data is shown in 
figure 3-5 (reference 7). In this figure the total rms slope 
of the large-scale structure is 

S ,= 0.1763 . , 

and the small-scale parameters are 

kE = 0.5 

kL = 20.0 

Near 6=0, the large-scale roughness predominates and 
equations (11), (12), and (13) can be used to predict . 

The required value of S , the total root mean slope, can be 
estimated from the work of Schooley (reference 18) , Cox 
and Munk (reference 19) , or Krishen (reference 7) , who 
have related the mean sea slope to wind. Figure 3-6 has 
been taken from Krishen (reference 7). In this figure the 
rms slope is 3^ in degrees where S = tan 3^ • For large 
values of 9 the return is largely due to small-scale sea 
structure and is given by equations (18) and (19) . Unfor- • 
tunately , W(t) , the roughness spectral density of the 
surface, has no simple functional form (reference 7).. It 
must be evaluated from rather slender experimental evidence 
(reference 20) . 

As wind rises and the sea becomes rough, an upper 

limit exists for the height of a wave of fixed length 

whether in the gravity or capillary range. For portions 

of the ocean where the wave height is limited, the roughness 

- 4 

spectral density W(t) assumes the form W(t) = Bt 

2 2 2 

(reference 21) where t = p + q , P > are the 
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Figure 3-5. — Comparison of NASA/JSC 0.4 GHz data with the 
composite scattering theory. (Krishen, 1971) 


3-2 



Windspeed, knots 


Figure 3-6. - Values of 0^ [s=tan 3^] as a function of windspeed 

(Krishen, 1971) . 




radian wave numbers in x and y directions on the surface. 

The value of B for the capillary waves given in reference 
- 2 

21 is 1.5 X 10 . Thus, for rough seas, a composite surface 

model again suffices, with returns due to the two roughness 
structures contributing to the radar cross section over a 
wide range of angles, including small values of 6 . 

i 

Regarding the matter of a theory that can be used to 
calculate ‘^oHV ’ satisfactory one currently 

exists. Most theoretical treatments of electromagnetic 
wave scattering from a rough surface indicate negligible 
depolarized scatter into the plane of incidence, although 
experimental measurements show that this scatter component 
is present. Some investigators, such as Rouse (reference 
22), have postulated that depolarization is largely a volume 
scattering effect. Rouse's analysis, however, includes 
parameters that are not physically meaningful in an inspec- 
tion of the surface . 
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4.0 ATMOSPHERIC LOSSES 

The computation of the backscattering cross sections 
involves corrections for path losses for transmission and 
reception through the intervening atmosphere. These losses 
appear as constants and in the equation 

[see appendix B, equation B-9). Atmospheric losses are not 
included in the JSC production processed S-193 data where 
Li = Lz = 1. It is therefore important to include (where 
atmospheric losses are significant) proper corrections in 
the o^ computation. In the band (S-193 operates in 
this band) of microwave frequencies, the radar energy is 
absorbed by atmospheric oxygen, water vapor, and rainfall. 
Rainfall can also cause scattering of the microwave energy. 

The loss due to the atmosphere has been studied 

intensively in literature (references 15, 23, 24, 25, 26, 

27, and 28). Experimental data has also been gathered to 

verify theoretical models. In general, the atmospheric 

loss depends on the atmospheric temperature, pressure, and 

water vapor density. Large variations normally occur over 

long intervals and are associated with major changes in 

air mass type at the observing site. LeFande (reference 27) 

divided the atmosphere into a series of 110 spherical shells 

of exponentially increasing thickness to a height of 30 

kilometers. His theoretical results for 60 percent relative 

humidity at 60°F (which gives a ground-level water vapor 

3 

density of 7.5 g/m ) are shown in figure 4-1. The 

attenuation per unit path length increases with frequency. 
For the range of frequencies up to 100 GHz, there is a water 
vapor peak at 22.3 GHz and oxygen absorption peak at 60 GHz. 



Figure 4- 
vapoT 



. — Absorption profiles for various altitudes for water 
concentration =7.5 g/m^ (60 percent RH at 60®F) . 

(R. A. LeFande, 1968). 


The final results of LeFande’s (reference 27) calcula- 
tions giving a total attenuation due to atmospheric oxygen 
and water vapor are shown in figure 4-2. The formulation 

S * ‘ . 

was developed from classical laws on electromagnetic 
attenuation in gases and a hnown oxygen, water vapor, 
pressure, and temperature distribution of the troposhere. 
Experimental data gathered thus far confirms the attenua- 
tions shown in figure 4-2 (see references 24 and 25). 
Another computation by Haroules and Brown (reference 23) 
yields the r suits of figure 4-3 for vertical incidence 
(0 = 0). The agreement between figures 4-2 and 4-3 is 
evident. 


The two-way attenuation can be expressed in terms of 
an exponential law over propagation paths where pressure 
and atmosphere composition are uniformly distributed. The 
attenuation for this case is (reference 15) : 


H 




exp 1-2 


f 


a^dH^ 
a o 


( 21 ) 


where is the attenuation coefficient per meter, Rg is 

the distance between the S-193 antenna and ground resolution 
cell and is the Skylab altitude. For the U*S. Standard 

Atmosphere, the data needed to calculate the value of the 
integral of equation (21) has been given in reference 28. For 
the U.S. Standard Atmosphere, equation (21) can be evaluated as 


0/0364047R 


L^L 2 = exp 


0 


H. 


( 22 ) 


at the frequency of 13.9 GHz. 
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Figure 4-2. — Total absorption versus frequency for paths through 
the atmosphere at various antenna elevation angles (water vapor 
concentration = 7.5 g/m^). (R. A. LeFande, 1968) 
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Figure 4-3. — Total vertical path attenuation by atmospheric gases 
(water vapor concentration = 7.5 gm/m^>. (Haroules and Brown, 1968) 
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Under normal circumstances where no heavy precipitation 
is present, the S-193 radar energy suffers a two-way 
attenuation of approximately 0.16 dB at vertical incidence 
and approximately 0.24 dB at the highest scan, angle of 48®. 
It should be noted that the mission requirements for Sensor 
Performance Evaluation data-takes were such that no data 
was to be taken in the presence of heavy precipitation and 
cloud cover. The data analyzed in this report belongs to 
this category. In the presence of clouds and heavy rain, 
the ’ s will be used with proper recognition given to the 
atmospheric conditions. 

Considerably more attenuation is caused by the presence 
of citjuds and rainfall over the observation site. The rain- 
fall attenuation for frequencies from 4 to 100 GHz has been 
calculated theoretically by Oguchi (reference 29) . His 
results were interpolated by LeFande (reference 27) . These 
results are shown in figure 4-4. No resonances have been 
noted. Stafford (reference 30) computed the attenuation 
constant due to a cloud deck from 6,000 to 12,000 feet with 
moderate- to-heavy rain below. His results are given in 
table 4-1. A review of various experimental and theoretical 
results (reference 31) indicates that attenuation because of 
rainfall cannot be ignored at 13.9 GHz. 


For accurate measurements of from ground scenes, ■ 

it is necessary to consider false returns from heavy 
precipitation. The backscatter from heavy rain or clouds 
can cause more return than the ground, under certain ground 
conditions, and antenna look angle. Typical backscatter 
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Figure 4-4. —Theoretical rainfall attenuation versus 
frequency for various rain rates. 

(R. A. LeFande, 1968) 
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TABLE 4-1. - 13.9 GHz SIGNAL ATTENUATION BY CLOUDS AND RAIN 
(two-way-path attenuation, in dB) 



Rainfall rate, inches/hour 

Angle from 

Drizzle 

Light rain 

Moderate 

Moderate -to- 

Heavy rain 

Very heavy 

nadir (degrees) 

.001 mra/hr 

.04 mm/hr 

rain .16 ram/ hr 

heavy . 39 mm/hr 

.63 mm/hr 

rain 1.6 mm/hr 

0 

2.9 X 10”^ 

1.0 X 10'^ 

6.2 X 10’^ 

2.6 

4.8 

10.2 

12.5 

3.0 X 10'^ 

1.0 X 10'^ 

6.3 X 10"^ 

2.6 

* 4.8 

10.4 

20 

3.2 X 10"^ 

a 

1.1 X 10"^ 

■ 6.7 X 10'^ 

2.8 

5.1 

11.0 

32.1 

3.6 X 10“^ 

1.2 X 10"^ 

7.6 X 10"^ 

3.1 

5.8 

12.5 

35 

3.6 X 10'^ 

1.2 X 10"^ 

7.7 X 10'^ 

3.2 

5.9 

12.7 

43.2 

4.3 X IQ-l 

1.4 X 10'^ 

9.1 X 10"^ 

3.8 

7.0 

15.1 

52 

5.5 X 10'^ 

1.8 X 10'^ 

1.2 

4.8 

8.9 

19.2 















V 


from rain is givei\ in table 4-II (reference 31) •. The 
reflectivity parameter is defined as: 


N 



i = l 


where N is the number of scatterers per unit volume and 
is the droplet diameter. The combined backscattering 
effects of a unit .^volume are > 

N 

i=l 

where a. is the backscatter cross section of the ith 
1 

scatter. 

The reflectivity of uniform rain is expressed in terms 
of radar cross section per unit volume ri(or Eo^) in 
table 4- II. The rainfall rates are given in millimeters 
per hour (mm/hr) . 

The a ’s over rainfall will not be used to determine 
o 

precision/accuracy of S-193 Scatterometer . The uncertainties 
due to the backscatter caused by the rainfall will prevent 
drawing any conclusions on the performance of the sensor. 
Mission requirements were stated so that data would be 
gathered over Sensor Performance Test sites under almost 
clear conditions (less than 50 percent cloud cover) , 
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TABLE 4-II. - REFLECTIVITY OF UNIFORM RAIN, METER^/METER^ (n OR 


* 


dB m ^ 

transmit frequency. 

GHz 

.. ^' R 
CdB) 

Type 

S 

3.0 

C 

5.6 

X 

9.3 

Ku 

15.0 

-3 

Heavy cumulus clouds 
4 gm/m'^ 

-118 

-108 

. 

-98 


14 

Drizzle, 0.25 mm/hr 

-102 

- 91 

-81 

-69 

23 

Light rain, 1 mm/hr 

- 92 

- 81.5 

-72 

-60 

32 

Moderate rain, 4 mm/hr 

- 83 

- 72 

-62 

-50 

41 

Heavy rain, 16 mm/hr - 

- 73 

- 62 

-53 

•43 


*Assumes drop diameter <<A 
**Tropospheric attenuation not included 





5.0 S-193 SCATTEROMETER BACKSCATTER DATA ANALYSIS 

, ■ — • ' - - - ■■ ■ . 

The data which will be used to determine precision/ 
accuracy (P/A) will be reviewed in this section. The range 
of the values of which should be examined for determining 

precision/accuracy will be determined by examining data from 
various SPE ground sites. Selection of data for the determi- 
nation of the sensor performance will be based on the anal- 
ysis of data from typical ground scenes. Parameters which 
influence the data will be outlined. Correlation of back- 
scattering cross section data with the ground location will 
be discussed. 

5.1 LOW CT DATA 
0 

Scatterometer data was collected with antenna scanning 
the deep space in several radiometer/scatterometer operating 
modes. There were three lunar calibration (LC) passes during 
Skylab 2 and 3 missions. Analysis of this data is given in 
reference 32. However, the relationship between the precision/ 
accuracy performance and the deep space data was not discussed 
in this reference. It is this aspect that will be discussed 
in this section. 

Deep space has been assumed to be an excellent "no back- 
scatter" target. Therefore, for Vg and the only inputs 

to the scatterometer receiver were the receiver noise and the 
radiometric temperature of the deep space. The average scat- 
terometer power (Vg’3 and scatterometer noise power (Vj^’) 
when normalized to account for integration time, time con- 
stant, and filter and amplifier gains should be equal. Using 



this argument and the equation: 


(IT). 

<V ' > ^ 


where < > denotes average, (IT)j^, (IT)g denote the integra- 
tion times for noise and signal, , Fg scatterometer 

filter gains for noise and signal, and S ’ °s are the 
scatterometer gains for noise and signal, respectively, 

Martin Marietta Company (MMC) personnel recommended that the 
integration times for. the scatterometer noise be changed in 
the production data processing program (reference 32) ; The 
new integration times (table 5-1) were used in the NASA S-193 
Production Processing Program.", Approximately half of the cor- 
rected signal and noise power minus noise power values were 
negative. This is to be expected since Vg' and are 

randomly varying about the same average values . Computation 
of using the Production Data Processing equations in 

this case would lead to a negative value (see appendix B) which 
is not possible. Therefore, these computations are suppressed. 


When signal plus noise power minus noise power is posi- 
tive, the positive values of P^^/P^ for deep space provide 
an opportunity to examine the variance of for zero 

returned power. To accomplish this, these Pp/Pp values 
were used to calculate’ the backscattering cross sections. 

The program assumes that the sensor is looking at the earth 
(in this case an "earth' which absorbs all the energy and 
returns none) . This hypothetical concept allows the compu- 
tation of the angle of incidence and range needed for evalu- 
ating No computation is done for the case where Pp/Pp 

would be calculated to be negative or zero. The extreme / 
values of and approximate average which fesult from such 

a computation are given in table 5-II. data does not show 

any particular dependence on polarization or roll/pitch angle. 
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TABLE 5-II. - DEEP SPACE S-193 SCATTEROMETER DATA 


Mission/ 
day of year/ 
time GMT 
start/stop 

Mode 

Positive P|^/P-j’ 

Extreme values 
of (dB) 

Approximate 
average of 
a data 

CdB) . 

Maximum 


SL2 

165 

15:42:11.5 

15:44:10.7 

CTC R/S 

-33.67 

-48.06 

-39.5 

SL2 

165 

15:44:34.8 

15:45:32.1 

ITC R/S 

-38.14 

-51.46 

-41.0 

SL2 

165 

15:49:11.8 

15:49:59.6 

ITNC R/S 

-37.83 

-52.47 

-40.5 

SL2 

165 

15:52:0.4 

15:54:41.8 

CTNC R/S 

- 

-33.14 

-51.72 

-38.5 

SL3 

224 

15:55:16 
15:57 :16 

CTNC R/S 

-33.97 

-51.63 

-39.5 

SL3 

224 

15:53:56 

15:54:52 

ITNC R/S 

-39.77 

-48.45 

-41.93 


5 = 4 




^ducirility of 

r,,-j -p^*. r4g TB FCH 



























TABLE 5- II. - DEEP SPACE S-193 SCATTEROMETER DATA (Concluded) 


Mission/ 
day of year/ 
time GMT 
start/stop 

Mode 

Positive Pj^/P.p 

Extreme values 
of (dB) 

Approx iina te 

average of 

a data 
0 

CdB) 

Maximum 

Minimum 

SL3 
2 24 

15:49:35 

15:50:32 

ITC R/S 

-36.82 

-50.90 

-41.9 

SL3 

224 

15 : 47 :10 
15:49:10 

CTC R/S 

-34.07 

-47.90 

-40.33 

SL3 

224 

15:41:16 

15:42:10 

CTC /S 

-34.58 

-48.77 

-41.4 

SL3 

254 

13:58:24 

13:59:17 

CTNC R/S 

-34.04 

-50.43 

-40.2 

SL3 

254 

13 : 53:05 
13:54:01 

CTC R/S 

-34.17 

-41,64 

-39.26 

SL3 
2 54 

13:4^-44 

13:49:57 

ITNC R/S 

-38 . 24 

-50. 07 

-41.8 

SL3 

254 

13:57:06 
13 : 58 : 00 

ITC R/S 

■ 

■ 

-34.00 

■ 

-44.72 

-37.00 
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The data in table 5-II shows that below -33 dB the back- 
scattering cross section precision decreases considerably. 

From the accuracy standpoint, a ground site for which is 

-33 dB and below, the S-193 Scatterometer essentially predicts 
a ”no return” target, The accuracy and precision are influ- 
enced not only by the S-193 Scatterometer but by the analog- 
to-digital (A/D) converter. It is interesting to examine the 
uncertainty caused by A/D converter alone in the very low 
signal range. 

i • 

The uncertainty is caused by a finite voltage range for 
which only one output count number will result (see figure 
5-1). In computing a the noise voltage is subtracted from 
the signal voltage. The difference between signal and noise 
voltages will be the same under each of the following three 
conditions (figure 5-1) : 

• When signal and noise voltage (Vg^^Qj^) and noise 
(VNnom^ are in the middle of the voltage range for 
a particular count, respectively 

• When signal and noise voltage ^Snom 

signal count and noise voltage = ^Nnom 

noise count 

• When signal and noise voltage = ^snom 

signal count and noise voltage +1/2 

noise count 

This will assume importance when signal + noise voltage is 
so low that the difference between signal + noise and noise 
counts are no longer accurately related to the difference 
in actual voltages . The uncertainties caused can be examined 
for a particular case, e.g., when Gg - 1 (highest gain) 
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Figure 5-1. — Effect of A/D conversion on measurement uncertainty 

at low signal voltages. 





and 


^N* 


Under these conditions (see appendix B) 

= constant (Vg - Vj^) (24) 


Now a maximum difference of (VJ ■ V*' . ) will result in 

^ Smax Nmin 

only output . ^ Using equation (24) the 

uncertainty in the measurement of in dB (on the higher 

side) will be: 


A' = 10 log 


10 


V” - V" . 
Smax Nmin 


V' - V" 
Snom Nnom 


(25) 


On the other hand, (VU' . - V'' ) difference causes the same 

V Smin Nmax 


nominal difference 
uncertainty in a (on the lower side) will be: 


■ ''Nnom^ output. 


Hence, the 


B 


= 10 log^Q 


V" - V" 
Smin NmaxJ 


VI4 

Snom 


- V" 
Nnom 


(26) 


> V" 
Nmax 


As a consequence, the real value of 
(a^ + A') to (Oq - B') but if 
only the upper limit (a^ + A’) can be deter- 
mined. Equations (25) and (26) were used to find the range 
of uncertainties for the Ava pass data. At GMT 18:59:16.901 


provided V' . 

^ Smin 

in dB could range from 
V" 


^ . < V" 

Smin ~ Nmax 


the value of 


oVV 


is -37.43 dB. However, due to A/D con- 


version characteristics, the actual 


can range from a 


maximum of -35.85 dB to a minimum of -39.95 dB. At GMT 
18 : 59 : 1 .639 the reported is -29.46 dB for VV polariza- 

tion. The actual a will be within a range of (including 
extremes) -29.19 to -29.75 dB. Computations similar to these 
make it obvious that essentially the data below -33 dB will 
suffer an uncertainty of 2 dB or more because of A/D conver- 
sion alone. This uncertainty reduces to approximately 0.5 


dB at a 
of all 


o 


of -29 dB. It is because of this that the use 


0 ^ data below approximately -33 dB should be avoided 
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Although this section primarily referenced the deep space 
data, all very low values of are being discussed. One 

such example is the 6 = 52.5° (approximately) data over Hurri- 
cane Ava. The spacecraft was in a solar inertial mode. 

Because of this a pitch angle variation was introduced (figure 
5-2). S-193 was operated in CTNC/right-only mode. Thus only 

the 0° .Doppler filters were used. The response curves for 
the 0° Doppler filters are shown in figure 5-3 (reference 33) . 
Computations show ,that beyond approximately 18:59:16 GMT the 
sensor was receiving only a small amount of energy mainly 
through the side lobes of the antenna. (A complete discussion 
of why this happens can be found in reference 33.) The values 
of for the Ava pass, beyond 18:59:16 are given in table 

5-III. The data is quite similar to that of deep space. 
Polarization dependence is lost. For this data the uncer- 
tainty due to A/D conversion alone is approximately 4 dB . 

The variance in the data is due to the poor S/N ratio at the 
input to the antenna. For this data the range of half stan- 
dard deviation is from -36.67 to -32.96 dB, and therefore 
corrections to this data are useless. 

There are many instances where the value suddenly 

drops very low (approximately -37 dB or below). During some 
of these periods, the scatterometer was in a standby mode 
(i.e., transmitter off), thus no scattered signal was 
received, only noise. (Note: the scatterometer was put in 

the standby mode before the radiometer; however, the data 
processing is continued up to the time radiometer was put in 
standby mode.) Two such example^ are: 

• EREP pass 11 14:49:29.415 GMT after this time 

are -37.47 and -42 . 60 dB . 
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TABLE 5 -II I. - a VALUES FOR 6-52.5® EREP AVA PASS 

o 

BEYOND 18:59:16 GMT 


T ime GMT 

a 

0 

Polarization 

18:59:16.9 

-37.43 

VV 


-36.01 

HV 


-30.98 

HH 


-36.50 

VH 

18:59:32.145 

-36.01 

VV 

< 

-37.40 

HV 


-33.03 

HH 


-34.10 

VH 

18:59:47.989 

-37.31 

HV 


-29.56 

HH 


-36.48 

VH 

19:00:3.233 

-36.76 

HV 


-29.98 

HH 


-36.45 

VH 

19:00:17.895 

-35.83 

VV 


-35.81 

1 ■ , HV 


-39.37 

j HH 


-39.14 

VH ‘ 
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PITCH IN DEGREES 



Figure 5-2, — Skylab pitch angle variation as a function 

of time (EREP Ava pass). 


5-11 
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49.800 kHz 


49,900 kHz 


50,000 kHz 


50,100 kHz 


50,200 kHz 


Figure 5-3. — Zero degree filter bank response curves for 
low center frequency (LCF) , middle center frequency 
(MCF) , and high center frequency (HCF) 

Doppler filters. 



• EREP pass 20 14:53:1.45 GMT after this time 

are -43.87 and -50.46 dB> 

Unfortunately, this invalid data has been used in determining 
the average and the standard deviation of (product 

S062-8) . It is therefore recommended that caution be exer- 
cised in using product S062-8 for data segments where 
plunges to low values. 

5.2 OCEAN a DATA 
o 

5.2.1 EREP Pass 5, Gulf of Mexico 

The data over the Gulf of Mexico provides a better oppor- 
tunity to study the precision/accuracy of the S-193 Scatterom- 
eter. This is because of the homogeneity of the electrical 
and surface roughness properties ; This aspect of the data 
will be addressed here. 

Extensive ground truth was collected by NOAA during the 
Skylab data take. The location where ground data was gath- 
ered by the NOAA aircraft is shown in figure 5-4. This data 
was taken nearly simultaneously with the EREP data. To assign 
the data to the proper location on the ground, the S-193 Scat- 
terometer data from GMT 18:2:24.155 to 18 : 4 : 24 : 104 was sorted. 
Thus four sets (I through IV) were located at four ocean loca- 
tions. The locations for VV data are shown in figure 5-4, 
along with the Shylab ground track. Since the 48° pitch angle 
was not achieved, this data is displaced from the data corre- 
sponding to other four-pitch angles. The backscatter data 
corresponding to these locations is given in figures 5-5, 5-6, 
5-7, and 5-8. The ocean surface wind data is given in 
table 5- IV. This data ’was discussed with D. Ross (NOAA/Miami) 
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TABLE 5- IV. - OCEAN SURFACE WINDSPEEDS FOR EREP PASS 5 


(a) Data Set: I, Pass 5 


Date 

Time 

GMT 

Latitude 

Longitude 

Windsreed 
^ measi’red at 
A/C altitude 

Wind at 
20 meters 

. 

73 06 05 

18: 51:01.2 

26.104 

-93.439 

16.0 



18:52:51.0 

26:164 

-93.494 

16.0 



18:55:51.0 

26.263 

-93 . 598 

20.0 



18:57:39.0 

26.318 

-93.664 

17.0 



18:58:51.1 

26.356 

-93.708 

18.0 



19:03:39.1 

26.301 

-93.647 

15.0 



19:06:35.5 

26.290 

-93.598 

20.0 



19:09:39.1 

26.362 

-93.784 

17.0 


. 

19:12:39.1 

26.433 

-93.851 

16.0 



19:15:04.9 

26.488 

-93.950 

17.0 



19:19:40.9 

26.510 

-94.021 

IS . 0 



19:52:24.2 

26.7.57 

-93.878 

2.0 



19:54: 12 . 2 

26.829 

-93.785 

12.0 



19:57 :12. 2 

26.944 

-93.625 

17 . 0 



20:00:12.2 

27.076 

-93.444 

16.0 


Average 



15.8 

14.22 
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TABLE 5-IV. - OCEAN SURFACE WINDSPEEDS FOR EREP PASS 5 (Continued) 


(b)L Data Set: II, Pass 5 


Date 

■' ' - 

Time 

(MI 

Latitude 

Longitude 

Windspeed 
measured at 
A/C altitude 

Wind at 
20 meters 

73 06 05 

18:26:21.6 

25.477 

-92.933 

17.0 



18 : 27 :33. 6 

25..461 

-92.889 

19.0 



18 : 29: 21. 6 

25.439 

-92.823 

12.0 



18:32:21.6 

25.488 

-92.829 

18.0 



18:35:21.6 

25.587 

-92.933 

18.0 



18:37:09.6 

25.642 

-92.904 

18.0 



18:38:57.6 

25.697 

-93 .054 

13.0 



18 : 40: 47. 4 

25.757 

-93.115 

17.0 



18 : 43: 49. 2 

25.856 

-93.214 

16.0 



18 :45:37. 2 

25.917 

-93.268 

17.0 



18:48:01.2 

25.999 

-93.345 

17.0 


Average 

16.54 

14.89 


5-15 


Tllb 
















TABLE 5-IV. -OCEAN SURFACE WINDSPEEDS FOR EREP PASS 5 (Continued) 


(c) Data Set: III, Pass 5 


Date 

Time 

GMT 

Latitude 

Longitude 

Windspeed 
measured at 
A/C altitude 

Wind at 
20 meters 

73 06 05 

18:00:28.1 

24.714 

-92.230 

15.0 



18:01:40.1 

24.758 

-92.269 

16.0 



18:02:52.1 

24.796 

-92.307 

13.0 



18:05:16.1 

24.879 

-92.384 

18.0 



18:06:28.1 

24.922 

-92.422 

15.0 



18:08:17.2 

24.983 

-92.483 

14.0 



18:10:41.9 

25 . 065 

-92.560 

15.0 



18:13:41.9 

25.164 

-92.659 

. 19.0 



18 :15:33. 5 

25.225 

-92.719 

17.0 



18:17:57.5 

25.307 

-92.802 

•17.0 



18:19:45.5 

25.367 

-92.862 

16.0 



18:21:33.5 

25.433 

-92 . 922 

14.0 


Average, 

15.7 5 

14.175 
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TABLE 5 -IV. - OCEAN SURFACE WINDSPEEDS FOR EREP PASS 5 (Concluded) 


(d) Data Set: IV, Pass 5 


Date 

Time 

GMT 

Latitude 

Longitude 

Windspeed 
measured at 
A/C altitude 

Wind at 
20 meters 

73 06 05 

16:20:06.8 

24.774 

91.708 

17 



16:21:22,4 

24.788 

91.791 

17 



16:22:35.0 

24.642 

91.868 

17 



16:23:53.6 

24.565 

91.95 

17 



16:25:09.2 

24.494 

92.021 

17 



16:28:56.0 

24.472 

92 ,098 

13 



16:39:42.8 

24.510 

91 .'917 

4 



16:34:36.2 

24. 565 

91.846 

14 



16:36:29.6 

24.593 

91.978 

7 



16:39:36.9 

24.543 

92.038 

8 



16:49:39.9 

24.554 

92.071 

15 



16:50:55.5 

24.549 

91.994 

12 



16: 52:11.1 

24 . 543 

91.917 

11 



16:54:04.5 

24 . 505 

91,9 5 

12 



16:55:57.9 

24.527 

92.075 

11 



16:59:44.7' 

24 . 505 

92.071 

11 



16:05:08.7 

24.423 

91.989 

10 



16:26:08.8 

24.401 

91.978 

12 



16 :32 ;46 .6 

24,384 

92.016 

12 


; 1 6 :45: 24. 5' 

24.478 

92.131 

13 


Average 

12.5 

11.25 


5-17 











1 . 

ATA SS: I 
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• NOAA AIRCRAFT DATA • 
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Figure 5-4. — Location of f ield-o£-view (FOV) for S-193 Scatterometer W 
data and NOAA aircraft for Skylab 2, EREP pass 5 (Gulf of Mexico). 











BACKSCATTERING CROSS SECTION (cr„) IN dB-> 



0 : 20 30 


ANGLE OF INCIDENCE IN DEGREES’^ 


Figure 5-7. — Backscattering cross section as a function of 6 

for data set III (pass 5). 
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Figure 5-8. — Backscattering cross section as a function of 6 

for data set IV (pass 5) . 
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and the winds were scaled to 20 meters from the ocean surface 
according to his recommendation. The wind direction was not 
uniform for the small resolution cells sensed by the aircraft 
sensors. The weather maps from the U.S. Weather Bureau show 
an average ocean surface wind velocity of 12 knots from 90° 
in this region (table 5-V was compiled by J. Carney/LEC) . 

The data of tables 5-IV and 5-V shows the homogeneity of the 
ocean wind field to within 4 knots (standard deviation). The 
a data also shows’ consistence, in that the dependence on 6 
and polarization is as expected at this wind velocity. The 
values of a and o are equal to within 0.5 dB (rms) . 

The precision/accuracy of this data will be reported in sec- 
tion 6 of this report. 


The homogeneity of the dielectric properties can be 
assessed by analyzing the L-band S-194 Radiometer data. For 
this vertical-looking sensor, the effects of ocean surface 
wind velocity (which varies from 9 to approximately 16 knots 
within one sigma standard deviation) are negligible. The 
measured antenna temperature varies from 105.9 to 105. 0°K 
over a period from GMT 18:03:32.35 to 18:04:36.56. These 
radiometric temperatures have a constant offset because of 
sun angle correction. This correction will reduce all tem- 
peratures by a few degrees. The standard deviation for 184 
data samples is less than 0.3°K. The approximate reflectivity 
R(o) can be calculated using the equation: 




ATM 


where ■ I 2 is the transmittance of the atmosphere, 


■A 


(27) 


is the 


radiometric antenna temperature, 
perature , 


Tg the sky background tem- 


Tg the ground temperature, and tlie radio- 


metric temperature of the atmosphere. For T^ = 4°K, L 2 


0.99, T = 300°K (from table 5-V) , T. = 105°K, T 


g 


ATM 


zs 
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TABLE 5-V. - SURFACE AND ATMDSPHERIC CONDITIfflC FOR SL-2 TARGET SITES 


Site 

number 

Site name 

Skylab 

pass 

number 

Cloud 

cover 

Visibility 

Pressure 

(MB) 

Air 

ten 5 )erature 

(°F) 

Dew 

point 

Average 
winds and 
direction 

Significant 
wave height 
(feet) 

Water 

tenq)erature 

("F) 

750598 

Gulf of 
Mexico 

1 

8/10* 

N/A 

N/A 

82.4 

23 

12 knots 
from 80° 

2 

73.4 

750598 

Gulf of 
Mexico 

5 

6/10* 

N/A 

N/A 

80.6 

24 

12 knots 
from 90° 

2 

75.2 

750598 

Gulf of 
Mexico 

8 

6/10* 

N/A 

N/A 

82.4 

24 

11 knots 
from 90° 

2 

73.4 

750598 

Gulf of 
Mexico 

11 

4/10* 

N/A 

N/A 

82.4 

24 

14 knots 
from 130° 

3 

73.4 

750233 

Great Salt 
Lake Desert 

5 

Clear 

20 miles 

1,022.0 

74 

21 

5 knots 
at 210° 
SSW 

— 



*These cloud covers were reported only for some portion of the ground site. The cloud cover over all the 
ground site area used was less than 50 percent. The duration of the winds was long (fully developed seas) . 



taken 1.5°K for U.S. Standard Atmosphere; the value of 
jR(o) I is equal t‘o 0.812. Furthermore, the variability of 
1 r(o) I for 0.9°K change in is negligible. It is inter- 

esting to compare this value of lR(o) j v/ith those measured 
and quoted in literature. From Von Hippel (reference 34) the 
value o:.‘ dielectric constant is: 

e = 76.7 (1 + j 0.157) 

for water at 25'’C at a frequency of 3 GHz. The value of 
)R(o)| for this case is: 

1r(o) I = 0.80 

This value of 1r(o)| at S-band shows good agreement with 
that predicted from S-194 L-band Radiometer measurement since 
a slight decrease is expected as the frequency increases. 

The S-193 Radiometer data for 6 =; 0.9° varies from 
125.31 to 127. 96°K for GMT from 18:2:24.155 to 18:4:24.104. 
Taking Tg = 10°K and T^ = 297°K, = 0.98, = 4.8°K 

for U.S. Standard Atmosphere lR(o)| = 0.777 for = 127°K. 

It is interesting to compare this with the value of jR(o)j 
calculated from e = 55 (1 + j 0.55) quoted in reference 4 
for sea surface at X-band. For this e, |R] - 0.78. 

This |R(o)t is slightly higher than that calculated 
at 13.9 GHz. This is in accord with the dependence of e 
on frequency. 

No dielectric constant measurements of water with salt 
concentrati ons have been quoted in literature at 13.9 GHz. 
Paris (reference 35) has compiled a table of aqueous sodium 
chloride dielectric constants measured by various investiga- 
tors. At 10 GHz the value of the dielectric constant (ref- 
erence 36) for aqueous sodium chloride at 21°C is: 

e = 65(1 + j 0.44) 
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and for 23.7 GHz at 20®C (reference 37) 

e = 42.97 (1 + j 0.74) 

A linear interpolation between these frequencies yields a 
value of , , 

£ = 59 (1 + j 0.66) 

for the dielectric constant of aqueous sodium chloride at 
20.5°C for a frequency of 13.9 GHz. The reflection coeffi- 
cient using this value of the diele'ctr'ic constant is: 

|R(o) 1 = 0.779 

which is reasonably close to the predicted value of 0.777 using 
S-193 Radiometer data. 

' I 

The following conclusions can be drawn by analyzing the 
S-194 and S-193 Radiometer data: 

• The value of the reflectivity |R(o)1 predicted using 
L-band S-194 Radiometer data shows a reasonably good 
agreement with those calculated at a slightly higher 
frequency (3 GHz) (considering', of course, the fre- 
quency dependence) . 

• The value of |R(o)| calculated using S-193 Radiometer 
data is reasonably in agreement with values calculated 
near this frequency. 

• jR(o) j shews proper frequency dependence behavior 
from S-193 and S-194 Radiometer data. 

• The S-194 Radiometer data shows that the dielectric 
properties of the Gulf of Mexico were homogeneous. 

The S-193 Radiometer also shows homogeneity. The 
interpretation (taking into account the frequency 
dependence of the dielectric properties) of S-194 


REPRODUCIBILITY ; 



and S-193 Radiometer data vtas necessary. The L-band 
frequencies are sensitive to salinity changes (ref- 
erence 38) compared with *the K^-band frequencies which 
are sensitive to the sea state (reference 39). Hence, 
it was deemed proper to analyze data at both frequen- 
cies to arrive at the stated conclusion of reflectivity 
homogeneity of the Gulf of Mexico for the data examined 

5.2.2 ERBP Pass 8, Gulf of Mexico 

The ITNC mode was exercised on June 11, 1973, to gather 
radiometer and scatterometer data over the Gulf of Mexico. 

The location of the FOV for the S-193 scatterometer vertical- 
transmit, vertical -receive data is shown in figure 5-9 for 
GMT 15:21:6.286 to 15:22:50.991. The wind measurements were 
taken with the LTN 51 airboine sensor by NOAA. The aircraft 
was flown at 200 feet altitude and corrections made to cal- 
culate winds at a height of 20 meters. The scatterometer 
data corresponding to the three locations shown in figure 
5-9 is given in figures 5-10, 5-11, and 5-12. The a 
versus 9 figures indicati3 that the general appearance of 
dependence is as expected from theoretical and previous 
experimental data considerations. The polarization dependence 
is also normal. 

The homogeneity of the dielectric properties (or reflec- 
tivity) can be examined by reviewing the S-194 and S-193 
Radiometer data. For GMT 15:21:9.93 to 15:22:59.07 the min- 
imum radi ometric antenna temperature for S-194 is 9Q.3°Kwith 
a maximum oi 92.1°K. This small range of values shows that 
no significant nonhomogeneities were present. The average 
temperature was approximately 91°K. Using equation (27) 
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W-174 


Figure 5-9. — Location of f ield-of -view (FOV) for S-193 Scatterometer W 

W 


data and NOAA aircraft data for Skylab 2, EREP pass 8 (Gulf of Mexico) 






ANGLE OF I 

Figure 5-10. — Backscattering 

foi data set 






BACKSCATTERING CROSS SECTION (o ) IN dB-*- 



Figure 5-11. — Backscattering ci 

for data set 
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EREP PASS 8 S-T93 ITNC 
13.9 GHz 
O VV 

• HH 
AVH 


AVG. LAT: 24.122 

AVG. LONG: -86.65 

DATA SET: II 

AVERAGE WIND VELOCITY: 12.3 KNOTS 






I I t 



BACKSCATTERING CROSS SECTION (a„) IN dB-»- 



Figure 5-12. — Baclcscattering cross section as a function of 6 

for data set III (pass 8) , 
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and table 5-V, the value of reflectivity ClR(o)|) is 0.85i). 
This is in agreement with the values given in section 5.2.1. 
The S-193 Radiometer temperature for 0® incidence shows a 
variance of approximately 2°K with an average value of 
approximately 125®K. The reflectivity computed from equa- 
tion (27) is jR(o)l = 0.768. This is in excellent agreement 
with EREP pass 5 data. 

The wind measurements reported by the U.S. Weather 
Bureau at NASA/JSC (see table 5-V) are quite in agreement 
with NOAA detailed measurements (see table 5 -VI) taken almost 
simultaneously (±1 hour) . The noteworthy aspect of the wind 

t 1 

fields for this pass is that these are very close to -the 
winds present at the time of EREP pass 5 over the Gulf of 
Mexico. This will permit comparisons of this data and pro- 
vide insight into the precision/accuracy performance of the 
S-193 Scatterometer . The ocean wind variations for the data 
shown in figures 5-10, 5-11, and 5-12 are also small (approx- 
imately 5 knots standard deviation). The homogeneity of the 
ocean wind field is very much like that of EREP pass 5. 

The preceding paragraph leads to the conclusion that 
the dielectric constant of the ocean surface is very nearly 
the same as for EREP pass 5 over the Gulf of Mexico. Fur- 
thermore, a fairly high degree of homogeneity was shown for 
the ocean surface for the EREP pass 8. 

5.2.3 EREP Passes 11 and 20, Gulf of Mexico 

Data from passes 11 (SL-2) and 20 (SL-3) offers an 
opportunity to study the performance of the scatterometer 
in intrack contiguous (ITC) mode. During pass 11 the data 
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TABLE 5-VI. - OCEAN SURFACE WINDSPEEDS FOR EREP PASS 8 


(a) Data Set: I, Pass 8 



15:48:10.8 
15:48:48.8 
15:49:22.8 
15:49:58.8 
15:50:36.6 
15:51:48.6 
15:52:24.6 
15:53:00.6 
15:53:36.6 
15:55:24.6 
15:57:12.7 
15:57:48.7 
15:58:24.7 
15:59:00.7 
15:59:36.7 
16:00:50.5 
16:01:26.5 
16:02:02.5 
16:02:38,5 
16:03:14.5 
16:03:50.5 
16:04:26.5 
16:05: 02. 5 
16:05:38.5 
16:08:14.5 
16:06:50.5 
16:0.7: 28.3 
16:08:04.3 
16:08:40.3 
16:11: 06.1 


Latitude 

Longitude 

24 

395 

87. 

001 

24. 

412 

87. 

023 

24 

428 

87 

045 

24. 

45 

87. 

061 

24. 

467 

87. 

083 

24 

505- 

87 

127 

24 

521 

87 

144 

24 

543 

87 

122 

24 

56 

87. 

187 

24 

62 

87 

248 

24 

703 

87 

347 

24 

73 

87 

38 

24 

763 

87 

413 

24 

796 

87 

446 

24 

829 

87 

479 

24 

895 

87 

55 

24 

922 

87 

578 

24 

961 

87 

610 

24 

988 

87 

649 

25 

021 

87 

676 

25 

054 

87 

709 

25 

087 

87 

742 

25 

129 

87 

77 5 

25 

158 

87 

808 

25 

186 

87 

841 

25 

219 

87 

869 

25 

252 

87 

907 

25 

280 

87 

940 

25 

312 

87 

968 

25 

444 

88 

.099 


Windspeed 
measured at 
A/C altitude 











TABLE 5-VI. - OCEAN SURFACE WINDSPEEDS FOR EREP PASS 8 (Continued) 


(a) Data Set : I , Pass 8 


Date 

Time 

GMT 

Latitude 

Longitude 

Windspeed 
measured at 
' A/C altitude 

Wind at 
20 meters 

73 06 11 

16:12:19.9 

25.510 

88.165 

8 



16:12:55.9 

25.543 

88.198 

4 



16:13:31.9 

25.571 

88.231 

. V . 6 



16:14:07.9 

25.604 

88.264 

6 



16:15:19.9 

25.664 

88.330 

7 



16:15:55.9 

25.697 

88.363 

' 7 



16:16:31.9 

25.730 

88.391 

7 



16:17 :43. 9 

25.796 

88.456 

7 



16:18:99.9 

25.829 

88.484 

11 



16:19:31.9 

25.895 

88.550 

6 


Average 

11.9 

10.71 
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TABLE 5-VI. - OCEAN SURFACE WINDSPEEDS FOR EREP PASS 8 (Continued) 

I 

(b) Data Set: II, Pass 8 


Windspeed 
measured at 
A/C altitude 



15:22:55.2 

15:23:31.2 

15:24:07.2 

15:24:43.2 

15: 25: 21. '0 

15:25:57.0 

15:26:33 

15:27:09.0 

15:27:45 

15:29:33 ’ 


15:30:09.0 

15:30:45.0 

15:53:09.0 

15:53:45.0 

15:34:21.0 

15:34:57.0 

15:35:33.0 

15:36:09 

15:36:45 

15:37:21 

15:37:57 

15:38:33 

15:39:09 

15:39:45 

15:40:21 

15:40:57 

15:41:33 

15:42:09 

15:42:45 

15:43:21 


Latitude 

24. 

005 

. 24. 

027 

24. 

044 

24. 

066 

24. 

088 

24. 

104 

24. 

126 

24. 

148 

. 24. 

164 

24. 

170 

24. 

164 

24. 

153 

24. 

126 

24. 

115 

24. 

109 

24. 

088 

24. 

082 

24. 

104 

24. 

109 

24. 

115 

24. 

126 

24. 

131 

24. 

142 

24. 

153 

24. 

159 

24 . 

170 

24. 

186 

24. 

203 

24. 

225 

24. 

241 



86.660 

86.682 

86.699 

86.721 

86.737 

86.759 

86.776 

86.797 

86.797 

86.726 


86.704 

86.682 

86.583 

86.5,61 

86.539 

86.539 

86.517 

86.528 

86.556 

86.583 

86.611 

86.638 

86.666 

86.693 

86.721 

86.748 

86.770 

86.792 

86.814 

86.836 


Average 





















TABLE 5-VI. - OCEAN SURFACE WINDSPEEDS FOR EREP PASS 8 (Concluded) 


(c) Data Set: III, Pass 8 


Date 

Time 

cm 

Latitude 

Longitude 

Windspeed 
measured at 
A/C altitude 

Wind at 
20 meters 

73 06 11 

14:23:23.8 

23.5 

84.633 

19 



14:23:59.8 

23.494 . 

84.694 

20 



14:24:35,8 

23.4B3 

84.749 

13 



14:25:11.8 

23.478 

84.809 

18 



14:25:49.6 

23.472 

84.864 

14 



14:26:25.6 

23.461 

84.919 

10 



14:27:03.4 

23.456 

84.974 

10 



14:28:51.4 

23.439 

85.128 

19 



14:38:27.4 

23.286 

85.853 

21 



14:39:39.4 

23.264 

85.946 

11 



14:40:51.5 

23.242 

86.034 

11 



14:44:27,5 

23.192 

86.287 

16 



14:46:15.5 

23.181 

86.424 

10 



14:48:03.5 

23.187 

86.345 

13 



Average 

15 

13.5 
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was collected in VV and then in VII polarization mode, while 
W polarization was exercised during pass 20. The ground 
track and locations where data was taken for the pitch angles 
of approximately 39.5® and 42° are shown in figure 5-13. 

The ocean surface conditions at the time of pass 11 are 
summarized in table 5-V. Aircraft data with 13.3 GHz Scat- 
terometer was also acquired during this pass. The data from 
pass 11 is shown in figure 5-14.’ The polarization and angle 
dependence of the backscattering cross sections are predic- 
table from mathematical models. Detailed correlations with 
ground and aircraft data will be presented in section 6. 

The ocean parameters for Skylab-3 data-takes over the 
Gulf of Mexico were compiled by J. Carney (LEC/ASD) and are 
given in table 5-VII. This data was collected by the U.S. 
Weather Bureau. The winds did not exceed 18 knots in the 
area for which the data will be analyzed in this repoft. 
Typical data is shown in figure 5-14. It should be noted 
here that in the ITC mode the data is taken over a large 
area for a plot of a versus 0. This makes it possible 
to have relative changes in a from angle- to-angle due 
to entirely different ground locations. Fortunately, for 
passes 11 and 20, the variations in ocean surface wind 
velocity were insignificant. This is also reflected in the 
plots of figure 5-14 since the dependence on 9 is as 
expected from a homogeneous rough target. 

One unique aspect of the data for passes 11 and 20 is 
the land/water interface caused by Yucatan Peninsula (figure 
5-13). At GMT 14; 49 ‘.0.08 of pass 11 the FOV was closest to 
the land/water interface. If one assigns -34.8 dB 
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TABLE 5-VII. - SURFACE AND ATFDSPHERIC OM)ITIONS FCR SL-3 TARGET SITES 


Site 

number 

Site name 

Skylab 

pass 

number 

Cloud 

cover 

Visibility 

Pressure 

(MB) 

Tein)erature 

C“F) 

Dew 

point 

Winds 

Significant 
wave height 
(feet) 

Water 

temperature 

C°F) 

750598 

Gulf of Mexico 

13** 

Overcast* 

[precipi- 

tation) 

N/A 

1,019 

81 

75 

8 knots 
from 130° 

2 

86 

750598 

Gulf of Mexico 

16** 

3/10 

N/A 

1,018 

87 

74 

10 knots 
from 90° 

2 

85 

750598 

Gulf of ^^exico 

20** 

5/10 

N/A 

1,016 

82 

76 

12 knots 
from 130° 

4 

87 

746508 

Tennessee/ 

Indiana/ 

North Carolina 

7 

Clear 

— 

— 

— 

— 

Below 18 
knots 

— 

— 

— 

South Atlantic 

7** 

* 

Less than 
4/10 

— 

' 

— 

— 

Below 18 
knots 

— 


750233 

Great Salt 
Lake Desert 

12 

Scattered 
at 11,000 
feet 

35 miles 

1,014 

83 

47 

4 knots 
from 110° 




*Data used in this report for precision/accuracy analysis did not include areas of overcast/precipitation^ 

**Fully developed seas. 




GULF OF MEXICO 
13.9 GHz 

A EREP PASS n 

AVERAGE OF VV DATA 

FROM GMT 14:47:32 TO 14:48:30 


O EREP PASS 11 

FROM GMT 14:48:52.816 
TO 14:48:54.807 
VH DATA FOR ONE SCAN 

• EREP PASS 20 

FROM GMT 14:49:36.086 
14 * 49*38 24 
VV’OATA FOR ONE SCAN 


ANGLE OF INCIDENCE (e) IN DEGREES- 


Figure 5-14. - Backscattering cross section as a function of 9 

for EREP pass 11. 
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to water at 47° incidence angle (see figure 5-13) and -15.92 
dB to the land surface as measured by S-193 Scatterometer , it 

i • . 

is interesting to calculate the value of 

14:49:0.08. Taking the value of a water-illuminated area of 
60 percent v/ith land the remaining 40 percent of the resolu- 
tion cell, the predicted value of is -19.8 dB. The 

measured value of a 


oVH 


oVH 

at GMT 14:49:0.08 was -21.7 dB . 


For pass 20, the closest point to the land/water interface 


was at GMT 14 : 51 : 0*. 089 . The ocean and land 


oVV 


for 


0 47° were taken, respectively, as -19.09 and -9.79 dB. 

as measured by S-L93 Scatterometer. Assuming one-half beam 
was illuminating the ocean surface (figure 5-13) , the calcu 


lated value of a 


oVV 


is -12.32 dB . The measured value for 


Oovv pass 20, GMT 14:51:0.08 is -13.96 dB. 


The reasonable agreement between the calculated and 
measured values of backscatter ing cross sections at the 
land/water boundary leads to two conclusions: 


• The FOV computation is reasonably accurate as speci- 
fied in TR o24 (reference 40). 

• The response of S-193 Scatterometer to land/water 
interfaces is reasonably accurate. 

For the analysis of precision/accuracy , the ' s over 
the ocean will be used from these two passes. 


For pass 11 from GMT 14:48:7 to 14:49:41, the S-194 Radi 
metric a..cenna temperatures range from 89.1 to 91.9°K. The 
average value is 90.6°K. The reflectivity is fairly constant 
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with an average o£ |R(o)| = 0.835. The S-193 Radiometric 
antenna temperatures also show homogeneity of the surface 
dielectric properties, Skylab pass 20 shows a minimum of 
90.4°K and a maximum of 92.3°K for the S-194 measured radio- 
metric antenna temperature (GMT 14:49:58 to 14:51:42). The 
S-193 Radiometer-acquired data also confirms the uniformity 
of the surface reflectivity. '' 

5.2.4 Other a Data Over the Gulf of Mexico and Atlantic 
0 

5. 2.4.1 EREP pass 13, Gulf of Mexico . A CTC radiometer/ 
scatterometer (pitch and roll offsets equal to zero) mode was 
exercised during EREP pass 13. Figure 5-15 shows the field- 
of-view plot for a data segment. The data to be considered 

in analyzing precision/accuracy will be HH data from GMT 
17:27:01 to 17:29:00 and VV data from GMT 17 : 25 : 36 . 84 to 
17:26:21.3. During these times, negligible precipitation 
was present. The overcast precipitation (table 5-VII) was 
primarily centered over f ield-of -view (FOV) around GMT 17:25:00. 
The data is for low windspeeds (below 10 knots). The theore- 
tical models are most appropriate for comparisons since tan- 
gent plane approximations are valid in this case. The mean 
radiometric temperature over all angles of incidence for S-193 
is 133. 8®K with a standard deviation of 2.4°K (GMT 17:27:00 
to 17:29:00). This is encouraging, since it verifies the 
uniformity of surface wind field as well as the homogeneity 
of surface reflectivity. 

5. 2.4. 2 EREP pass 16, Gulf of Mexico . Two modes were 
exercised during pass 16 over the Gulf of Mexico. The intrack 
noncontiguous (ITNC) radiometer/scatterometer , W mode was 
exercised from GMT 16:04:50 to 16:07:21. Part of this data 
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Figure 5-15. — S-193 Radiometer/Scatterometer 
f ield-o£-view, pass 13. 
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was taken over land. This part was excluded from the analysis 
of backscatter precision/accuracy. The average ocean surface 
winds were approximately 10 knots. The second mode exercised 
was the CTNC left radiometer/scatterometer . Only HH polari- 
zation data was gathered. Figure 5.-16 shows the FOV of the 
sensor. The data averaging was performed on only the ocean 
data; data gathered over land, (last scan) was not considered. 
The homogeneity of the surface properties was once again 
checked (as in section 5 . 2 . 4 . 1) , Similar results were 

.f 

obtained for the surface reflectivity. 

r 

5. 2. 4. 3 EREP pass 7, Atlantic . The CTC scatterometer- 
only mode was exercised over South Atlantic from GMT 14:44:3.6 
to 14:44:24,188. The ocean surface winds were below 18 knots 
(exact velocity not known). The pitch offset was 29.4®. The 
angle of incidence varied from 33.11 to 35.46°. This aspect 
of the data is very important, since the variability in the 
data caused by 6 variation is almost minimum for this mode 
with the pitch offset. Since no significant atmospheric 
effects predominated (cloud cover less than 40 percent) , the 
variance of this data was studied as a possible measure of 
the precision. The variation in angle of incidence (2.35°) 
and polarization mixing (see appendix B) were taken into 
account in the data interpretation. The location of the data 
was around latitude-13 . 9° and longitude -37.0°. 

5.3 LAND a DATA 
0 

The Sensor Performance Evaluation sites were chosen on 
the basis of homogeneity in roughness and surface dielectric 
properties. Great Salt Lake Desert (GSLD) and uniform crop 
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or forest areas were considered prime sites. In this section 
a discussion of the data will be presented. The detailed com- 
parisons with theory and other experimental data will be given 
in section 6.0. 

5.3.1 Great Salt Lake Desert Sqatterometer Data 

5. 3. 1.1 EREP pass 5 . During ,EREP pass 5, a CTC zero- 
pitch, zero-roll offset mode was exercised over the GSLD site. 
Scatterometer data corresponding to horizontal- transmit , hori- 
zontal-receive (HH) and radiometer H-polarization was acquired. 

The coverage for six selected scans (not successive) is 
shown in figure 5-17. The finite period of measurement 
extends the ground coverage. This extension is shown for 
the scatterometer by the dashed area. The center of main 
beam intersection with the ground is shown as ”S” for scat- 
terometer and "R" for the radiometer. There is approximately 
40 percent overlap between the scatterometer- and radiometer- 
sensed area at a particular pitch angle. In general, the 
instantaneous resolution cell on the ground is elliptic in 
shape, but for the CTC mode with zero-degree pitch, zero- 
degree roll offsets, the differences between circles and 
actual shapes are negligible. Because of this, the instan- 
taneous scatterometer FOV has been shown as a circle in fig- 
ure 5-17. The numbers associated with the center point of 
selected cells show the GMT data acquisition time in seconds. 

The FOV for the L-band radiometer operating at a wave- 
length of 21 cm is also shown in figure 5-17. The radiometer 
energy received by the antenna is sampled at a rate of three 
per second. At the nominal altitude of 440 km the antenna 
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\ 

receives energy from approximately a 60-mile radius circle 
on the ground. There is a 97 percent overlap between two 
successive measurements. 

The data for several scans is shown in figures 5-18 

and 5-19. From this data and photographs, it is obvious that 
the GSLD area shown in figure 5-17 is not homogeneous in sur- 
face roughness and dielectric properties. Ground data was 
also gathered over Great Salt Lake. It was a warm day with 
no overcast (table 5-V) . To select a uniform area, the 
aircraft-acquired data and Skylab S-193 and S-194 Radiometer 
data were reviewed. The 13.3 GHz scatterometer data showed 
good uniformity over the flight line shown in figure 5-17 
except over the Wildcat Mountain (see the figure 5-17 area 
marked R-6406B). This aircraft and spacecraft data will be 
compared in section 6. The laser profiler data shows most 
areas are smooth over the flight line. A typical smooth 
surface data is shown in figure 5-20. Power spectral den- 
sities and the autocorrelation function of the surface were 
also computed from the laser data. The results from two 
time segments are shown in figures 5-21 and 5-22. The laser 
data showed data corresponding to two scans from GMT 
17:57:45.875 to 17:57:47.611 gathered over the smoothest 
area. Furthermore, in this area, the Kirchhoff approximation 
could be used, and the correlation distance it) would be 
approximately 8 meters and rms height (h) of 1 meter. The 
values of h and t are typical and are suggested by the 
data shown in figures 5-20, 5-21, and 5-22. Accurate values 
of h and £ for the entire flight line have not been deter- 
mined. Even if these were determined, they would not neces- 
sarily be completely representative of the area covered by 
S-193 Scatterometer for the scans from GMT 17:57:45.875 to 
17:57:47.611. 
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(-roll) (+ROLL) 

-wangle of incidence in DE6REES-+ 


Figure 5-18. — versus 0 for Great Salt Lake Desert, 

pass 5 (I) . 
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BACKSCATTERIN6 CROSS SECTION {a) IN dB- 


I 



(-ROLL) ■ (+ROLL) 

ANGLE OF INCIDENCE IN DEGREES-*' 


Figure 5-19. — versus 9 for Great Salt Lake Desert, 

pass 5(11). 
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Figure 5-20. — Height of the surface in units for a segment of GSLD 

laser profiler data. 
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Figure 5-21 . 


Frequency versus rms 
profiler data 


amplitude for a segment of laser 
over GSLD. 
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Figure 5-22.- — Autocorrelation for a segment of laser profiler 

data over GSLD. 
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The S-193 Radiometer data was pldtted on the FOV scan 
plot (figure 5-23), Radiometer antenna temperatures were 
coded for five ranges. This representation showed unexpec- 
tedly low temperature ranges for the shaded areas shown in 
the figure. This was indicative of surface and subsurface 
moisture. S-194 Radiometer-sensed antenna temperature also 
dipped in this area (figure 5-24), It should be emphasized 
here that the footprint of S-194 is quite large, and the dip 
near GMT 17:57:47 does not necessarily mean that L-band-sensed 
radiometric temperature' is higher than that measured at the 
K^-band. Most significantly, both data in figure 5-24 show 
presence of moisture within the same ground area. Aircraft- 
acquired MFMR X-band data also showed similar results. The 
surface nonhomogeneity in dielectric properties presented a 
serious limitation. Only the data from GMT 17:57:45.875 to 
17:57:49.73 was selected for the comparisons needed to esti- 
mate accuracy. 

The data from the first scan, GMT 17:57:33.14 to 
17:57:34.88, and the last scan (figure 5-17), GMT 17:58:0.702 
to 17:58:2.443 when compared with the remaining data (figures 
5-18 and 5-19), shows that the area within these two scans 
gives a nearly specular return. 

5.3.1. 2 EREP pass 12. A CTC radiometer/scatterometer 
VV mode with pitch offset of 29.4® was exercised during pass 
12. In this mode the 0 variation was only approximately 
2.5® (0 varies from 31.8 to 34.3). Because of this, the 
variance in was considerably reduced. Figure 5-25 shows 

the S-193 Radiometric antenna temperature distribution over 
the sensed area. Once again, the presence of moisture in 
the shaded area was indicated. 
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Figure 5-23. — FOV and S-193 Radiometer antenna temperature 
distribution over GSLD (pass 5) . 
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Figure 5-24. — Comparison of S-194 and S- 

(nearly vertical inci 
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Figure 5-25. — FOV and S-193 Radiometer antenna temperature 
distribution over GSLD (pass 12) . 
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For the area where no moisture was indicated (after the 
fifth scan), the maximum and minimum were -10.78 and 

-7.22, respectively, This was indicative of some degree of 
uniformity. Since the Kirchhoff approximation is not appro- 
priate for angles of incidence greater than approximately 
25° (unless small structure is absent)', no adequate compari- 
sons could be made between theoretical models and this data. 
Furthermore,, this data could not be compared with pass 5 data 
because of different incidence angles and polarization. 

5.3.2 Other Land Sites 

Three land sites covered during Skylab missions 2 and 3 
indicated a fairly high degree of homogeneity. These were 
(1) EREP pass 7, Tennessee/Indiana area; (2) the EREP pass 
13 Colorado/Kansas Area around latitude 37.4° and longitude 
-102.35°; and (3) the area in the vicinity of latitude 37° 
and longitude -97°, EREP pass 15 in Kansas. Although.no 
attempt will be made to model these sites for accuracy deter- 
mination, some insight into the precision can be gained from 
this data. 

5.3. 2.1. Tennessee/Indiana, EREP pass 7 . An ITC, HH- 
polarization mode was operated in the vicinity of 37° latitude 
and -85° longitude. The radiometer data indicated homogeneity 
of the surface. The scatterometer data was highly con- 

sistent. The area had no cloud cover. 

5. 3. 2. 2 Colorado/Kansas, EREP pass 13 . A CTC Radiom- 
eter/Scatterometer mode was exercised with pitch-offset of 
15.6° and roll-offset of 0°. 
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The S-193 Radiometer H-polarization antenna temperature 
averaged approximately 272° with a standard deviation of 
approximately 11° for data from GMT 17:19:50 to 17:20:40. 

This is expected, since the site consisted of dry grass lands, 
sand, and wheat fields and was relatively smooth. 

The scatteroiheter backscatter (a has an average value 

O Xii* 

of -5.95 dB and the standard deviation range is from -8.24 to 
-4.46 dB. It should be noted that at 0 - 17° scatterometer 
backscatter is relatively less sensitive to the surface rough- 
ness than at other angles. The averages and standard devia- 
tion for radiometer and scatterometer data were taken over 
260 data samples. 


5. 3. 2.3 Kansas , EREP pass 15 . A CTC R/S, HH-polari- 
zation mode with roll-offset 0° and pitch-offset 29.4° was 
exercised in the southeastern corner of Kansas. This area 
is mostly smooth With a fairly small area covered by surface 
water. Other than being used as pasture, parts of the area 
are used to grow wheat and milo. The S-193 Radiometer meas- 
ured average antenna temperature from GMT 16:37:39 to 
16:38:20 of 280°K with a standard deviation of 5.7°K, The 

a had a total variation from a minimum of -12.35 dB to a 

oHH 

maximum of -7.11 dB . The average value of a is -9.099 

oHH 

and the standard deviation range is from -10.17 to -8.24 dB. 
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6-0 PRECIS IQN/ACCURACY ESTIMATES 

The parameters which influence the backscattering cross 
sections belong to one or more of the following categories: 
sensor, intervening medium, and the ground scene. Consid- 
erable effort was expended to select data over uniform 
scenes with minimum atmospheric attenuation (clear skies or 
less than 50 percent cloud cover with no rain) . This made 
it possible to study the variation of caused predomi- 

nantely by the S-193 Scatterometer system. Obviously it is 
very difficult, if not impossible, to find a completely 
uniform site (including the intervening medium) and, there- 
fore, the values of precision/accuracy are the worst case 
estimates (or pessimistic upper bounds) . 

6.1 PRECISION ESTIMATES 

Precision is given in terms of one standard deviation 
computed from data for a given mode and polarization 
(assuming a homogeneous ground scene, including intervening 
medium). The precision estimates are given in table 6-1. 
From the (dB) data the mean value was computed by con- 

verting the dB values into numbers, averaging, and finally 
converting back to dB values. For the standard deviation, 
a^'s (not in dB) were used in the following formula. 

p (standard deviation) = 

where n is number of data values used. The standard devi- 
ation range was computed from [mean ~ p] The 
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precision was expressed as the larger of the two values of 
{[mean ± p] dB minus [mean a^] dB}. The range of the 
values of measurements is also given in table 6-1. 

In the interpretation of the results given in this table, 
the following remarks should be kept in view: 

• For the EREP passes 5, 8, and 11 over the Gulf of 
Mexico, only 17° angle of incidence data was used. 

Out of the five angles at which data was collected, 

^ the 17° angle shows the least variation with ocean 

surface winds/sea state. This has been established 
theoretically and experimentally. The dependence 
of a on polarization is as expected. 

• For the CTC, pitch offset = 29.4°, roll offset = 0°, 
mode there are variations due to two factors which 
contribute to the data variance. First, the angle 
of incidence varies by about 2.4° within a scan. 

At an angle of incidence' of 33° this could cause a 
variation up to 2 dB for ocean winds up to 15 knots 
(pass 7, GMT 14:44:3.6 to 14:44:24.188) (refer- 
ences 7, 41, and 42), Second, the antenna motion is 
such that the polarization states for the received 
and transmitted signals are not horizontal or 
vertical but in between. This effect, however, is 
small for this mode since the maximum angle of cross- 
track motion is 11°, and the backscattering cross 
sections for W and HH polarizations are approximately 
the same for moderate ocean windspeeds. Despite 
these variations because of angle of incidence and 
polarization, the data shows small variation 

(table 6-1, EREP passes 7 and 15). This is indeed 
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TABLE 6-1. - PRECISION ESTIMATES FOR SKYLAB 2 AND 3 SCATTEROMETER 


EREP pass/day 
of year (DOY) 



18:07:10 Gulf of 
Mexico 


Mode/polar- 

ization 


Tmc/w 


mc/tti 


mc/vH 


Oq DATA 


Average Nimber Mean deviation range 

angles of .of '^o 

incidence degrees samples (dB) 


0.98 0.67 1.27 


1.36 1.12 1.S9 


-14.9 -15.37-14.47 



Variation Precision 
CdB) I better than 
(dB) 


Maxinun 


0.47 1,58 0.31 


0.85 1.76 0.27 


-15,84-14.43 0.47 



14:44:36 14:44: South CTC/W 

24.188 Atlantic 

P=29.’4) 


CTCA« 

(R=0, 

P*29.4) 


14:27: Teai^see ITC/ffi 

6.821 Indiana 

North 
Carolina 



8/162 

1 

15:20:45 

11/165 

14:47:42 


14:48:30 



riNC/W 

17.7 

13 

-0.36 

-1.86 

0.74 

-3.11 

1.34 

ITNC/m 

17.8 

13 

-0.04 

-1.47 

1,03 

-2.96 

1.43 

rrac/vH 

17.7 

13 

-15.44 

-16.78 

SBl 

-18.5 

-13.94 




-14,02 -14.34-13.73 


-2.56 2.7 2.0 


-14.6 -13.58 0.32 


"The standard deviation is large and therefore [mean a - p] is either negative or yields a value vMch is even smaller than the minnnum 
tf in the data set. “ 









































































































































TABLE 6-1. - PRECISION ESTIMATES FOR SKYLAB 2 AND 3 SCATTEROMETER (Continued) 

a DATA 
o 


EREP pass/day 
of year (DOT) 

Off 

Area 

Mode/polar- 

izaticm 

Average 
angles of 
incidence degrees 

Nunber 

of 

sanfiles 

Mean 

Oo 

(dB) 

a Standard 
0 

deviation range 
(dB) 

°0 

Variation 

(dB) 

Precision 
better than 
(dB) 

Ftob 

To 

Minium 

Maximal 

Minimm 

Maximal 

13/216 

17:25: 

36.84 

17:26: 

21.3 

Gulf of 
Mexico 

CIC/W 

(P=0“, 

R?0<*) 

J 

Minus 

roll 

11.0 

21 

5.6 

2.92 

7.25 

0.65 

7.77 

2.68 

9.26 

21 

8.41 

7.79 

8.95 

5.39 

9.88 

0,62 

7.53 

21 

11.1 

10.44 

11.66 

9.15 

11.96 

0.66 

5.77 

21 

12.80 

12.32 

13.23 

11.78 

13.56 

0.48 

4.01 

21 

14.17 

13.52 

14.73 

13,3 

15.11 

0.65 

2.29 

21 

14.98 

14.33 

15.59 

13.97 

15.93 

0.65 

Plus 

roll 

1.2 

21 

-15.19 

14.64 

15.67 

14.62 

16.29 

0.55 

2.14 

21 

14.62 

14.26 

14.95 

14,03 

15.38 

0.36 

3.79 

21 

13.44 

13.08 

13.77 

12.82 

14.07 

0.36 

5.53 

21 

12.07 

11.64 

12.46 

10.85 

12.56 

0.43 

7.27 

21 

10.37 

10.14 

10.59 

9.98 

10.73 

0.23 

9.02 

21 

8.38 

8.22 

8.54 

8.08 

8.65 

0,16 

13/216 

17:19:50 

17:20:46 

Colorado/ 

Kansas 

CTC/HH 

(P=1S.6® 

R=0") 

17.5 

264 

-5.95 

-8.24 

-4.46 

... 

-9.53 

-2.05 

1.49 

15/217 

16:37:39 

16:38:20 

Ksnisas 

CTC/HJ 

(P=29.4'’, 

RfO") 

From 

32.69 

to 

34.07 

223 

D 

-10.17 

-8.24 

-12.35 
^ 

-7.11 

1.07 

16/220 

16:6: 

24:006 

16:7: 

21:179 

Gulf of 
Mexico 

ITNC/R/S 

W 

49.59 

14 

-22.22 

-23.56 

-22.0 

-24.8 

-21.0 

1.34 

43.0 

16 

-20.1 

-21.42 

-19.07 

-21.65 

-18.95 

1.41 

31.41 

16 

-14.0 

-14.46 

-13,61 

-14.65 

-13.4 

0.46 

16.44 

16 

1.5 

1.33 

1.73 

1.15 

1.87 

0.23 

0.85 

12 

13.8 

13.54 

14.06 

13,26 

14.13 

0.26 
























































































TABLE 6-1. - PRECISION ESTIMATES FOR SKYLAB 2 AND 3 SCATTEROMETER (Concluded) 

0 DATA 
0 


EREP pass/day 
of year (DOY) 

(MT 

Area 

Mode/polar- 

ization 

— 

Average 
angles of 
incidence degrees 

Number 

of 

saii5>les 

Mean 

CdB) 

a Standard 
0 

deviation range 
(dB) 

a 

0 

Variation 

(dB) 

Precision 
better than 
(dB) 

From 

To 

Minimum 

Maximum 

Minimun 

MaxijTiuin 

16/220 

16:07:39 

16:09:16 

Gulf of 
Mexico 

CTNC 

L(R/S) 

HH 

49 

20 

-29.74 

-31.15 

-28.68 

-32.7 

-28.45 

1.41 

41,4 

24 

-23.36 

-26.09 

-21.74 

-28,75 

-19.97 

2.73 

30.4 

20 

-17.47 

-19.66 

-16.02 

-21.91 

-15.23 

2.19 

15.94 

20 

0.87 

0.4 

1.3 

0.08 

1.53 

0.43 

0.28 

16 

13.96 

13.66 

14.24 

13.63 

14.59 

0.3 




encouraging. All the comments made thus far also 
apply to EREP pass 13, GMT 17:19:50 to 17:20:46 when 
a CTC, pitch offset = 15.6®, roll offset = 0® was 
exercised. 

• For the ocean data, the windspeeds remained 

fairly constant within the areds viewed by the antenna 
during the times used in table '6-1. However, no 
accurate measurements of wind direction were avail- 
able. is a function of the wind direction. 

Figure 6-1 shows data taken with 13.9 radiometer/ 
scatterometer aircraft system. This data was col- 
lected by NASA/Langley Research Center. Though the 
variations in wind direction during Skylab passes 
over the ocean were relatively small, it is clear 
that up to 1 dB variation in could have resulted 

from wind direction alone. 

• The data used in computing the precision estimates 

in table 6-1 was carefully correlated with the 
ground scene. There were several land/water bound- 
aries causing sudden changes in data. In these 

instances, the data was sorted according to land or 
water site. All data within the time intervals shown 
in table 6-1 was used except for pass 16, day of the 
year (DOY) 220 GMT 16:07:39 to 16:09:16. From this 
data, two sets of four data values were dropped 

corresponding to GMT 16:8:10.072 to 16:8:12.525 for 
the highest angle, and GMT 16:9:8.427 to 16:9:9.137 
for the lowest angle. These data were abnormally 
different from the remaining data. The lowest angle 
data was dropped because the sensor FOV was on land. 
However, no particular reason was determined for the 
behavior of the highest angle data. 
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BACKSCATTERING CROSS SECTION o 



Figure 6-1. — Radar cross section versus relative wind-heading (NASA/LaRC) . 





• One aspect of table 6-1 warrants special explanation. 
It is the angle of incidence which has been used in 
sorting the data. The pitch and roll girabal angles 
are used to determine the angle of incidence. The 
gimbal angles did not stay constant during data 
runs. In modes which included 48° command angles, 
the flex harness stiffness prevented the antenna 
from reaching this angle. The, attained angle ranged 
from 43 to 46°. The analysis of gimbal angles is 
given in referepce 43. It appears that variations 

of the pitch and roll gimbal angles ranged up to 
0.5° for SL-2 and SL-3 during a pass. The variations 
for individual data takes analyzed in table 6-1 are 
smaller. However, any variation in these angles 
causes a corresponding variation in 9. Thus, in 
the presence of a homogeneous scene, the will 

show some variation since it is dependent on the 
angle of incidence. This is therefore one reason 
for variance of caused by the system. 

• The variance of the values given in table 6-1 

is the sum of variances caused by the sensor and 
the ground scene, since these two effects are inde- 
pendent. Because of this the precision estimates 
given in this table should be considered upper 

bounds. Each scatterometer is a result of 

o 

several independent measurements (see appendix A) . 
Therefore, the precision of is expected to be 

good. 

The values of the precision upper bound from table 6-1 
have been plotted in figure 6-2. Out of 39 values only 
five are greater than of equal to 2 dB, Por the remaining 
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PRECISION UPPER BOUND (dB) 



AVERAGE BACKSCATTERING CROSS SECTION (ao) IN dB ► 


Figure 6-2. — Precision upper bound as a function of (cr^) from table 6-1. 



34 sets of data, the precision upper bound is 1.5 dB. The 
trend for these 34 data points (figure 6-2) is such that 
precision upper bound is lower for greater than 0 dB 

and higher for less than 0 dB. This is expected since 

the signal- to-noise ratio deteriorates for low a^. From 
the data analayzed, it can be concluded that the precision 
upper bound is 1.5 dB for S-193 Scatterometer for a 
range of 18 to -30 dB. 

I 

6.2 ACCURACY ESTIMATES 

The remarks made in section 6.1 are also true for the 
accuracy estimates. Accuracy will be infaFied by two methods. 
The first is the . ’iparison of S-193 data with other 

experimental data; the second compares theoretical values 
of calculated using ground truth with S-193 data. 

6.2.1 Comparison of S-193 Data with 

Other Experimental Data 

During the period of Skylab missions 2 and 3, data was 
gathered over selected sites with 13.9 GHz radiometer/ 
scatterometer aircraft-borne system;. In figures 6-3 and 
6-4 the comparison of data for nearly the same ocean surface 
winds is given. The data was taken within 1 hour of the 
Skylab overpass. The data shows excellent consistency. 

The two types of data (Skylab and aircraft) show a differ- 
ence of approximately 2 dB. It should be noted that 13.9 GHz 
radiometer/scatterometer aircraft Mission 247 data was 
gathered over an area part of which was under cloud cover 
and moderate shower activity. The aircraft - scatterometer 
data was processed at NASA/Langley Research Center. 
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13.9 GHz AIRCRAFT 
RAD/ SCAT 

HH, MISSION 238, 
FIT 20 6/5/73 
UPWIND 13 KNOTS 


13.9 GHz S-193 

SCATTEROMlTER 

HH, 

PASS 5, 

DATA SET. II 
WIND 15 KNOTS 
FROM 90° 

HH, PASS 5 
DATA SET I 
WIND 14.22 KNOTS 
FROM 90* 



ANGLE OF INCIDENCE IN DEGREES 


Figure 6-3. — Comparison of 13.9 GHz radiometer/scatterometer 
aircraft data and S-193 Scatterometer data (HH polarization). 
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i THfi 

T-n TiAAVa 


BACKSCATTERING CROSS SECTION IN dB 



0 5 10 15 20 25 30 


ANGLE OF INCIDENCE IN DEGREES — ► 

Figure 6-4. — Comparison o£ 13.9 GHz radiometer/scatterometer 
aircraft data and S-193 Scatterometer data (VV polarization) . 
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The second comparison was made with S-193 altimeter- 
acquired data. The S-193 altimeter is a high signal- 

^o-noise system (18 to 30 dB, depending on mode) . It was 
therefore felt that the performance of S-193 Scatterometer 
could also be verified by a comparison such as is given in 
figure 6-5. Details lof the computations for the altimeter 
a are given in reference 44. These data were not taken 
simultaneously, but over nearly the same ocean windspeeds . 
Only a limited number of values were available for the 
calculated from S-193 altimeter data (through NASA/Wallops 

I 

Space Center) . The S-193 Scatterometer data is close to 
the altimeter (at the most a difference of 2.8 dB) . 

It is not intended here to provide a number for the accuracy 
using figure 6-5, but rather to verify general agreement 
between the two kinds of data gathered under nearly similar 
ocean surface conditions. 

For several years, a number of aircraft missions have 
been flown by NASA/JSC to study the dependence of radar 
return on such parameters as local windspeed, wind direction, 
and the spectrum of the sea. In figure 6-6, the comparison 
of the data from one of these missions (13.3 GHz NASA/ JSC 
aircraft mission 119) and S-193 Scatterometer is presented. 
Mission 119 was conducted over the North Atlantic in 1970. 
Once again the three data sets are within 2.5 dB. The maxi- 
mum difference between the aircraft and Skylab data is from 
0 to 5° angles of incidence. ' 

Nearly simultaneous data was obtained with 13.3 GHz 
scatterometer over the Gulf of Mexico during EREP pass 11. 
When the two sets of data (from aircraft 13.3 GHz scatterom- 
eter and S-193 Scatterometer) were compared, a 15 dB constant 
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Figure 6-5. — Comparison of altimeter mode 2 a data and S-193 

Scatterometer data. ° 
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BACKSCATTERING CROSS SECTION (Oo) IN dB 




13.9 GHz S193 
PASS 11 

• 

AVERAGE OF 13 CYCLES 
FROM 14:47:36.90 
TO 48:26.868 


WIND SPEED 14 KNOTS 
POLARIZATION VV 


13.3 GHz SCAT 

□ 

MISSION NASA/ JSC 

119 FLIGHT 5 LINE 3 RUN 1 

FOREBEAM VV 


WIND SPEED 11.0 KNOTS 


13.9 GHz S193 

A 

PASS 5 
DATA SET IV 

WIND SPEED 11.25 KNOTS 
VV POLARIZATION 


□ 

A • □ 


f . g 


15 20 25 

ANGLE OF INCIDENCE IN DEGREES 

Comparison of 13.3 GHz aircraft and S-193-acquired 
Scatterometer data. 



difference was noted. It has been suggested that the 15 dB 
higher 13.3 GHz data was due to new system calibration, 
which was not perhaps reflected in the data processing. 
Further evaluation could not be pursued due to the lack of 
time. However, it was. felt that a comparison of -15 dB) 
from 13.3 GHz with S-193 data could verify the general 

dependence on the angle of incidence. The resulting compari- 
son is shown in figure 6-7. The dependence of S-193 

data on the angle of incidence compares very well with that 
of the aircraft- acquired data. 

Additional data sources were not found in the literature 
at or close to 13.9 GHz frequency. This is why no further 
comparisons could be made. 

6.2.2 Comparison of S-193 o^ Data with 
Theoretical Results 

In section 3.4 the applicable theories to the scattering 
from ocean surfaces were discussed. The composite surface 
model was selected here for the comparison with the experi- 
mental data. The data from three Skylab passes is 

shown in figure 6-8. The average wind velocity for these 
data is approximately 14 knots. The total mean square slope 
(S ) was calculated for waves longer than approximately 

o 

1 foot from the following equation (reference 45) : 

si = 0.01 (0.8 + 0.08V) (28) 

o 

where V is ocean surface wind velocity in knots. 
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BACKSCATTERING CROSS SECTION (go) IN dB 



0 10 ?.0 30 40 50 60 70 


ANGLE OF INCIDENCE IN DEGREES ► 


Figure 6-7. — Comparison of aircraft-acquired data (Oq -15 dB) 
at 13.3 GHz and Skylab 13.9 GHz scatterometer data. 


6-17 




20 30 40 50 6 


ANGLE OF INCIDENCE IN DEGREES 

Figure 6-8. — Comparison of theoretical results using composite 

model and S-193 data. 



2 

The total mean square slope (S ) for all waves (high 
frequency and gravity waves) was computed from the following 
equation (reference 45) : 

= 0.01 (0.3 + 0.264V) (29) 

The high frequency waves (includes capillary and high 

frequency gravity waves) are sensitive to wind velocities 

much more than the low frequency gravity waves (reference 46) . 

The total root-mean-square slope for these high frequency 

2 

waves was computed by subtracting the values of S from 

2 ® 

S . For the 14-knot wind velocity, this yielded the total 

2 

mean- square slope (S,) of 0.0192. Wu has given the following 

^ 2 

relationship between and the root mean square height 

(E) (reference 46) , 

= 0.01 (0.4 In Z + 3.38) (30) 


Equation (30) is valid up to the ocean wind velocity 
of approximately 7 meters per second. In this report it 
will be used up to 7.2 m/sec. For the Gaussian distributed 
high frequency waves, the correlation length (L) can be 
calculated from 




(31) 
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For the 14-knot ocean surface vrind, equations C-^) 
through (31) yield the following values; 

= 0.04 
kS = 0.087 
kL = 1.233 


For e = 55 + j30 and R = 0.78, the backscattering 
cross sections were computed from the following equation 



where is the distance from S-193 antenna to the illumi- 

nated area and the Skylab altitude. 


The first factor on the right side gives the atmospheric 

losses for the clear atmosphere. C<^qYy) is the backscatter 

h 

from equation (11) and from equation (19) for 

large and small-scale ocean surface roughness, respectively. 
The results are shown in figure 6-8. The theoretical and 
experimental values are within 3 dB. This result is quite 
encouraging. Comparison of theoretical for HH polariza- 


tion combination and S-193 a 


oHH 


data showed agreement to 


within 3 dB for a 


range of 14 to -28 dB. 


The comparison of theoretical results and S-193 a 

2 ^ 

data is based on the relationship between S^ (mean square 
slope) and Z (the root mean square height of the high 
frequency waves) given by Wu (reference 46). However, 
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experimental verification of the dependence of E. on the 
ocean surface wind velocity has not been reported. The 
theoretical results should therefore be considered only 
approximate. It is interesting to note that the best com- 
puter fit to the data set I of piss 5 (VV polarization) was 
obtained for the following roughness parameters: 

= 0.032 
kZ = 0.11 
kL = 2.25 

A comparison of the theoretical and experimental results 
is given in table 6- II. 

TABLE 6- II. - COMPARISON OF THEORETICAL AND EXPERIMENTAL 
Oq data for best COMPUTER FIT 


0 

EREP pass 5 
polarization 
VV (dB) 

Theoretical 
results (dB) 

0.919 

13.21 

12.52 

17.1 

1.42 

0.85 

32.17 

-11.9 

-12.76 

43.7 

-16.22 

-16.3 

50.1'y 

-17.65 

-18.56 


The best computer fit is within 1 dB of the value of 
the S-193 data. The set of roughness parameters is 

also quite reasonable. From this it would seem that the 
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capillary and small -gravity structure has a larger correla- 
tion length than what is predicted from Wu's (reference 46) 
analysis. It should be noted that the direction of wind 
relative to the Skylab ground track also influences the 
backscattering cross section (see figuie 6-1) . The effect 
of wind direction was not included in Wu's experiments and 
has not been reflected in the accuracy analysis given in 
this report. 

The gently undulating profile of the Great Salt Lake 
Desert satisfies the requirements of the tangent plane 
approximation. The backscattering cross sections can 
therefore be calculated from equation (11) of section 3.4. 
The comparison of theoretically computed values of for 

= 0.08, e = 3.4 + jO.6 and S-193 data is shown in fig- 
ure 6-9. This dielectric constant is typical of sandy 
surfaces (reference 47) . The difference between the two 
sets of data '^hown in figure 6-9 is less than 3 dB . Also 
shown on the same figure is the aircraft- acquired 13.3 GHz 
data. The three sets of values are very consistent 

regarding the dependence on the angle of incidence. 


From the analysis and comparisons presented in this 
section, it can be concluded that the accuracy of is 

better than 3 dB. Furthermore, this estimate of accuracy 
is expected to be valid for from -28 to 18 dB. The 

signal- to-noise ratio increases for higher cTq's- This is 
why cy^'s up to the saturation limit of 18 dB are included. 
It is not possible to further tighten the limits of this 
estimate by theoretical modeling, since the ground truth 
required for such a calculation is not available (correla- 
tion distance, root mean square slope for high frequency 
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O 13.9 GHz SCATTEROMETER 
GREAT SALT LAKE 
PASS 5, VV 
GMT 17:57:48.946 TO 
17:57:49.715 
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V 13.3 GHz 
A/C DATA 
VV POLARIZATION 
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0 2 4 6 8 10 12 


ANGLE OF INCIDENCE 0 IN DEGREES 


Figure 6-9. — Comparison of S-193 Great Salt Lake data with 
theoretical results and 13.3 GHz data. 
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capillary and gravity-waves, and exact atmospheric attenu- 
ation) . In view of the limitations imposed on the mathe- 
matical models, it is felt that the accuracy of S-193 
data is certainly better than the upper bound of 3 dB. 

6.3 COMPARISON OF S-193 SCATTEROMETER PREFLIGHT 
PERFORMANCE AND PRECISION/ACCURACY 
ESTIMATES 

An error analysis of the scatterometer system was per- 
formed by Moore (reference 48) using measured preflight 
sensor parameters. The following assumptions are involved 
in his analysis: 

• There are two basic types of errors: bias or 

systematic errors, both known and unknown, and 
random errors. 

• Known bias errors can be calibrated out of the 
system with data processing. 

• Unknown bias errors, such as switch insertion loss 
and short term temperature variations, etc., are 
fixed during any one measurement sequence but change 
randomly from one measurement sequence to the next. 

• Random and unknown bias errors are independent and 
will be root-sum-squared (rss) to provide a most 
probable measure of their total effect. 

• Those components which are not temperature-controlled 
will be monitored to determine their temperature so 
that the error introduced will be of the known bias 
type. 
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• Losses before any amplifier in a chain will be lumped 
together. 

• Rotary joint insertion loss variation due to position 
occurs at a. rate slower than the radiometer chopping 
rate and will be rejected. For the scatterometer 
calibration, it will be considered an unknown bias 
error. 

• Intermediate frequency bandpass filters will be 
assumed to be ideal rectangular bandpass filters 
and bandwidth as specified. They will otherwise be 
subject to drift and variation in their other param- 
eters such as insertion loss, bandwidth, etc. 

• RF unknown bias errors will be rss together to com- 
pute a single unknown bias error. 


The following error sources have been identified by 
Moore (reference 48) . 

• Transmitted power error 

The mean loss in the transmit path is 1.323 dB known 
bias with an unknown bias error of ±0.059 dB . This 
variational error is attributed to the circulator 
unit whose long term drift is slow but will not be 
able to be calibrated inflight and the unknown bias 
error due to switch indexing. The rotary joint 
error is also assumed to be cyclic in nature, i.e., 
a bias error whose peak level is known but there 
belag no way to calibrate its affect except statis- 
tically by data processing. For these reasons, 


PASS w 
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the variations in the transmit path have been rss 
together and considered an unknown bias error with 
random distribution. 

• Calibration path loss from transmitter to circulator D 
(see appendix A) 

The attenuation in the path is 105.87 dB with an 
unknown bias error of ±0.028 dB. 

• Transmitter leakage 

Leakage from the transmitter may enter the scatterom- 
eter receiver through circulator D. The total isola- 
tion is 120 dB minimum with an unknown bias error 
of ±10 dB. 

• Rotary joint leakage 

The rotary joints have a maximum radiation leakage 
of 95 dB below the power level transmitted through 
them. Assuming the leakage to eminate from the 
most undesirable spatial point and the leakage 
opening to act as an isotropic radiator, the free- 
space loss at 13.9 GHz from the rotary joint will 
be 45.02 dB/foot. The rotary joints are about 
1 foot behind the dish edge. It will be assumed 
that the dish will also act as an isotropic antenna 
in the near field behind the dish. Therefore, the 
leakage signal due to two rotary joints is 95 dB 
+45 dB -3 dB (two joints) equal to 137 dB below the 
transmitted power through the joint. 

• Coupling between channels 

The dual -channel rotary joints leave 95 dB maximum 
leakage between channels. The output of the scatte- 
rometer transmitter couples directly into the 
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down converter input as seen in figure A- 5. Since 
there are two dual-channel rotary joints, there are 
two paths of leakage, i.e., -95 dB/path +3 dB/two 
joints = -92 dB leakage into the down converter. 

The expected power level of the scatterometer cali- 
bration signal at the output of the TDA is -105,87 dB 

before TDA +30 dB TDA gain = -75.87 dB. 

• Loss from circulator E to TDA input (see figure A- 5) 

The loss is 6.369 dB with an unknown bias error of 
±0.058 dB. 

• Loss/gain of down converter i.f. ampere power splitters 

The gain of the i.f. amplifier is 20 dB ±0.1 dB. The 
uncertainty is a long term known bias. The power 

splitters have a loss of 7.035 dB ±0.42 dB , but these 

uncertainties were known before flight. 

• Scatterometer processor errors (dynamic range error) 

The scatterometer processor error is 2.52 percent and 
distributed as follows : 

Mixer variation 82.9 percent 

Gate isolation 15.7 percent 

50 MHz ampere variation 1 . 4 percent 

Total 100.0 percent 

Component variations after the three-range-gate 
amplifier section are reduced by the gain of the 
amplifier in use since the variation only affects 
the finite range of computed after the coarse 

selection o’f the range of . For example, the 
range of will be known to lie between +10 to 

-10 dB, -10 to -30 dB, -30 to -50 dB exactly 
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by selecting the proper range-gate amplifier. The 
variational error will apply only to the resolution 
of the 20 dB range chosen. 

• Scatterometer error 

The error variations discussed thus far have been- 
shown as unknown bias errors which, while constant 
for a measuring sequence, are distributed randomly 
from one sequence to another. The best estimate of 
scatterometer accuracy is a root-sum-squared estimate 
of these errors as follows: 

Receiver signal error ±0.055 dB 1.2 percent 

(circulator unit) 

Transmitted power error ±0.059 dB 1.3 percent 

Calibration path error ±0.028 dB 0.6 percent 

Calibration error, ±0.03 dB 0.7 percent 

(undesired signals) 

Signal path loss error ±0.058 dB 1.3 percent 

Processor errors ±0.11 dB 2.52 percent 

rss = 3.46 percent 

or ±3.46 percent of the calibration signal. Besides 
this unknown bias uncertainty, one should consider the 
signal power contributed by undesired signals in 
relation to the desired signal. 
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Undesired (dB) 


Desired 

CdB) 

Differ- 

ence 

(dB) 

Error 

percent 

Switch isolation - 120 

dB 

-106 

-14 

±3.98 

Rotary joint leakage - 

137 dB 

-106 

-31 

±0.07 

Coupling between channels - 

-75.87 

-16.13 

±2,36 


92 dB maximum > 

The antenna gain error (squared) is ±12.1 percent. 

Total systematic (unknown 

2 2 

bias) error (rss) = (3.46 percent + 3.98 percent 

■ ■ '7 2 

+0.07 percent + 2.36 percent 

o 1/2 

+ 12.1 percent ) 

= ±13.4 percent 

The total unknown bias error estimate of ±13.4 per- 
cent would correspond to approximately ±0.55 dB. In 
addition to these systematic errors there are several 
sources of random errors. General Electric personnel 
(reference 48) have evaluated two error sources. 

These are: 

Statistical error due to noise ±12.2 percent 

Scatterometer processor ±2.3 percent 

(differential channel error) 

A significant error source is due to antenna pointing. 
Accurate value of the antenna pointing error is not available. 
However, ground tests conducted during Skylab missions indi- 
cate a pointing accuracy of 0.3° (reference 43). This angular 
error can cause the to be computed for incorrect angles. 

The error in computating due to Rg (the range) and 9 
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(the angle of incidence) are indeed small (less than 0.05 dB) 
for the highest 0) . But the error due to assigning an 
incorrect angle to a value could amount to a significant 

error depending on the ground scene and the angle of inci- 
dence. Maximum errors Will result for very flat surfaces. 

For actual rough ground surfaces this error is small (approx- 
imately 2 dB in the worst case) . 

To meet the constraints of the shroud envelope, the 
focal length to diameter ratio of the S-193 antenna was 
reduced. This factor, plus the limitations in the antenna 
feed and microwave switching network, resulted in low isola- 
tion between the vertical and horizontal antenna polarization 
ports. Based on the estimates of cross coupling from the 
University of Kansas, the ratio of power received in the 
desired polarization and the cross -polarization for the 
S-193 Radiometer was only approximately 10 to 13 dB . 

For the radiometer an attempt to make a first order 
correction for this mixing of energy from two polarizations 
was performed in production data processing. For the 
scatterometer the condition is more complex. Small errors 
are generated in the vertical- transmit , vertical-receive 
(VV) and horizontal- transmit , horizontal-receive (HH) modes. 
However, in the cross-polarized modes, vertical-transmit, 
horizontal-receive (VH) and horizontal- transmit, vertical- 
receive (HV) , the situation is extremely difficult. 

Cross-polarized return signals are generally 10 to 
15 dB below the level of like-polarized return signals 
(see data given in section 5). Since the antenna pro- 

vided only approximately 20 dB of isolation on transmission or 
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reception between polarizations, extraneous-like polarized 
return signals typically 5 to 10 dB below the power o£ the 
desired cross -polarized signals are introduced into the data. 

An undesired signal (5 dB lower in power) added to a desired 
signal of the same frequency can cause errors as large as 
+1 to -3.5 dB depending upon the phase relationships between 
the two. The stocl^astic nature of the return ' signals may 
contribute to increasing this error. Consequently, an 
uncertainty must be assigned to the cross -polarized scatterom- 
eter data. 

The preflight error estimates lend enough evidence as 
to the suitability of the results derived from the flight data. 
Considering the preflight performance, the precision better 
than 1.5 dB and accuracy better than 3 dB, as estimated in this 
report, are reasonable. It should be emphasized that for 
short segments of data (over a test site) precision as well 
as accuracy will approximate the preflight values given in 
appendix A. 
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7.0 CONCLUSIONS 

A precision estimate was generated by computing the 
standard deviation of the backscattering cross sections over 
selected uniform surfaces. A worst-case estimate of 1.5 dB 
was found. For the estimation of accuracy two methods were 
adopted. The Srl93 Scatterometer was compared to the 
existing microwave data at or near 13.9 GHz frequency. 
Theoretical values of were also computed using a com- 

posite scattering model and ground data. The comparisons 
showed the accuracy estimate to be better than 3 dB . It 
should be emphasized that these estimates apply to data 

for VV and HH polarization for the backscattering cross 
section range of 18 to -28 dB. For ' s greater than 18 dB, 
saturation of scatterometer further degrades the precision/ 
accuracy. For below -28 dB the analog- to-digital con- 

verter and the system noise combined degrade the accuracy 
and precision of the measurement. 

The dependence of backscattering cross sections on 
angle of incidence and surface roughness was found to be 
consistent with the theoretical and experimental results 
quoted in the literature. From the detailed analysis of 
S-193 Scatterometer data presented in section 5.0, it is 
possible to conclude that the quality of data is more 

than adequate for most applications for which the sensor 
was designed. In fact, significant applications have been 
reported in thu area of ocean and land remote sensing 
(reference 19} . For applications requiring an accuracy of 
more than 1 dB, the da .a must be evaluated carefully for 
specific data segments to ascertain the data accuracy. 
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One advantage of the S-193 Scatterometer precision/ 
accuracy analysis was in evaluating the algorithms and 
production data processing programs. ^As a result of data 
analysis which was undertaken for precision/accuracy 
determination, several anamolies were found. These 
anamolies were investigated and resulted in several modi- 
ficatio'ns to the NASA S-193 production processing programs 
The S-193 Scatterometer data was reprocessed for analysis 
by principal investigators. 

For the evaluation of the scatterometer algorithm, a 
method was developed to model the sensor on a time-sharing 
computer system. The effects of polarization isolation 
were studied using this model. With the existing polariza 
tion isolation of 20 dB , errors could result in the cross- 
polarized backscattering cross section. The magnitude of 
these errors depends on the polarization characteristics 
of the scattered energy from the target. 
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1 . 0 THE S-193 RADIOMETER/SCATTEPQMETER OPERATIONAL MODES 

The radiometer and scatterometer can operate in various 
scanning and polarization modes jointly and separately 
(reference 50) . A summary of these modes is given in 
table A- I and briefly explained in the following sections. 

1.1 INTRACK NONCONTIGUOUS (ITNC) MODE 


This mode is used for a joint radiometer and scatterom- 
eter operation. In this mode, only the pitch angle is varied. 
A resolution cell on the ground (figure A-1) is seen by the 
radiometer and scatterometer at approximately the following 
pitch angles: 0°, 15.6°, 29.4°, 40.1°, and 48°. The complete 

scan cycle time in this mode is 15.25 seconds. The roll angle 
is always zero. 



iMglW'e AT, Intrac'k noncontiguous flTNCj scan mode 
(n is the normal to the surface at point A^) . 
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TABLE A- I . 


NOMINAL S-193 RADIOMETER/SCATTEROMETER MODES 


Operation 

Scanning mode choice 

Polarization choice 

Pitch/roll angles 

Radiometer/scatterometer 

1. Intrack noncontiguous 
mode [ITNC} 

1. Scatterometer W, HH, VH, 
and HV 

Radiometer V and H 

2. One polarization combina- 
tion CW or HH or HV or 
VH) for scatterometer 
and V or H for radiometer 

0°, 15.6% 29.4% 
40.1% and 48“ 
(pitch) 

Radioraeter/scatterometer 

1. Crosstrack contiguous 
(CTNC) left/right 

2. CTNC, left 

3. CTNC, right 

Same as for ITNC 

Same as for ITNC 
(roll) 

Radiometer/scatterometer 

■ 

1. Intrack contiguous 
(ITC) mode 

1. One polarization combina- 
tion for scatterometer 
(W or HH or VH or HV) 
and V or H for radiometer 

Same as for ITNC 
(pitch) 

Radiometer/ scatterometer 

1. Crosstrack contiguous 
CCTC) left/right 

1. W or HH for scatterom- 
eter and V or H for 
radiometer 

+11“ to -11° 
(roll) 

Scatterometer only 

1 . CTC 

1. V and H radiometer data 

+11° to -11° 
(roll) 

Scatterometer only 

1. CTC 

1. Scatterometer data for 
W and HH 

+11° to -11° 
(roll) 




In reviewing the S-193 Radiometer/Scatterometer Skylab- 
acquired data, it was determined that some scan angle posi- 
tions in this mode were different from the nominal prelaunch 
values. The angles, 0®, 15.6°, and 29.4°, are not markedly 
different. However, the 40.1° and 48° angles show noteworthy 
change. In particular, the last angle remains, for most part, 
within 46 to 47°. The 40.1° angle is within 1.5° of the 
nominal value. 

On the Earth Resources Experimental Package (EREP) pass 
40 the antenna gimbals malfunctioned. The ITNC mode was not 
used subsequently. 

1.2 CROSSTRACK NONCONTIGUOUS (CTNC) MODE 

In this mode, the roll angle is varied identically to 
the intrack noncontiguous mode, and the pitch angle remains 
zero. The motion of the field- of- view (FOV) is shown in 
figure A- 2, where it can be seen that individual cells are 
viewed from only one antenna position. Because of the motion 
of the antenna in the pitch direction, the cells lie on a 
curved arc. There are three forms of this mode — left scan, 
right scan, and left/right scan as shown in the figure. The 
outermost cell is viewed at approximately 52° (corresponding 
to 48° gimbal angle) and the innermost cell at approximately 
0° at all times. The total scan time for a complete cycle 
is 15.25 seconds. The selection of polarizations is given in 
table A-T. 


The S-193 data for the CTNC mode shows that the antenna 
scan angles are approximately the same for 0°, 15.6°, and 
29.4° angles. The 40.1° pitch angle reaches only approximately 




Figure A- 2. — Crosstrack noncontiguous (CTNC) mode. 


5RljTPEODUGIBILITY OF Tii. 

PAXJE !S POOF 



37.5° for the Skylab -acquired data. The right scan extends 
up to approximately 43° instead of the nominal value of 48°, 
and the left scan extends up to approximately 46° instead of 
48°. Some oscillation in the antenna pitch angle is also 
noticeable at each dwell angle. However, actual antenna 
angles were recorded. Skylab-4 mission antenna scan motion 
was also variable. 

1,3 INTRACK CONTIGUOUS (ITC) MODE 

The pattern is similar to the intrack noncontiguous 
mode (figure A- 3), except that the antenna is scanned much 
faster and there is no dwell at any antenna pitch angle. 

The entire inflight path is eventually scanned at all inci- 
dence angles with this process. 



Figure A- 3. — Intrack contiguous (ITC) mode. 
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The scan cycle time is chosen so that ;>t the vehicle 
velocity the resolution cell at incidence angle 48° overlaps, 
the previous cell by approximate!) 25 percent, the 40.1° 
cell overlaps its predecessor by less than 20 percent, etc., 
down to the 0° incidence angle case where gapping rather than 
overlap occurs. The complete cycle for one scan takes 
approximately 4.0 seconds. 

f i 

As the vehicle progresses on successive scans, the entire 
path is viewed at- 48° and less, except for gapping at lowest 
angles. Table A- I gives the selection of modes and 
polarizations. 

In the ITC mode, the starting angle was about 43° [the 
nominal prelaunch value was 48°) during the Skylab-2 and 3 
missions. Since the Doppler filters are centered around 48°, 
the scatterometer data recorded for 43° is highly attenuated. 
Corrections for the Doppler filter attentuation have been 
implemented into the NASA/Data Systems and Analysis Division 
S-193 processing program. Other angles are also slightly off. 
The difference increases with increasing pitch angle. How- 
ever, no correction to the scatterometer data is needed at 
the angles other than the highest angle (approximately 43°) . 

During the Skylab-3 mission, a malfunction occurred in 
the antenna gimbals. The pitch gimbal was disabled as a fix. 
Consequently, no data was gathered in the ITC mode after the 
fix. 
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1.4 CROSSTRACK CONTIGUOUS (CTC) MODE 


This mode contains three submodes and further selection 
of polarizations (see table A-I). It provides a side-to-side 
linear scan covering '±11.375° and a turnaround to repeat. As 
can be seen in figure A- 4, this is a mapping mode. To compen- 
sate for the satellite forward velocity which could cause 
skewing of the pattern perpendicular to the flightpath, the 
pitch gimbal is scanned bacKwards slightly as the roll angle 
oscillates between its limits. Measurements are made for 
every 1.896° of beam center motion, ranging from -11.375° to 
+11.375° in roll. The total time of one cycle is 4.24 second. 
The pitch offset angles for this mode can be chosen as 0°, 
15.6°, 29.4°, or 40.1°. The roll offset angles can be chosen 
from 0°, +15°, -15°, -29.4°, and +29.4°. Either pitch or 
roll offset angle is selectable. 


FLIGHT DIRECTION 


148 km 
C^8 NMi) 

i 


t 


BEAM CENTERLINE 


4-4 - I - - 1 f i I M - - “4 - + 


-MEASUREMENT INTERVALS 
■185 km 'uioo NMi ^ 


Figure A-4. - Crosstrack contiguous (CTC) mode. 


A-7 



A study of Skylab-acquired data in the CTC mode revealed 
that the scan extends only up to a total of approximately 20 . 6 '^' 
instead of 22.75°. The repeatability of the timing sequence 
also differs from that indicated by figure A-4. 

Because of the antenna malfunction during the Skylab-3 
mission, the pitch gimbal was pinned at 0°. Consequently, 
the ground scans are not parallel for the Skylab-4 (SL-4) 
mission. The roll angles are also different in SL-4 data. 

t 

2.0 SUMMARY OF SCATTERQMETER SUBSYSTEM PARAMETERS 

The nominal scatterometer parameters, in addition to 

those described in 2.1, are as follows (reference 50): 

w 

Data rate (radiometer/ 10.66 kbps (effective) Biphase 

scatterometer) 


13.9 GHz 
TWT 

20 W (minimum) 
1.5 dB (maximum) 


Transmitter 

Frequency 

Output tube 

Peak power of tube 

Power losses to antenna 
input port 

Pulse width 
Pulse shape 

Pulse repetition frequency 


5.05 milliseconds for all scan 
angles in all modes 

100 ysec rise time, maximum 
100 ysec fall time maximum 

125 pps for all scan angles 
in all modes 
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Receiver 


Center frequency 
First i.f. 

Second i.f. 

System noise 

Second i.f. bandwidths 


Pitch angle (degree) 

0.0 

15.6 

29.4 

40.1 

48.0 

Number of i.f. filters per 
pitch angle 

Signal plus noise integra- 
tion times (milliseconds) 

Noise integration times 
(milliseconds ) 

Measurement precision 


Integration rate 


13.9 GHz 
500 MHz 
50 MHz 

1,200°K (maximum) 

Function -of pitch angle only 


Minimum i.f. bandwidth (kHz) 

68.4 
66:6 
61.0 
54.7 

47.5 

3 


41 (see table A-II) 


6 (see table A-II) 


Maximum standard deviation of 
0.0708 y, where y is mean of 
measurement, at 52° incidence 
angle, a = 30 dB . 

2.207 milliseconds wide, turned 
on 5.3195 milliseconds after 
start of transmit pulse for all 
scan angles except 48°, for 
which the width is 2.351 milli- 
seconds, turned on 5.61 milli- 
seconds after start of transmit 
pulse 


Detection 


Square law device 



TABLE A-II. - SCATTEROMETER INTEGRATION TIMES 


•1 = Maximum gain curve 
2 = Second highest gain curve 
3= Third highest gain curve 
4 = Lowest gain curve 


^Revised upon discussions with 
General Electric engineers 


Mode 

Integration 
time per 
pulse 

CTC scatterometex 
only 

2.187 ns 

ITC 48" 

1 

2.353 ms 

ITC less than 48* 

2.187 ns 

CTC radiometer/ 
scatterometer 

2.187 ms 


NC radiometer/ 
scatterometer or 
scatterometer only 
lo.o* 

2.187 ns 

IS. 6" 

1 

2.187 ms 

29.4° 

2.187 ms 

40.1° 

2.187 ms 

GO 

o 

0 

2.358 ms 


Nuaber of 
pulses 
integrated 



Total scatterom- 
eter (signal ♦ 
noise) integra- 
tion time (I.T's) 


17.496 

15.309 

13.122 

10.935 


21.177 

18.824 

16.471 

14.118 


19.683 

17.496 

15.303 

13.122 


30.618 

28.432 

26.244 

24.057 



1 ms for first pulse + 20.488 

4.872 ms for next four pulses 


50.301 

48.114- 

45.927 

43.740 


85.293 

83.106 

80.919 

78.732 


124.659 

122.472 

120.285+ 

118.098 


139. 968 
137.781 
135. 594 
133.407 


174.122 

171.769 

169.416 

167.063 



Scatterometer 
noise inte- 
gration time 
(I.T.)N 


6.813 ms 


N/A 


26.582 



10.22 


10.22 





















































Dynamic range 


-66.2 to -131.2 dBm, measured 
at antenna output terminals 
(65 dB range overall) 


3,0 THE S-195 SCATTEROMETER SYSTEM MATHEMATICAL MODEL 

3.1 GENERAL DESCRIPTION 
/ 

t 

The S-193 documents by General Electric Corporation 
describe the scatterometer system one part at a time. From 
these subsystem details a functional block diagram was pre- 
pared for the purpose of developing a sensor mathematical 
model. The sensor has been broken into a series of "elements" 
which represent accurately a particular component character- 
istic. By combining one or more of these elements, a component 
is simulated. The simulated components are combined into a 
subsystem. Together, the subsystems represent the system. 

The functional block diagram of the S-193 Scatterometer 
system is given in figure A- 5 . 

The input to the antenna is the vertically and horizon- 
tally polarized scattered power (Py, P^^ , respectively) . 

Passive microwave vertically and horizontally polarized powers 
emitted by the surrounding scene are denoted by Pj^y and 
, respectively. These are added to the radar return 
power. The cross -polar ization effects caused by the antenna 
and orthomode transducer are accounted for by introducing 
leakage factors and (figure A-6). 

per unit cros s -polar ized power (i.e., horizontal) present when 
the antenna is switched to receive vertically polarized 
signals. X^^^ can be defined similarly. 
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Figure A-5. — S-193 Scatterometer functional block diagram (concluded). 


































Figure A-6, ~ Cross -polari zation model. 


The following equations were used to take the cross- 
polarization into account: 




(A-l) 

*-*’h * *’nh^ 


(A-2) 


In arriving at equations (A-1) and (A- 2), the following 
assumptions were made: 

# The antenna is a linear element. 

• The ratios of the mixing of the polarization of the 
antenna remain constant and have been correctly 
measured. 
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• There is zero correlation between the vertical input 

t » 

signal, the horizontal input signal, the vertical 
noise received, and the horizontal noise received. 

The losses in the circuits (figure A- 7) were modeled 
using the following equation: 


P. P (L* - 1) 
1 , n 


(A- 3) 


P^ » Input power 



P. 

X 


Output power 




L* 

= Element 
Loss 



) _ Passive element noise 

.n power at temperature T^^ 


Figure A-7. — Loss model. 


In equation (A 3) P^ is the passive noise power due to 
the eleiiieit (loss L*) being at temperature Tj^ . The noise 
powers were computed using the Boltzmann's equation 


P 


N 


ii5 


n 


(A- 4) 
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watts 


where k is the Boltzmann constant (1.38 x 10 
sec/°K) , T the element temperature °K, and the noise 

bandwidth . 

The leakage factor has been modeled (figure A-8) using 
the following equation: 

Pq = P^(1 - X) + XPy^ (A- 5) 

where P^ is the leakage power into the element. 



Leakage 

factor 

X 

f 

P^ = Power in 

P = Power out 
0 


i 

II 

X 

IX 

Leakage power into device 


Figure A-8. — Leakage model. 


The amplifier gains have been treated by computing out- 
puts as the product of input and the gain. 

The equation used for the antenna voltage standing wave 
ration (VSWR) is as follows (see figure A-9) 
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(A- 6) 


P 


B 


4V 


V 


(1 + V ) 


2 


tVy - 1)^ 

CV^ + 1)^ 


where is the power in the cross -polarized channel leaked 

into the channel under consideration. 



ANTENNA 

VQWD 


P^ = Input power 

V C) n IV 

Pb 


' ' ■» 

Output power 


I P^ = Cross-polarized power 
Figure A-9. — VSWR model. 

The basic elements o£ the functional block diagram 
include the following types: 

• Cross polarization 

• Dissipative loss 

• Discontinuities (reflections) 

« Imperfect isolation (leakage) 

• '■"ystem noise 

• Gains (constant and temperature dependent) 

• DC offset 

• Integrator drift 
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• Square law detection 

• Analog-to-digital conversion 

In the scatterometer model the radiometric nois^ tempera- 

t . . ' 

tures have been converted to the noise powers. The advantage 
of treating the scatterometer in this manner is that the noise 
power for the system is adequately simulated. The output 
noise power from the scatterometer system depends not only on 
the system thermal state but also on the radiometric tempera- 
ture of the surrounding scene. 

The DC offset and the integrator drift were taken into 
account by adding voltages calculated from the following 
expression : 

TT 

(IT) ‘drift + (DC offset) 

where IT is the integration time and TC is the time constant. 
Proper integration times and time constants have to be taken 
into account for noise plus signal, noise, and calibration 
data. 


The computer model automatically selects paths and 
parameters to compute outputs corresponding to each S-193 
mode and submode. Amplifier gains, filters, integrators, 
command angles, scan modes, etc., are selected, depending on 
the desired output. The calibration data used in the mathe- 
matical model was taken from references 50 and 51. Tempera- 
ture dependence of the filter gains has been programmed. 
Since the DC offset term is only slightly dependent on the 
temperature, a constant temperature of 77°F was assumed. 
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This can cause a maximum error of 0.9 millivolts in the 
output based on actual results obtained from the computer 
model. The list of parameters used in the mathematical 
model is given in table A-III. This table also relates the 

t. 

symbols in the program to those in the scatterometer func- 
tional block diagram. The element losses, leakages, and 
gains have been assumed constant, as reported in reference 51 

The sensor simulation also included the analog-to- 
digital converter. Complete calibration data was not avail- 
able for the analog- to-digital converter. There are only 
20 calibration points instead of 1,024. In the analog-to- 
digital converter input, the voltage values are taken as 
each of the 10 bits is switched individually from 0 to 1. 
After the most significant bit reaches 1, each bit in turn 
is switched from 0 to 1 until a full count of 1,023 is 
obtained. In the simulation program since only 20 points 
were available, an interpolation technique, developed by 
Akima (reference 52) , was used to convert each input voltage 
to counts. Three A/D curves (reference 51} were used in the 
scatterometer simulation model. 

3.2 SCATTEROMETER COMPUTER MODEL AND RESULTS 

The analytical models presented in section 3.1 of this 
appendix were coded into FORTRAN (or XTRAN) statements for 
each functional block of the scatterometer model. These 
statemei r , subroucines, subroutine calls, selection logic, 
and other statements were combined into a computer program. 

A list of t?ie scatterometer model subroutines and a brief 
description of each are given in table A-IV (reference 53). 
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TABLE A-III. - LIST OF SCATTEROMETER PARAMETERS 


Computer 

symbol 

Mathematical 

symbol 

Value 

Comments 

PV 

Py 

1.0 X 10"^° 

1 . 

PH 

Ph 

1.0 X .10"^^ 

1 . 

TNV 

Pnv 

1.0 X 10^ 

1. 

TNH 


1.0 X 10^ 

1. 

PS CAT 


2.0 X 10^ 


IFIL 


1 

2 . 

IGAIN 


1 

3. 

I CMC A 


1 

4. 

IMODE 


1 

5. 

LAH 

^AH 

7.8412 X 10’^ 


LAV 

^AV 

6.0145 X 10'^ 


TNANT 

^NANT 

2.3665 X 10^ 

7 . 

LSAV 


1,055626 


LSAH 

T ^ 

^ AH 

1.072161 


TND 


2.9995 X 10^ 

6. 

LSD 


1.048409 


TNRFO 

^NRFO 

2.9775 X 10^ 

6. 

LSC 


1.027845 


LC 

^C 

1.23288 X 10"^ 


TNSCATX 

"^NSCATX 

2.8015 X 10^ 

6 . 

LSGV 

1 . . 

^*GV 

1.0218550 
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TABLE A- III. ~ LIST OF SCATTEROMETER 
PARAMETERS (Continued) 


Computer 

symbol 

Mathematical 

symbol 

Value 

Comments 

TNG 

T 

^NG 

2.991 X 10^ 


VH 

VSWR^ 

1.263 


TNLWG 

”^NLWG 

2.3665 X 10^ 


LSWGV 

T * 

^ WGV 

1.0 


LSWGH 

T * 

^ WGH 

1.0 


VV 

VSWRy 

1.262 


LGV 


4.39 X 10'^ 



gv 



TNE 

T 

NE 

2.9865 X 10^ 

6. 

LSGH 

T * 

gh 

1.01935 


LGH 


1.72 X lO'"^ 


TNALTX 

"^NALTX 

o 

X 

)-> 

o 

6. 

L12C 

^12c 

6.3895 X 10"^ 


TNL12C 

'^NL12C 

2.764 X 10^ 

6. 

TNLIWG 

"^NLIWG ■ 

2. 764 X 10^ 

6. 

LSIWG 

L* 

^ IWG 

1.051404908 


LSRJT 

T 

^ RJT 

1.051961874 


TNLRJ 

T 

NLRJ 

2.70 X 10^ 


LRJ 

'^RJ 

1.874 X 10‘^° 


TNL2WG 

'^NL2WG 

2.9775 X 10^ 

6. 

LS2WG 

^"*2WG 

1.037265659 
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TABLE A-III. - LIST OF SCATTEROMETER 
PARAMETERS (Continued) 


Computer 

symbol 

Mathematical 

symbol 

Value 

Comments 

TNLWS 

"^NLWS 

2.9775 X 10^ 

6. 

LSWS 

T * 

WS 

1.009252886 


LCC 

^CC 

1.90546 X 10"^ 


TN6 5 

T 

^N65 

2.9775 X 10^ 

6. 

TNWGS 

^NWGS 

2.9775 X 10^ 

6 . 

LIWGS 

^IWGS 

6.3095734 x 10‘^ 


L2WGS 

^2WGS 

5.6234 X lO"”^ 


Ld 

^d 

9.1877 X 10'^ 


TNF 

^NF 

2.992 X 10^ 

6. 

LSEANT 

L* 

EANT 

1.051869 


TNHOT 

”^NHOT 

3.922 X 10^ 

6. 

TNCOL 

'^NCOL 

3.183 X 10^ 

6. 

TNLSF 

"^NLSF 

2.992 X 10^ 

6. 

LSF 

L‘f 

1.028532 


LEA 

^EANT 

3.5495 X 10"^ 


VLPF 

VSWRj^pp 

1.1451 


LSLPF 

I 

LPF 

1.059863656 


TNSYS 

"^NSYS 

1.250 X 10^ 

8 . 

GTDA 

^TDA 

7.70312 X 10^ 

9 . 

TNLGPF 

"^NL6PF 

2.9775 X 10^ 

6. 
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TABLE A-III. - LIST OF SCATTEROMETER 
PARAMETERS (Continued) 


Computer 

symbol 

Mathematical 

symbol 

= Value 

Comments 

LS6PF 

^ 6PF 

1.462117 


TNLBPF 

"^NLBPF 

2.893 X 10^ 

6 . 

LSBPF 

T * 

^ BPF 

1.0 


GMA 

^MA 

4.7315 


TNLPD 

’^NLPD 

2.893 X 10^ 


LSPD 

^*PD 

1.0 


GIFPOl 

^IFPOl 

1.045 X 10® 

10. 

GIFP02 

^IFP02 

1. 

10. 

GIFPR 

^IFPR 

1. 

10. 

TNFILT 

"^NFILT 

2.7685 X 10^ 

6. 

LSPDFIL 

T * 

^ PDFIL 

1. 

11. 

TNLPCM 

'^NLPCM 

2.7685 X 10^ 

6. 

LSPCOM 

T * 

" PCOM 

1 . 

12 . 

GIF 

^IF 

1. 

10. 

GPLDET 

^PLDET 

1. 

13. 

IT 

^T 

1.74122 

14. 

TNWS 

'^NWS 

2.9775 X 10^ 

6. 

TN40 

^N40 

2.9775 X 10^ 

6. 

TC 

TC 

3.30 X 10^ 

15 . 

DRIFT 

DRIFT 

27 X 10"^ 

16. 

DC 

DC 

1 X 10’^ 

17 . 
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TABLE A-III. - LIST OF SCATTEROMETER 
PARAMETERS (Concluded) 


COMMENTS 

1. 

Inputs to simulated sensor 


2. 

LCF = 1 , MCF = 2 , HCF = 3 


3. 

1-+-Maximum gain, 4-^Minimum gain 


4. 

l-»-48% 2->40°, 3->29°, 4-^15°, 5-»-0° 


5. 

l-»-ITC; 2^ITNC; 3 , 4 , 5->-CTNC ; 6-»-R/S 

CTC; 7^S CTC 

6. 

Housekeeping data 


7. 

Computed from housekeeping data 


8. 

System noise temperature 


9. 

Gain of TDA 


10. 

Gain at i.f. 


11. 

Loss in power divider and filter 

(dB) 

12. 

Loss in power combiner 


13. 

Gain of square law detector 


14. 

Integration time 


15. 

Time constant 


16. 

Drift rate 


17. 

DC offset of integrator 
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TABLE A-IV. - SCATTEROMETER MATHEMATICAL MODEL SUBROUTINES 


Subroutine 

Purpose 

ACPT 

Acts as a driver routine to accept changes of 
parameter values (if any) , typed in from the 
demand terminal. i 


Input : 

Changes of program parameters. 


Output : 

Terminal display/printout of current 
value . 

ADCON 

Provides calibration data to convert an analog 
voltage to number of counts. 


Input : 

Analog voltage and filter number. 


Output ; 

Number of counts. 

GAINT 

Interpolates the filter gain ratio as a func- 
tion of temperature. 


Input : 

Command angle, temperature, filter, 
and mode number. 


Output : 

Gain ratio. 

GAINS 

Automatically selects gain value to be used 
as a function of polarization, mode, and input ' 
pcv:er P^ or P^^ . 


Input : 

Signal power, polarization mode, and 
current gain. 


Output : 

Gain factor. 
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TABLK A-IV. - SCATTliROMIiTHR MATHEMATICAL MODEL 
SUBROUTINES (Continued) 


Subroutine 

i 

Purpose 

AKITRP 

Uses the 20 input calibration points from A/D 
converter as a basis to interpolate the com- 
plete output calibration of A/D converter. 


Input : 

Calibration points and counts, 
analog input voltage, number of 
calibration points. 


Output : 

Number of counts. 

CX 

Computes 
element . 

the power output for a leakage 



PO = PI + (PX - PI)*AL 


Input : 

Input power (PI) , power through 
leakage "element" (PX) , the leakage 
(AL). 


Output : 

Power output from leakage element. 

CL 

Computes 
element . 

power at the output of a lossy 


PO = [PI + PN*(AL - 1)]/AL 


Input : 

Input power (PI) , thermal noise 
power of the "element" (PN) , and 
element loss (AL) . 


Output : 

Power output from lossy element. 
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TABLE A- IV, - SCATTEROMETER MATHEMATICAL MODEL 
SUBROUTINES (Continued) 


Subroutine 

Purpose 

VSWR 

Computes 

boundary. 

the output power at a reflection 


PO 

= PI*4*V/(V + 1)**2 


•• 

+ PR*[(V - 1)/(V + L)]**2 


Input : 

Input power (PI) , VSWR (V) , power 
presented to the function by the - 
thermal radiation of the components 
following the function (PR) . 


Output : 

Power at the "output" side of the 
reflection boundary. 

CPOL 

Accounts 
antenna . 

for the cross-polarization of the 


PI = (PV + PNV)*(1 - LAV) + LAV* (PH + PNH) 


P2 = (Pj^ + PNH)*(1 - LAH) + LAH*(PV + PNV) 


Input : 

Signal and noise powers, cross- 
polarization leakage. 


Output : 

Power at the output of the antenna 
for each polarization channel. 

PNOISE 

Converts the element temperature and effective 
noise bandwidth to noise power 



PN = (1.38*10’^^)*TN*BW 


Input : 

Element temperature and bandwidth. 


Output : 

Noise power. 
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TABLE A-IV. - SCATTEROMETER MATHEMATICAL MODEL 
SUBROUTINES (Continued) 


Subroutine 

Purpose , 

FILTER 

Table lookup to filter bandwidth. 

Input: 'Command angle, scan mode, and filter 

number. 

Output; Filter bandwidth. 

FILG 

Performs the normalization of the Doppler 
filter gains. 

Input: Signal power, command angle, scan 

mode, and filter number. 

Output: Signal power relative to the middle 

center frequency filter of the 0° 
Doppler filter set. 

DETTC 

Performs logical selection of time constant. 

Input; Scan mode , polarization, command 

angle . 


Output: Time constant. 

DETIT 

Performs logical selection of the integration 
t ime , 

Input: Scan mode, polarization, command 

angle, and signal gain. 

Output: Integration time. 
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TABLE A-IV. - SCATTEROMETER MATHEMATICAL MODEL 
SUBROUTINES (Concluded) 


Subroutine 

Purpose 

DDRIFT 

Determine voltage drift rate via a logical 
selection and table lookup. 

Input: Filter number and time constant. 

Output: Voltage drift rate. 

DDC 

Simple logical lookup for DC offset value. 


Input: Signal gain, filter number, and 


time constant . 


Output: DC offset value. 
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A complete listing of the S-193 Scatterometer computer 
model is given in table A-V. Several funs were made to 
illustrate the usefulness of the computer model. The results 
of the data runs are given in figures A- 10 through A- 13. 

The scatterometer output, calculated using the sensor 
mathematical model, for CTC radiometer/scatterometer mode is 
shown in figure A- 10. The input power level has been varied 
from -135 to -70 dBm. This graph has been drawn for the 
vertical polarization case. The radiometric brightness 
temperature of the sensed scene was assumed to be 4®K and 
270°K. Figure A-10 illustrates the effect of radiometric 
brightness of the target on the output of the scatterometer 
system. As expected, the brightness temperature is of signif- 
icance only for the lowest signals (or the highest gain curve). 
The linear behavior of the curves at high power levels is 
expected, since the model does not simulate nonlinear behavior 
at the dynamic range extreme. Note, however, that the output 
limits at 5 volts when the A/D converter hits maximum counts. 
The dashed part of the curves in figure A-10 is the expected 
output- versus input-power relationship since a linear system 
is simulated. 

To illustrate the effect of integration time, time 
constant, and filter characteristics, two input versus output 
plots have been drawn in figure A-11. These plots correspond 
to CTC and ITNC modes. In these computer runs the radiometric 
temperature of the surrounding scene is assumed to be 

270°K. The input versus output curves for CTC and ITNC modes 
are not the same. In the S-193 Scatterometer data processing, 
the difference in the input versus output is properly accounted 
for . 
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TABLE A-V. - SCATTEROMETER MATHEMATICAL MODEL 
COMPUTER PROGRAM 


'LIST 

1 

SCfiT SMD C8LLED SMfilH MFIY 7 1975 

L. 

IMPLICIT REfiL 

3 

COMMON /RNNOUT/ PINTCS) 

4 

COMMON /BIGOUT DIGD<6) 

5 

INTEGER IPOL 

i ■ 6 

COMMON /'KK/'PV.PH.TNV.TNL.PSCftT.IFlL. IGHIN. ICMBR. IM0DE.BW.>. 

-y 

LNH, LNV. TNNNT. LSfiV. LSRH. TND. LSD. TNRFD. LSC . LC> TNSCTX. LSb 


3 

TNG. VH> TNLW6. LSW6V. LSW6H. VV. LGV. TNE. LSGH. LGH. TNRLTX. L 1£ 

9 

TNL 1 MG » LS 1 MG . TNLR J . LSR JT . LSR JR . LR J . TNL2WG . LS2WG. TNLWS . L 

SMS. 

10 

LCC. TN65. TNWGS. LI MGS. L2WGS.LD. TNF. LSER. TNHOT. TNCDL. TNLS 

F . 

1 1 

LSF . LEH. VLPF . LSLPF . TNSYS . GTDR. TNL6PF . LS6PF . TNLBPF . LSEPF I 


1£ 

TNLPD.LSPD.GIFPOl.GIFPOe.GlFPR.TNFILT.LPDFIL.TNLPCM.LSP 

COM. 

13 

GIF. GPLDET .IT. TNWS . TN40. TC 

14 

COMMON^ KK 1 ^-DR I FT . DC 

15 

PV=1.E-10 1- 

16 

PH=l.E-10i 

17 

TNV*100. 

13 

TNH*100. 

19 

PSi;:RT=£0. 

£0 

IFIL=l 

3L 

IGRIN=1 

£2 

1CMDR=1 

£3 

IMDDE=i 

£4 

X : BWIN SET TO 0 

£5 

LRH=. 078412 

£6 

LRV=. 060105 

£7 

TNRNT=236.65 

£3 

LSRV=1. 055626 

£9 

L;SRH= 1.072161 

30 

TND=299.95 

31 

LSIi= 1.048409 

3£ 

rNRFD-297.75 

33 

LSC= 1 . 027845 

34 

LC=. 000123288 

35 

TNSCTX=<£80. 15 

36 

LSGV= 1 . 0£ 1 855 0 

37 

TNG=299. 1 

38 

VH= 1 . 263 

39 

TNLMG=236.65 

4 0 

LSMGV=l. 

41 

L$MGH=1. ’ 

4£ 

VV= 1 . £62 

43 

LGV=4.39 E-4 

44 

TNE=£98.65 

45 

LSGH=1 . 01935 

46 

Li5H=1.7£ E-4 

47 

TNRLTX==300. 

48 

L1£C=. 063895 

49 

TNL1£C=£76.4 

50 

TNL1MG=£76.4 

51 

LS1MG=1. 051405- 

5£ 

LSR JT*1. 051962 

53 

LSr JP»1.34?132 

54 

TNLRJ*£70. 

55 

LRJ= 1.874 E-10 
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TABLE A-V. - SCATTEROMETER MATHEMATICAL MODEL 
COMPUTER PROGRAM - Continued \ 


56 

TISL2W6-297.75 


57 

LS2WG»1. 037265 


58 

TNLWS*297.?5 


59 

LSU)S*1. 009253 


60 

LCC=1. 90546 E-07 


61 

TN65-297.75 


62 

TNWGS=297.75 


63 

117 LI WGS*6 . 3 095734 E-07 


64 

118 L2WGS=5.6234 E-07 


65 

LD*9.1877 E-4 


66 

TNF=299.2 


67 

LSEB=1. 051869 


68 

TNHDT»392.2 


69 

TMCDL*318.3 


70 

TMLSF=299.2 


71 

LSF=1. 028532 


72 

LEF(=3.5495 E-3 ' 


73 

VLPF*1. 1451 


74 

LSLPF»1. 059864 


75 

TNSYS=1250. 


76 

GTDfl«77 0.312 


77 

TNL6PF*297.75 


78 

LS6PF* 1.462 177 

• 

79 

TNLBPF=289.3 


80 

109 LSEPF-1. 


31 

GMR-4.7315 


• 32 

TNLPD=289.3 


83 

no LSPii-i. 


84 

Ill GIFP01»1. 064E10 


85 

112 GIFP02*!. 


■ 86 

113 6IFPR-1. 


37 

TNFILTa276.85 


88 

LPDFIL*!. 0 


89 

TNLPCM*276.85 


90 

114 LSPCOM*!. 


91 

115 G1F=1. 


92 

116 GPLLET=1. 


93 

IT=174. 122 


94 

TNIjJS=297.75 

. 

95 

TH4 0=297. 75 


96 

TC=33. 0 


97 

DRIFT=27. E-06 


93 

DC=4.9 E-03 


99 

% : MAIN PROG ENTERS HERE. 


too 

7. ; BCPT CONTROLS RECEIPT OF 

TERMINAL INPUT CHANGES TO Pfl 

1 RflMETERS. 

101 

V.l CBLL ftCPT 


102 

7. WRITE <6»^> 


103 

7: WRITE <6. ♦> '• INITIBL POWERS 

— PV*PH. . . 

104 

105 

106 

107 

WRITE <>6f ion PV»PH 


108 

101 FDRMRT<1H PV« '’»E15r7*' 

PH» dE15.7> 

1 09 

DO lOD lP0L«lf6 


110 

IF aP0L.GE.5> lCMDfl=5 


111 

IF dPDL.EQ.S) IFIL-2 


112 

IF <IPDL.EQ.6> IFIL-l 
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TABLE A-V. - SCATTEROMETER MATHEMATICAL MODEL 
COMPUTER PROGRAM - Continued 


U3 

114 

115 

116 
117 

■ 118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 
ISO 
181 
182 

183 

184 

185 


CRLL FILTER<IFIL.I MODE . I CtlDfl f By> 

CALL PHDISEaMU)S.BU.PNWS> 

CfiLL PMDISEah40.BW»PH40> 

C8LL PNDISE<TNV»BUpPNV> 

CRLL PNDISEaHH.BU.PMH> 

CfiLL PNOISEariflNTpBU.PHflMT> 

CfiLL PNOISE <TMDfBW*PMD> 

CfiLL PNOISE (TNRFOtBU.PNRFD) 

CfiLL PNOISE <TH8 fBU*PNG) 

CfiLL PNOISE <TNLU6.BW.PNLUG) 

CRLL PNOISE <TNE.BU.PNE) 

CfiLL PNOISE <TNfiLTX.BUI«PNRLTX> 

CfiLL PNOISE <TNL12 C»BWpPNL12C> 

CfiLL PNOISE <TNL1UIG.BW»PNL11JG) 

CfiLL PNOISE <TNLRJ.BW.PNLRJ> 

CfiLL PNOISE <TNL2UG*BU*PNL2WG) 

CfiLL PNOISE <TN65.BW.PN65> 

CfiLL PNOISE ONUIGSfBUiPNyGS) 

CfiLL PNOISE <TNFiBU.PNF> 

CfiLL PNOISE (TNHOTpBW.PNHOT) ' 

CfiLL PNOISE aNCOL»BU.PNCDL> 

CALL PNOISE <TNLSF.BW»PHLSF> 

CfiLL PNOISE aNSYSiBUfPNSYS> 

CfiLL PNOISE <INLBPFpBUpPNLBPF) 

CfiLL PNOISE «,TNLPD*BW.PNLPI>> 

CRLL PNOISE aNFILTpBW.PNFlLT) 

CfiLL PNOISEcTNSCTXpBWpPNSCTX) 

CfiLL PNOISE OHLPCtl.BU.PNLPCri) 

IF aP0L.GE,3> PV-0. 

IF <.IpaL.GE.3)PH»0. 

CALL CPUL <PV . PNV. PH . PNH » LfiV . LfiH . P 1 1 P2> 

CALL CLtPlpPNANT.LSflV.PVfiF) 

CALL CL<P2.PNftNTiLSAHpPHflF) 

CfiLL CL<PNRFDiPND»LSD»PCD> 

CALL CX<PCDpPNRFO»LD*PCD> 

CfiLL CUPCD.PNRFDpLSC.PCC) 

CALL CX(PCC.PNPFDiLC»PCC> 

IFaPOL.ECi.EJGO TO 2 
IF<IP0L.E«.4> GO TO 2 

7, : IP0L»2 IMPLIES HORIZONTAL 

CfiLL CL <PCC * PNG . LSGV s PLV> 

CfiLL VSUR<PHfiF»PLV»VH»VSWRH> 

CALL CU:vSURHiPNLU)G.LSUGHtP.'> 

CfiLL CL(PpPN6.LSGV.PP> 

CfiLL VSUIR<PVfiF»PP.VV.VSIORV> 

CALL CL (.VSWPY , PNLWGp LSWGVf PUG'> 

CfiLL CL<-PW6»PNG.LSGV»PLG> 

CALL CX<PLGpP»LGV.PVH> 

60 TO 3 

2 CALL CUPCC.PNGpLSGHpPLH) 

CfiLL VSWR<PVfiFf PLH pVVpVSWRV) 

CfiLL CL^VSI»)RViPNLWG.LSW6V.P> 

CfiLL CL<PpPN6f LSGH*PP> 

CRLL VSIJR<PHAFpPP.VH.VSWRH) 

CfiLL CUVSWRHpPNLWGpLSWGH.PWG? 

CfiLL CLc;PWGiPNG.LSGHpPCG> 

CfiLL CX<PCGpPpLGHpPVH'> 

3 CONTINUE 

CfiLL C.X<PNflLTXpPNSCTXpLMGpPNL12C > 

CfiLL CX<PNL12 CpPNSCTXpL12CpP12DBC> 

I F < I POL . GE . 5) CfiLL CX <PNflLTX p PSCHT p LUG p PNL 1 2CJ 
IF<IPDL.GE.5> CfiLL CX <PNL12CpPSCflTp L12CpP12DBC> 
CALL CL (P laiiBCi PNLIUG. LSIUGp PL 1UG> 

CALL CHPLHJGpPHLRJpLSPJTpPRJ> 

CALL cmPRJpPNL2WGpLS2UGpPL2UG> 

CfiLL CL<PL2 WGpPNUISpLSWSpPWS> 

PGP-PUS 

IFlIP0L.E0,3> PCP-PNUS 
IF».IPDL.EGI.4') PCP-PNUIS 
CfiLL CXilPNeSpPCPpLCCpPGSO 
IF<IPDL.LE.4.) P65C - PN40 
CALL CL (PVH'PNRFDpLSCpP40DB) 

CALL CX(P40DBpP65CpLCpP40DB) 





TABLE A-V. - SCATTBROMETER MATHEMATICAL MODEL 
COMPUTER PROGRAM - Continued 


CALL CL<P40DB.PMD»LSD»PCD> 

CALL CX<PNU)GS.PL2UGfLlWGS»PUGS> 

CRLL CX(PHUGS>PUGS>L2U6S»PUGS> 
IF<IPDL.UE.4.>PWGS - PNRFD 
CRLL CX<PCD.PWGS.LD»PCD> 

CRLL CL<PCD.PNE.LSER»PCE> 

PMHC- <PNHD7+PNCDL) /'2. 

CRLL CL<PNHC.PNLSFfLSF»PLF> 

CRLL CX<PCE»PLF»LEfi»PCE) 

CRLL VSWR<PCEf PNRFD »VLPF»PLPF> 

CRLL CLCPLPF.PNRFOiLSLPF.PLPF> 

PNOIS-PLPF+PNSYS 

PGTDR = GTDRePNDIS 

CRLL CL<PGTDR.PNL6PF.LS6PF.P6PF> 

CRLL CLCPGPF.PNLRJ.LSRJR.P2RJ> 

CRLL lX<P2RJ.PNL1U6.LRJ.P2RJ> 

CRLL CLCP2RJ.PNLBPF.LSBPF»PBPF> 
GPBPF»GMR»PBPF 

CRLL CL05PBPF.PNLPD.LSPD.PPD> 

PXFP01»6IFP01*PPti 

PIFPR-GIFPR*PIFP01 

P I FP 02*G 1 FP 02*P I FPR 

CRLL CLCP1FP02.PNFILT.LPIIFIL.PPDL> 

CRLL CL CPPDL.PMLPCM.LSPCDM.PPCL> 

CRLL GRINS CIPDL. IGRIN.GIF. PPCL> 
PIFG-PPCL*GIF 

CALL FlLGClCMIiR. IMDDE. IFIL.PIFG.PIFG> 

CRLL GRIHTCIFIL. IMDDE. ICMDR. THFILT. RRTID> 
PIF6 - PIFG*RttTID 

S : WRITE C6. '»> ' RRTID « '. RRTID 

CRLL DETIT CIMDDE. IPOL. ICMDR. I6RIN. IT> 

CRLL DETTC CIMDDE. IPDL. ICMDR. TO 
CRLL DDRIFTCIFIL.TC.DR1FT> 

CRLL DDCCIFIL.IGRIN.TC.DO 

IFCIP0L.LE.2> UIR I TEC6.*>' ICMDR- '.ICMDR.' < 
►IFIL-'. IFIL. ' ♦♦♦IGRIN-'. IGflIN 

: s WRITE<6.*> 'DC -'.DC. ' ♦♦♦ DRIFT-'. 

PPLD-P I FG* I T/TC4GPLDET 
PDfiD-PPLD+ClT*DRIFT CIT/TO *DC> 

PINTCIPDU - PDRD 

; I WRITEC6'»*> 'TC -'.TC.' ♦♦♦ IT -'.IT. 


« WRITE 
CONTINUE 
WPITE‘6.*> 
WRITE (6. ♦> 
WRITE C6.*> 
WRITE C6.*> 
WRITE C6.*> 
WRITE <6. ♦> 
CRLL RDCON 
WRITE C6.*> 
WRITE C6.*> 
WRITE <6. ♦> 
WRITE C6.*> 
WRITE <G.*> 
WRITE <6.*> 


CG.*> 'PINT C'.«IPDL.'>-'. PINT CIPDO 


'WERT- ' 
'VHDR - ' 
'VNDISEV- 
'VNDISEH* 
'VCRLl - 
'VCRL2 - 
CIFIL> 
'DWERT • 
'DVHDR - 
'DNVERT • 
'DNHDR - 
'DVCHLl ■ 
'DVCRL2 ■ 


.PINTC1> 
.PINTC2> 
'.PINT<3> 
'.PINT <4> 
'.PINT <5> 
'.PINTC6> 

' .DIG0<1> 
'.DIGD<2> 
'.DIED <3> 
'.01 GO <4> 
'.DIED <5> 
'.DIGDC6> 


SRLl'Si2.3 
JEDIT PNDISE 
''LIST 


I MR IN PRDG EXITS HERE. 


I • PNDISE' — COMPUTES NOISE POWER. 
subroutine PNDISE CTN.BW.PN> 

PH- Cl. 38 E-23>*TN*BW 

RETURN 

END 


SRU'Si .3 




TABLE A-V. - SCATTEROMETER MATHEMATICAL MODEL 
COMPUTER PROGRAM - Continued 


4 ED IT CPQL 



•^LlST 



1 

\ 

1 ^CPQL'^ -- COMPUTES PDLftRIZfiTlOM PDUER. 

2 


SUBROUTINE CPOL <PV> PNV» PH. PNH. L«V. LRH. P 1 . P2> 

3 


REAL LfiV.LftH 

4 


PI »< PV->PNV)*<l.-LfiV>+LflV»<PH+PNH> 

5 


P2 »<PH+PNH>*tl.-LHH)+ LfiH»<PV+PNV> 

6 


RETURN 

7 


END 

-'Q 



SRU^'Si .4 



!EDIT CL 



'LIST 



1 

54 

: 'CL' — COMPUTES POWER LOSS. 

2 


SUBROUTINE CL<TI.TE.fiL.TO> 

3 


TO » ai + TE ♦ <RL - l.>> FlL 

4 

54 

WRITE <6. ♦> ' CL- •'.TO 

5 


RETURN 

6 


END 

•'Q 


1 

SRU"Si.4 



•EDIT CX 



•'LIST 



1 

7. 

: 'CX' -- COMPUTES CROSSTfiLK. 

2 


SUBROUTINE CX <TI . TX.flL.TO> 

3 


TO - TI * <TX - TI> ♦ flL 

4 

54 

WRITE <6. ♦> ' CX- '.TO 

5 


RETURN 

€ 


END 

'0 



SRLI-'Si .3 



lEDIT VSUIR 



'LIST 



1 


: 'VSWR' — COMPUTES POWER DISCONTINUITY WITHIN MODEL. 

2 


SUBROUTINE VSWR <TI . TR. V. TO> 

3 


TD-T I *4 . ♦V/ <V+ 1 . > ♦♦2'tTR* < <V- 1 . > ✓ < V+ 1 . ) > ♦♦2 

4 

54 

WRITE <6. ♦> ' VSWR- '.TO 

5 


RETURN 

6 


END 

SRU-'SI .4 



lEDIT FILTER 


1 

54 

: 'FILTER' — DETERMINES FILTER BfiNDWIDTH. 

2 


SUBROUTINE FILTER<fiFIL.-RMODE.(RCMDH.BW) 

3 


DIMENSION RBW<3»5) 

4 


DflTR RBW/34e90. .64445. .64930. .71270. >70870. .73775. .7950 

5. 



5 

79330 

. . 80465. . 84570. . 87405. . 84990. . 89425. . 89750. . 91995. / 

S 


IFlL-ftFIL 

■ 7 


ICMDR-FiCMDfi 

3 


IMDDE-fiMODE 

9 


SO TO <l-l.2.2.2.1»l> .IMODE 

10 

1 

BW-RBWCIFIL. iCMDfi) 

11 


RETURN . 

12 

ic 

BW-BBWaFIL»5> 

13 


RETURN 

14 


END 

"'Q 



SRU'’S» .4 






V 


TABLE A-V, - SqATTEROMETER MATHEMATICAL MODEL 
COMPUTER PROGRAM - Continued 


!EDIT 

GftINS 

i 




1 


8 '•GfilNS:'’ ~ LOGICfiLLY SELECTS flUTOMfiTIC GHIN F8CTDR. 

a 


SUBRDUTIME GHINSaPOU.flGRIN.GIF.PPCL) 

3 


DIMENSION 6flIN<.4> 

4 


DRTft GfilN^'l. .. 1007135f. 0082846. . 0008044/ 

5 


I GAIN *«GflIN 

6 


GO TOO. 1> 2.2.3. 3;>.IPQL 

7 

1 

IF«IG81N-10.:) lu.20.20 

9 

10 

IF<PPCL.LE. .94) IGflIN-1 

10 


IF<PPCL.GT. .94) IGflIN-2 - • 

11 


IFCPPi:L.GT.9.333) IG8IN-3 

1£ 


IF<PPCL.GT, 113.4607) 1G8IN-4 

13 


GO ro <=' 

14 

2 

IGRIN-l 

15 


GO TO 9 

16 

3 

IG8IN-4 

17 


60 TO 9 

18 

20 

IGRIN « I6RIN -10 

19 


GIF -GRIN* IGRIN) 

ao 


RGRIN * IGRIN , i 

ai 


RETURN 

22 


END 

SRU'Si 

.4 


SEDIT 

FILG 


■'L 

1 


8 'FILG- — EQUALIZES GRIN BETWEEN FILTERS. 

£ 


SUBROUT I NE F ILG '^RCMDR. RHODE . RF I L . P I FG . PDF6) 

3 


DIMENSION R6AIN<3.5) 

4 


DRTR R6RIN/1. 15. 1.15. 1.10. 1.16. 1.15. 1.15. 1.16. <5.1.20 

5 

1 

.21. 1. 17. 1.21. 1. 05. 1. 02. 1. 06/ 

6 


ICMDR-RCMDfi 

7 


IFIL- RFIL 

3 


IMDDE-RMDDE 

9 


GOTO a. 1.2. 2.2. 1. 1) .IMODE 

10 

1 

PDFG-P I F64RGR I N H F I L . I CMPR) 

11 


1 WRITE<6.*) 'R6RIN<IFIL»ICMDR> -''.RGRINOFIL. ICMDR) 

12 


RETURN 

13 

2 

PDFG-PIFG-RGRIN'IFIL.S) 

14 


1 WRITE(6.*) RGRIN<IFIL.5) ■«' . RGRIN < IF IL. 5) 

15 


RETURN 

16 


END 

'Cl 


, 

:ru'ss 

.4 


•EDIT 

DETTC 


1 


8 DETTC' — DETERMINES TIME CONSTANT. 

2 


SUBROUTINE DETTC <AMDDE. IPOL.RCMDR. TO 



ICMDR-RCMDH 

4 


IMDDE-RMODE 

5 


IF * IP0L.GE.5) GO TO 20 

6 

10 

GO TO a 1.12. 12. 12. 12. 16. 17) .IMODE 

7 

1 1 

TC-4.0 

8 


RETURN 

8 

12 

TC-33. 0 

lb 


IF aCMDfi.EQ.5) TC-10.22 

11 


RETURN 

12 

16 

TC-10.22 

13 


RETURN 

14 

15 

17 

TC-4. 0 

16 


RETURN 

17 

20 

TC-10.22 

18 


RETURN 

19 


END 

. 'Q 



3RU- SI 

.4 



REPRODUCIBILITY OF 
PAOF TF ’ 
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TABLE A-V. - SCATTEROMETER MATHEMATICAL MODEL 
COMPUTER PROGRAM - Continued 


JEDIT 

DETIT 





1 

V. 

: ''DETIT^ -- DETERMINES INTE6RRT1DN TIME. 

•*% 

k. 


SUBROUTINE DETIT <RMODE* IPOL.RCMDRf R6RIN»fiT> 

3 


IMODE*RMODE 

4 


ICMDR=RCMDfl 

5 


flITPP-£. 187 

6 


IF<ICMDR.EQ. 1> RITPP*2.353 

7 


IF <IP0L.6E.5) GO TO £0 

8 


GOTO (:il»l£»12.1£.l£»16»17> » I MODE 


11 

RT=R I TPP» CIO. -BGR I N> 

lb 


IF<IPDL.GE.3> RT*£7.063 

11 


RETURN 

12 

1£ 

GO TO ''!1£1 » l££f 1£3> 1£4» i£5> > ICMDR 

13 

1£1 

RT=R I TPP* <75 . --RGR I N i 

14 


IF <IP0L.GE.3) RT=1£5.53£ 

15 


RETURN 

16 

1££ 

RT»R 1 TPP* <65 . -RGR I 

17 


IF<IP0L.GE.3) RT»1£5.53£ 

18 


RETURN 

19 

123 

RT-RITPP* <58-R6RIN) 

£0 

* 

IF<IP0L.GE.3> RT-57.990 

£1 


RETURN 

££ 



£3 

1£4 

RT=RITPP*<40.-RGRIN) 

24 


IF<1P0L.GE.3> RT*57.990 

£5 


RETURN 

£6 

125 

RT=fiITPP* <24-RGRIN:> 

£7 


IF<IPDL.GE.3» RT*£4.094 



RETURN 

£9 

16 

RT=R 1 TPP* < 1 5-RGfl 1 N> 

30 


IF <IPDL.GE.3> RT=13.686 

31 


RETURN 

3£ 

17 

RT*RITPP*<9-RGRIN> 

33 


IF <IPDL.GE.3> RT=6.544 

34 


RETURN 

35 

£0 

RT* £0.488 

:36 


RETURN 

37 


END 

SRU'S: 

.6 


lEDIT 

DDRIFT 


1 

'4 

: 'DDRIFT'- ~ DETERMINES DRIFT RRTE. 

£ 


SUBROUTINE DDRIFT <RFIL»TCf DRIFT) 

3 


DIMENSION DRRTE<3»3) 

4 


IFIL»RFIL 

5 


DRTR DP.RTE'' . 163» ,£10» .£06» . 060f . 083» . 079> . 017» . 0'£7» . 031 

6 


IF<TC.EQ.4) ITC=r 

»' 


IF < Tr.EQ. lo.iea.) ITC*£ 

8 


IF <TC.EQ.33> ITC=3 

9 


DPIFT*DRRTE<IFILf ITO*. 001 

10 


RETURN 

1 1 


END 

-Q 



SRU-'S: 

.4 
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TABLE A-V. - SCATTEROMETER MATHEMATICAL MODEL 
COMPUTER PROGRAM - Continued 


1 {EDIT DDC 1 

1 


: 'DDC" — DETERMINES DC OFFSET VfiLUE. 

2 


SUBROUTINE DDC '^flFIL.FlGRIN* TC, DO 

3 


DIMENSION D0FFC3,4,3> 

4 


DRTfl DDFF/6.4,4.6,4.8,3.8,4. 0,4. 0*3. 5, 4. 0,3. 9, 5. 0*4. 0*3 

I 1 

5 

6 

.4,5. 0,4. 7, 3. 8*4. 4, 3. 9*3. 5*4. 4, 3. 8, 5. 0*4. 4, 3. 5* li 

6 

6 

. 1*4. 9*4. 4, 3. 8, 4.3, 3. 7, 3. 5, 4. 3, 3. 6, 5. 0i*4.3,3. 3/ 

7 


IFIL=RFIL 

8 


IGflIN>= RGRIN 

9 


IF<TC.EQ.4> ITC»1 

10 


IF<TC.EQ. 10.££) ITC=£ 

11 


IF <TC.EG!.33) ITC*3 

12 


DC*. OOl^DDFFCIFIL, IGfilN, ITO 

1-3 


RETURN 

14 


END 

1 ;PU'SI.4 1 

1 'EDIT 8DC0N 


1 


! 'flUCON' — RNRL06 VDLTR6E TO DIGITRL COUNT CONVERTER. 

p 


SUBROUTINE RDCON<RFIL> 

3 


COMMON /RNROUT /fiV V , RVH , RNV * RNH , RC 1 * RC£ 

4 


COMMON /D I GOU r /DVV , DVH* DNV , DNH , DC 1 * »C£ 

5 


RERL VOLTS <£0>* VOLTK£0> , VOLTE <£0) , VOLTS c;20:) 



RERL CNTS<£0> 

1 


DRTR VDLTl/-. 004, . 002* .007* . 0206* . 0373* . 0765, . 154* .310, 

1 . 8242 » 1.245 

*2.497*3.744* 4. 37, 4. 684, 4. 84, 4. 9£* 4. 958, 4. 984* 4. 992, 5 

1 . 003^- 

8 


VOLT£>' - .007»*001p. 0065 f • 0156? • 03‘;i? • 0737 ? • 15 1 ? . '309? I 

1 • iS^ Od 9 % 1 


1 . 

246 *2.493*3.746*4.37*4. 684 * 4.34,4.919*4. 959 * 4 . 975 * 4 . 995 * 

5.0 05/ 

10 


DRTR VDLT3/-. 003* . 003, . 0083, . 0209* . 0422* . 0797, . 1573* .31 

5 » . 623* 1.2 


44* 2. 50 1,3. 748* 4. 37* 4. 683. 4. 84 1*4. 923* 4. 956* 4. 977 *4. 98 

9,5.010/ 

11 


DRTR CNTS/0. . 1 . , 2. , 4. , 8. , 16. *32. *64. * 128. ,£56. ,512. *768 

. . 896. * 960. 


*992. , 1008. * 1016. , 1020. * 1022. * 1023./ 

12 


IFIL=RFIL 

13 


GO TO (1111*2222,3333), IFIL 

14 

nil 

DO 991 I =1,20 

15 

991 

VOLTSd;- = VDLTKl) 

16 


GO TO 999 

17 

££££ 

DO 992 I = 1,20 

13 

992 

VOLTSd) = V0LT2(I:) 

19 


GO TO 999 

20 

3333 

DO 993 I = 1*20 

21 

993 

VOLTS < I) = VDLTSCn 

22 

999 

CRLL flKITRP(.VDLTS.CNTS*RVV*DVV*20) 

23 


CRLL RK I TRP <VOLTS * CNTS* RVH * DVH * 2 0) 

£4 


CRLL RK I TRP < VOLTS , CNTS * RNV * DNV * 2 0) 

25 


CRLL RK I TRP (VOLTS* CNTS* RNH* DNH *20) 

£6 


CRLL RK I TRP (VOLTS * CNTS , RC 1 * DC 1 , 2 0) 

27 


CRLL RKITRP(VOLTS*CNTS*RC2*DC£*20) 

£8 


RETURN 

£9 


END 

I SRU'Ss.5 


} 
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TABLE A-V. - SCATTEROMBTER MATHEMATICAL MODEL 
rOMPUTBR PROGRAM - Continued 


I EDIT WKITRP 

■'U 


I '«KITRP' — INTERPOLftTIVE LINE SMOOTHING^ 

SUBROUTIME ftKITRP<XXtYY.XD#YDiH) 

REHL T <2> . M <5> p XX < 1 > . Y Y < I ) i XO . YO • X <6> I V <6> 

IF(N.6E.-3> GO TD 1 
WRITEClp 1000>N 

FORMAT OSH BKITRP - N MOST BE .GE. « N - pI5> 

STOP 

DO S I « SpN 

XXMXD - <XXvI-l>-XD>*<XX<I)-XD) 

IF<XXMXO.LE.O.> GO TO 33 
CONTINUE 

WRITE<fcH001)XX<l> »XX<N) tXD 

FORMATdH t' AKITRP XD MUST BE BETWEEN '»EI5,7t' AND 


'.E15.7, '' 

■ 

a 

X 

',E15.7> 

15 


URITE<6.4> 'PV - ',PV,' PH -',PH 

IG 


STOP 

17 

33 

CONTINUE 

18 

3 

K » 3 

19 

4 

X<K) - XX<I+K-41 

2D 


Y<K> ■ YY<I+t<-4> 

21 


K - K - 1 

as 


IF<K.LE.0> GO TD 5 

23 


IF<I+K-4.GE. 1) GO TD 4 

24 


IFCK.EQ.l) 60 TD 6 

25 


X (2) »XX < 1 > +XX <2> -XX <3> 

2G 


xci> - a.*xxa>-xx<3> 

27 


YCS-) ■YY<1)> <XC2>-XX<l>>*i:2.*<YY<l)-YY<2>>/<XX<t)-XX<2>> 

- <YY <2> X 



28 


-YY <3) > ^ <XX <a'> -XX <3> ) ) 

29 


Y(i>-Y<2:»4'<xa>-x<a).>4<<:Y <a>-YYa>>/<x<a>-xxa^> + <YY<i.> 

-YY <2> > ^ 

X 


30 


<XX C1>-XX<3>>> 

31 


GO TD 5 

32 

G 

X(l>-X <2>*X<3> -XX <3> 

3.3 


Y < 1 .) -Y C2> <X <3J -X <2> ) ♦ <2 . ♦ <Y <2> -Y <3> > / <X -X 0> > - <Y O) 

-YY <3> > /• 

H 


34 


<X<3>-XX<3>>> 

35 

5 

K » 4 

36 

7 

X<K> » XX<’I+K-4> 

37 


Y<K> » YVa+K-4) 

38 


K =* K + 1 K ' 

39 


IFCK.GT.ej GD TD 10 

40 


IFa+K-4.LE.N> GD TD 7 

41 


IFCK.EO.G) GD TD 9 

42 


X(5> - X <,4) ♦X 0> -X 0> 

43 


Y O) »Y (4) + (X O) -X <4) 1 ♦ <2. ♦ <Y (4 > -Y O) ) (X <4> -X O) > -7. 

44 


CY ';3> -Y 0> > 7 CX 0> -X <2) ) > 

45 


X<6)-t. >X<4')-X<2> 

46 


Y <6^ -Y <5i + <X <6> -X <5) > ♦ < <Y <5> -Y <4> > / <X <5> - X <4> > + tY ^4> -Y < 


<X < 4 > -X <3J > - <Y 0> -Y <2> > ^ <X <3> -X <S) > > 

GO TD 10 

X<G> -X <5> +X (Ai -X <3) 

Y (G-) -Y (S) + <X CG) -X <5) > • <2, ♦ <Y <5> -Y <4> > / (X <5> -X <A> > -X 
<Y «■> -V O) ) ^ CX C4> -X <3) ) > 

DO 11 K • 1,5 

M (K> - <,Y <K) -Y <K+ 1> ) / <X <K) -X 1 > ) 

DO 12 K » liS 

T - <ABS <M <F+3> -M <K*2> > *M 1? +ABS <M i K+ 1 > -M «.tO > <K+2) 

<ABS <M <K-t-3> -M CK+2) > +RBS CM <K+1> -M <K> » 

IF<C*'KKEQ.M<K+1>> .AND. (M<K'*-2> .EQ,n<K+3;>>> V. 

T <K> - 4-M CK+2> > /-a. 

CONTINUE 
PO » YO> 

PI « T<1> 

P2 - <3.*<Y<4>-Y<3>>/'<X<4>-XO'>>-2.*TCl-)-T<2>)^<X<4>-X<: 
P3 - a<l)+T<2>-2.*CY<4>-YO>>/'<X<4>-XO>>)/'<X<4>-XC3>') 


DX - XQ-X<.3> 

DXDX - DX*DX 

YO-P 0+P 1 ♦DX+P24DXDX+P3*nx*DXDX 

RETURN 

END 
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TABLE A-V. - SCATTEROMETER MATHEMATICAL MODEL 
COMPUTER PROGRAM - Concluded 


lEDIT 

1 

GRINT 

1 » 

1 


! 'GRINT" -- DETERMINE TEMPERRTURE BIRSED FILTER GRIN R 

FITID. 





SUBROUT I NE GR I NT <RF I L » RHODE » RCMDR f TEMP > RRT I D) 

3 


DIMENSION RRTID6<3.3.5> 

4 


DRTR RRTIDG/.882* 1. * .865* .932t 1. . .889it .865 j 1. > .839* 

5 


.977* 1. * .863* .933* 1. * .888* .837* 1. * .831 « .984* 1. * .862 

6 

1 

* . 944* 1 . > . 882* . 860* 1 . * . 826* 1 . 04* 1 . * . 851 * . 984* 1 . * . 875* \ 

7 


. 872* 1 ... 822* 1 . 12* 1 . * . 837* .971*1.*. 878* . 964* 1 ... 819/ 

8 


IFIL ■ RFIL 

8 


I MODE - RHODE 

10 


ICMDR - RCMDR 

It 


TEMPIN TEMP - 273.15 

1£ 


IF<TEMPIN-24.44> 10*20*30 

13 

10 

TMPINR =1.0 

14 


IFaEMPIN-*-10.) 12. £2. 32 

15 

12 

URITE<6.»> '■♦♦WRRNING^^SCRT TEMP LESS THRN -10DE6 CENT. 

1 TEMPIN *■'» y. 


16 


TEMPIN 

17 

32 

GO TO 8 

18 

22 

GOTO (4. 4. 5. 5*5*4. 4> * I MODE 

18 

4 

RRTID = RRTIDG<1* IFIL* ICMDR) 

£0 


RETURN 

£1 

5 

RRTIO * RRTIOGYl* IFIL.5) 

££ 


RETURN 

£3 

£0 

RRTIO =1. 

£4 


RETURN 

25 

30 

TMPINR =3. 



lF<TEMPIN-30.) 13.23*33 

27 

23 

GO TO (.6* 6. 7. 7* 7* 6* 6) * I MODE 

£3 

6 

RRTIO = RRTIOG <3* IFIL* ICMDR) 

£8 


RETURN 

30 

7 

RRTIO = RRTI06<3* IFIL.5) 

31 


RETURN 

32 

33 

WRITE<6.*v '■♦#(;IRRNING**Se:fiT TEMP EXCEEDS 30DEG CENT. TEM 

PIN 

y % 


33 


TEMPIN 

34 

13 

GO TO 9 

35 

8 

GO TO Cl. 1.2. 2.2* 1* 1) .IMODE 

36 

1 

RRTIO - RRTIDGCTMPINR* IFIL* ICMDR) 



GO TO 3 

38 

£ 

RRTIO = RRTIDGCTMPINR* IFIL.5) 

38 

3 

SLOPY » <1. -RRTIO) /34. 44 

4 0 


TEMPIN » 24.40 - TEMPIN 

41 


BYY * l.-CSLOPY»TEMPIN> , 

42 


RRTID ■ BYY 

43 


RETURN 

44 

9 

GO TO <.19.19.29.29*29*19*19) .IMODE 

45 

18 

RRTID = RRTIOGaMPINR* IFIL* ICMDR) 

46 


GO TO 39 

47 

28 

RRTID = RRTIOGaMPINR* IFIL.5) 

48 

39 

SLOPY = d. -RATIO) /5. 56 

48 


TEMPIN « TEMPIN -24.44 

50 


BYY » l.-(SLDPY4TEMPIN) 

51 


RRTID -BYY 

52 


RETURN 

53 


END 

SRU'Si 

.8 


!DFF 



IJSfiGE 

□N 05/07/75 RT 07850125 1 

SRU'Si 11.6 

ELAPSED TlMEi OOi 23108 

GOOD BYE. 
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Figure A-10, — S-193 Scatterometer volts output versus 
power input for two radiometric antenna temperatures. 
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Figure A-12. — S-193 Scatterometer volts output versus 
power input for a difference of 15 dB between 
vertical and horizontal input power. 
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The affect of the cross-polarization leakage on the 
S-193 Scatterometer output is illustrated in figure A-12. 

The input versus output has been drawn for two cases. In the 
first case the vertical incoming power (P^) is 15 dB higher 
than the horizontal incoming power (Pj^) . The output vertically 
polarized voltage is not significantly different from 
Py = case. The reason for this is that the leakage from 
horizontal is very low because of the lower magnitude of the 
horizontally polarized power. In the second case, the 
vertically polarized incoming power is 15 dB below the hori- 
zontally polarized power. Now the leakage, due to higher 
Pp^ , is significant. The vertically polarized output is the 
sum of two contributions — vertical incoming power and leakage 
from the horizontal port. This leads to an increase in the 
output power. 

There is another way of interpreting figure A-12. For 
the unlikely cases where the horizontally polarized power is 
15 dB above the vertically polarized power, the proper curve 
should be used to determine the incoming vertical power from 
the sensor output. If this is not done and the output is 
calculated using the characteristic curve where Py is large 
(15 dB or more) as compared with Pj^ , then an error will 
result. This resultant error will depend upon the magnitude 
of the input signal. The cross-polarization assumes signifi- 
cance since the difference in horizontally and vertically 
polarized incoming power is unknown for remotely sensed rough 
scenes. 

Serious errors can result from using one curve for the 
calculation of input power from the output measured voltage. 

It is interesting to note that the errors in estimating the 
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input power will decrease i£ the incoming vertically and 
horizontally polarized powers differ by a less amount which 
is the usual case. This is why only one set of curves has 
been used in the S-193 production data processing. This is 
illustrated in figure A-13, where Py and differ by 

5 dB for each plot. In figures A-11, A-12, and A-13, the 
radiometer temperature of the surrounding scene was taken to 
be 270°K. 

The scatterometer mathematical model presented in this 
report is not valid for the Skylab-4 (SL-4) mission. This 
is due to a drastic change in cross -polarization and antenna 
pattern. It is, however, possible to modify the present 
scatterometer computer model to reflect these changes. 
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DATA PROCESSING TECHNIQUES 



APPENDIX B 

DATA PROCESSING TECHNIQUES 
1.0 S-195 SCATTEROMETER DATA PROCESSING 

1.1 S-193 RAD IOMETER/ SCATTEROMETER PROCESSING FLOW 

The first step in the Skylab-acquired data processing is 
to duplicate the 2 8- track Earth Resources Experiment Package 
(EREP) tape. From this tape two 14- track tapes are developed 
(figure B-1) . These tapes contain Interrange Instrumentation 
Group Format A (IRIG A) time converted from the original time 
words on one of their tracks.. One of the 14-track tapes con- 
tains the S-192 sensor data and the other contains S-190A, 
S-191, S-193, and S-194 sensor data. From the latter 14-track 
tape (the S-193 data on this tape is on a frequency modulated 
subcarrier), a nine-track digital S-193 Radiometer/Scatterom- 
eter tape is prepared. 

This tape (also known as S061-4 product) is time edited 
and has data in raw pulse code modulated (PCM) counts. 

From S061-4 the raw processed tape (S061-2) containing 
data in PCM counts is generated directly (figure B-2). Cali- 
bration data and Skylab Best Estimate of Trajectory (SKYBET) 
Ephemeris Data tapes are used in data processing to generate 
computer-compatible tapes (CCT's, S061-1), tabulations 
(S062-A, B, C) , and plots (S063-D, 6). The output product 
data flow is shown in figure B-2. The detailed definition of 
each product is given in Earth Resources Production Processing 
Requirements for EREP electronic sensors document 
(reference 40). The production data processing program can 









Figure B-2. — S-193 Radiometer/Scatterometer output product data flow. 












also be used to generate stripcharts, oscillographs, and 
event records (S064-1, 2, 3). 

1.2 CALCULATION OF SCATTEROMETER BACKSCATTERING 

CROSS SECTION 

The geometry of a resolution cell for S-193 Scatterom- 
eter is shown in figure B-3. Point C represents the inter- 
section' of the antenna boresight line with the Earth's 
surface. The angle a between the boresight direction and 
the Z-local vertical (figure B-3) depends on the roll (4>^) 
and pitch (4)^) angles of the antenna. In terms of the pitch 
and roll angles, a is given exactly by 

2 2 2 

tan a = tan 4> tan <p (B-1) 

r p 

The range of values of and cj)^ extend up to a 

maximum of 48° (appendix A). The angle of incidence 6 is 
defined as the angle between the local normal to the Earth 
at the point C and the vector described by joining the SRylab 
S-193 antenna position S to C. 


The power per unit area at a range r from the trans- 
mitting antenna is given by 



PtG, 


47rr' 


(B-2) 


The power AP^ is incident on the surface of the Earth. 

P.J, is the power radiated by the transmitting antenna in watts. 
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Figure B-3. — S-193 Scatterometer resolution cell geometry. 


and G,p is the gain o£ the antenna in the direction of the 
resolution cell. This incident power is scattered in the 
space surrounding the resolution cell. The scattered power 
at the surface can be expressed as 




(B-3) 


where is the scattering cross section of the surface 

per unit area. AP^ travels back to the receiving antenna 
rthis can be the same as transmitting antenna but polariza- 
tion of reception may be different) . The power per unit 
area at the receiver is given by 


AP 


ar 


AP. 


4TTr ‘ 


(B-4) 


The receiving antenna has an effective receiving aper- 
ture or effective area . Thus, the received power from a 

unit surface area is 



AP. 


47Tr‘ 



(B-5) 


Combining equations (B-2) through (B-5), the power 
received as a result of scattering from unit surface is 
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The S-193 Scatterometer operates essentially in a con- 
tinuous wave (CW) mode. The carrier frequency of 13.9 GHz 
is pulse modulated (125 pps) with a 62.5 percent duty cycle. 
This assures small modulation sidebands and the transmitter 
wavelength (A.) can be assumed a constant. 


The receiving antenna aperture is given by 



(B-6) 


where is the gain of the receiving antenna. From equa- 

tions (B-5) and (B-6) 


AP = 
r 


!t^o35 

(4tt)^ r"^ 


(B-7) 


The total power received (P^^) from an illuminated area can be 
calculated from equation (B-7) by integrating over the 
reflecting area: 






/ 


G . G o 
t r o 


dA 


(4tt) a t 


(B-8) 


where dA is the element of the reflecting area. 


In evrluating (B-8) two approximations can be made: 

• For the ground scenes considered in this report, 
can be considered a constant over the area A. 
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• The error caused by assuming r = Rg over a resolu- 
tion cell for a = 48° , will amount to an error of 
0.02 dB in the value of (or a ). 

With the preceding assumptions the radar backscattering 
cross section is given by (reference 3) : 

1 

f(0,<j)) dA 
A 

(B-9) 

where 

L^yL^ = path losses through intervening medium for trans- 

mission and reception 

= one way antenna gain in the direction of the 
antenna bores ight 

f(6,ip) = two-way antenna gain pattern in any direction 

(figure B-3) specified by angles 6 and (p . 

The value of is to be related to the surface param- 

eters of the remotely sensed scene. For this purpose the 
location of the field-of-view (FOV) has been calculated for 

each value of a . The details of this calculation are 

o 

given in Earth Resources Production Processing Requirements 
for liRHP Electronic Sensors document (reference 40). The 
values of R and the local angle of incidence 6 are also 
calculated using the FOV program and the SKYBET tape. 


%C0i,<l>i) = 


(^”5^ • "e „ r 1 

^"2 ^^ 1 ^ 2 ^ P ^ 
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The approximate value of the antenna pattern integral in 
equation (B-9) is given by (reference 54) : 


/■ 


G^ce^,*^) fCe,(|.)da = sec a 


(B-10) 


In equation (B-10), is the nominal altitude of the 

Skylab (taken as 435 kilometers for calculating 

is the gain of the S-193 antenna. Equation (B-10) assumes a 

pencil beam antenna. The University of Kansas has performed 

integrations for the antenna pattern equation (B-10) . These 

values are dependent on the receiving and transmitting 

polarization states of the antenna. For SL-2 and SL-3, the 

2 2 

values of = G I H are shown below: 

C 0 0 0 


Transmitting antenna Receiving antenna 
polarization polarization 



V 
H 

V 
H 


V 
H 
H 

V 


6.778 X lo"^ 
7.146 X lo"^ 
6.745 X lo"^ 
6.745 X 10^ 


The ratio of (Pj^/P,p) is a function of the sensor 
operating modes and sensor output voltages. This ratio is 
given by: 


P 

P 


R 

1 


K. 


KrK,^, 


(IT)c 

nry^ 



AC (DF) 


V., - V. 


N 


(IT) 

ITT) 


N 


XW7 F,, G 


N 


N. 


V. 


(B-11) 
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In equation (B-11), 


K. 


K 


R 
IT 


scatterometer calibration path attenuation 
loss unique to the receive path 
loss unique to tl>e transmit path 

integration time with subscripts S, N, C to denote 
signal and noise, noise, or calibrate, respectively 
(mode, gain, and angle dependent) 


TC = integration time constant with subscripts S, N, C 

to denote signal and noise, noise, or calibrate, 
respectively (mode and angle dependent) 

~ scatterometer gain for signal, noise, calibrate, 
respectively 


Ps,P^,Fc = scatterometer filter gain for signal, noise, cali- 
brate, respectively (filter, gain, pitch angle, 
and Tjp dependent) 


AC(DF) = angle correction to account for filter attenuation 
encountered in ITC mode at 48° angle for the LCF 
filter because the antenna only- achieved 43°. For 
all other modes and/or angles, AC = 1 . At 48° 
ITC: 

AC = antilogj^Q 
+ 'a^DF^l 


[a^ + a^DF + a2DF^ + a^DF^ + a^DF^ 
for 0.44 < DF < 0.48 MHz 


AC = 1 for DF >. 0. 48 MHz 
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No calculation is made for DF < 0.44 MHz or for 
MCF or HCF filters. 

a^ = 693762.396 

= -7550134.85 

a^ = 32845532.2 

i 

■ a^ = -71395735,3 
a^ = 77541217.8 
a^ = -33662017.0 
DF = Doppler frequency 



[q(Tjp) 

X ITg/TCg] - 

X 

CO 
H 
1— 1 
1 1 

drift ] 

Vs = "'n - 

(T j p) 


[IT^ X 

drift ] 

- 

[qCTjp) 

X ITg/TCg] - 

[IT^ X 

drift] 


Tjp = internal processor temperature A012-193 

q = scatterometer voltage correction constant 

(filtor, gain, and time constant dependent) 

vg = measured signal plus noise voltage 

Vj^ = measured noise voltage (the next value 

following Vg for ITC mode, for other modes 
polarization and command angle of must 

match Vg ) 
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= measured calibrate voltage* (scatterom- 
eter calibration 1 or 2 is used, depending 

I 

on which has the greater ratioj 

drift = integrator drift correction, a function of 
the time constant and filter 


The values of the' parameters K^, Kj^, K,^, (IT) 


(IT)^, (TC)g, (TC)^,, (IT)^, drift, q, 


Vrp y 


'N 




G F 
^’C’ 'S’ 


(IT) 


N’ 

and 


are given in EREP Calibration Data document (reference 55) 
The voltage values of Vg, Vj^, and are contained in the 
measurements A063-193 and A064-193 depending on the scatterom- 
eter status . 


Relative voltages are also computed for the three compo- 
nents signal, noise, and calibration as follows: 

^ ^ [(TC)g/(IT)gl 

^ [CTC)j^/(IT)j^] 

V” = [Vc/Gcl X [(TC)^ X (IT)^] 

1.3 TIME TAGGING AND DATA SEQUENCE 


The S-193 Scatterbmeter operates in several scanning 
inodes and submodes. Details of these modes are given in 
appendix A. During each scan, data is recorded in a particular 
scciuence depending on the mode chosen. The integration times 
depend on the mode. The raw data products (processed using 


*Can be entered by control data. 
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NASA/Data Systems and Analysis Division (DSAD) computer 
program) tabulate the data in the same order as the acquisi- 
tion sequence. The times have been, properly scaled from 
Airlock Module Time (AMT) to GMT. The time used in the raw 
data products will henceforth be called the Data Stream Time 
(DST) . In equation (B-11) proper values of Vj^ are to be 
subtracted from Vg . For instance, if the receive polariza- 
tion is vertical, noise corresponding to the vertical receive 
channel following the signal should be used for . The 

production processing program has been developed on the basis 
of one pass processing. Consequently, the sequence in which 
the computations are done are not the same as the data acquisi 
tion sequence. However, the tabulated scatterometer data in 
product S062-11, reflects actual data acquisition times. 

When relating a particular measurement to the ground 
scene, it is necessary that the coordinates of the illuminated 
area be calculated at the time the measurement was taken. A 
study showed that corrections had to be done to arrive at the 
center of the measurement time (CMT) . Since each scatterom- 
eter measurement is collected for a finite period of time, 
CMT(tj^) should represent the center of this period, accurate 
to a 3-sigma confidence limit of 9 milliseconds; 

t = t - At 
m s 

where 

t^ = starting time of the measurement 
At = half of the data collection period 
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For ITNC and CTNC radiometer measurements: 

At = 130 msec for 48.0° and 40.1° command angles 
At = 66 msec for 29.4° and 15.6° command angles 
At = 31 msec at 0.0° command angle 

For ITNC and CTNC scatterometer measurements: 

At = 4 X (N - G) msec 

G = scatterometer gain setting number (D005A193) 
N = 74 for command angle 48.0° 

N = 64 for command angle 40.1° 

N = 57 for command angle 29.4° 

N = 39 for command angle 15.6° 

N = 23 for command angle 0.0° 

For TTC radiometer measurements: 

At = 18 msec 

For ITC scatterometer measurements: 

At = 4 X (9 - G) msec 

For CTC radiometer measurements : 

At = 18 msec for radiometer/scatterometer 
At = 31 msec for radiometer only 
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For CTC scatterometer measurements: 

At = 4 X (N - G) msec 
N = 14 for radiometer/scatterometer 
N = 11 for scatterometer only 

Data products S062-7, S062-11, S062-16 tab group two, 

S063-1, S063-2, S063-3, S063-4, S063-5, S063-7, and S063-8 

have tj^ on them. All other products have t^ . Product 

S061-1 also has the difference t - t used for each FOV 

s m 

calculation. 

The times given for the statistical data products are 
the times of the first measurement of the sample to be 
averaged. The angles and other data given on the "averaged 
scattering cross section" tabulations are the average value 
of the samples. 

1.4 TABULATED ANGLES FOR SCATTEROMETER DATA 

A particular S-193 measurement is taken for a finite 
amount of time. This time depends on the mode, type of data 
(radiometer or scatterometer), and roll/pitch angles. In 
ITC and CTC mode the antenna angles are varying during the 
measurement period. Before an accurate FOV calculation can 
be made for these modes, roll and pitch angles were inter- 
polated to t c corrected data measurement time (t ) . The 

m 

corrected roll/pitch angles (A^) were computed by using the 
equation (reference 40) : 

= A - (a' - A) (At)/(t’ - t^) 
for |a’ - A| <3° 
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where 


A = corrected roll/pitch 

A = uncorrected roll/pitch 

A’ = previous uncorrected roll/pitch' ' 

At = see paragraph 1.3 

tg = uncorrected roll/pitch time 

tg = previous uncorrected roll/pitch time 

The angles given in this report are the corrected angles 
A^ and corresponding angles of incidence. 

1.5 POLARIZATION LABEL 

The polarization labels applied to the production proc- 
essed data are the same as the data stream. These labels do 
not follow the normal convention in some modes. In literature 
the polarization is defined by the scattering geometry. The 
polarization labels for intrack modes (no roll), are correct. 
The crosstrack (zero pitch) polarizations should be relabeled 
so that vertical (V) is changed to horizontal (H) . In other 
modes where neither pitch nor roll is zero (for example, cross 
track contiguous mode with 15.0°) care must be exercised in 
interpreting the data. In the analysis presented in this 
report normal polarization convention has been followed. The 
polarization labels given in the production-processed data 
have been properly interpreted in the comparisons with theo- 
retical values of backscattering cross sections. 


B-16 



1.6 CORRECTION FOR 48® ITC SCATTEROMETER DATA 


The actual maximum attained angle for the ITC mode is 
approximately 43® instead of 48®. The sharp Doppler filter 
characteristic curve introduces large errors for the actual 
angle attained. These errors have been removed by involving 
proper correction [AC(DF)] factors in the production data 
processing program. The procedure for correcting data 
involves calculating the Doppler frequency using SKYBET data 
corresponding to the attained angle. The attenuation due 
to the Doppler filter was determined by interpolating the 
filter, characteristics. The scatterometer backscattering 
data was then corrected' for filter attenuation. The data 
analyzed in this report has been corrected for the effects 
of Doppler filter attenuation. 

1.7 MISCELLANEOUS SCATTEROMETER DATA 
PROCESSING REMARKS 

For the f ield-of -view calculations the SKYBET tape 
computations assume a "perfect Z-local vertical" vehicle 
attitude whenever this data was not available. 

The accuracy of the EREP pointing has been determined 
to be 0.7® per axis 3p ( p is the standard deviation). 

There are some data dropouts in the S-193 production- 
processed data products. Wherever it was important for 
Sensor ' erformance Evaluation, the raw data was used to 
calculate the scattering cross sections. 

Calibration data used in the S-193 data processing was 
taken from the S-193 acceptance test data. The range of 
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temperatures used during these tests were not as wide as were 
encountered during the Skylab data takes. Interpolations and 
extrapolations were doneto obtain the values for the tempera- 
ture-dependent variables. The effect of these is not serious 
since ratios are used in. the computation of . 

The antenna scan performance differed from that before 
launch. The integration times could also be different from 
those given as a result of the system acceptance test. For 
the case when no signal was received (deep space) , the 
average signal plus noise power density when equated to the 
average noise power density yielded a set of noise integra- 
tion times slightly different from those given in the 
acceptance test data. Details of the scan performance and 
noise analysis are given in S-193 Sensor Performance Report 
(reference 32) . The integration times for the noise were 
revised to reflect the new values as recommended in 
reference 32. 

The backscattering cross section is reported in decibels 
(dB) relative to 1. After the value of a is computed using 
equation (B-9) , the output value from production data proc- 
essing is 

The value is computed by calculating the average value of 
CTq and converting this to decibels.’ 

At the end of a S-193 Radiometer/Scatterometer data 
take, the scatterometer is switched to standby (STBY) position 
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In this configuration, the transmitter is shut off and the 
receiver is still on. The radiometer is switched to STBY 
(standby) approximately 2 seconds later. The scatterometer 
data during its STBY operation at the end of a radiometer/ 
scatterometer is invalid, since no valid signal is received. 

The average values of reported for such data takes is 

also in error. For this reason, average values of were 

computed using valid data avoiding also the use of default 
values in such a computation. 

The data processing equations used in production data 
processing program assumes a linear model for the sensor. 

The acceptance test data (reference 55) shows that for 
power received by the antenna in excess of -70 dBm and less 
than -115 dBm the system is nonlinear. 

2.0 AIRCRAFT DATA PROCESSING 

NASA/JSC, 13.3 GHz scatterometer underflight data was 
acquired during the SL-2 mission. This data will be used 
for comparison with spacecraft-acquired data. 

The 13.3 GHz scatterometer is a continuous-wave Doppler 
radar system, designed to measure reflectivity per unit area 
as a function of the angle of incidence (9). The scatterom- 
eter antenna illumiiiates a fan-shaped area (approximately 
120° along the aircraft flightpath) , and the data is gathered 
for vertical'- transmit , vertical -receive polarization states 
only. As a result of the forward motion of the aircraft, 
Doppler frequency shifts are introduced and the signal returned 
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by a ground resolution cell can be retrieved be bandpass- 
filtering at the corresponding Doppler frequencies 
(figure B-4). The returned energy may be separated using 
the Doppler equation as a function of incidence angle 

2V, 

fp = sin 0 (B-12) 

where 

fp = Doppler frequency 

= aircraft ground velocity 
X = wavelength of the transmitted powt^’" 

0 = angle o‘ ncidence 

The returned energy is received from alx angles of 
incidence simultaneously and is divided equally into two 
channels, one of which is 90° out of phase with the other. 

The data for each channel, detected by a direct-rf - to-audio 
conversion technique, is amplified and recorded on an FM 
tape recorder. The fore-and-aft beam data are separated by 
use of a sign sensing technique (reference 56). To calibrate 
the system, a ferrite modulator is used to provide an absolute 
power reference level of the transmitted signal. The 
versus 0 information is obtained by subtracting known system 
losses and aircraft attitude and velocity factors and com- 
paring the remainder with a reference signal level. 

The radar cross section per unit area is given by the 
equation 
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1 


2(4tt)^ Vh^ 

2 

h ' 

(B-13) 

where 

Pj, = transmitted power 

Pj^ = power received in the Doppler window defined by A£^ 

Gt»Gj^ ” transmitting antenna and receiving antenna gain, 

respectively, as a function of 9 (incidence angle) 
and \p (crosstrack angle) 

h = altitude of the aircraft 

ip = crosstrack angle (figure B-4) 



Equation (B-13) may be rearranged for computer calcula- 
tions as 

E. 

0^(0) = RC + 20 log h + 10 log V + 20 log ~ 

BW|5 2 

+ 10 log + R(D) - G^F' (6) + Z(9) 


(B-14) 


where 

RC = radar offset constant 

h = aircraft height 

E^ = average radar data at ith filter 
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Figure B-4. — 13.3 GHz scatterometer resolution cell geometry. 


* average reference data 
BW^ = bandwidth of ith filter 

R(D) = system rolloff 

Z(0) * any system errors which can be determined 

2 

GqF*( 9) = two-way antenna gain 

BWj^ = reference bandwidth 

The radar offset constant RC is computed from the 
following equation 
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RC = 10 log^Q 2(4tt)^ + FMC - 10 log^^ - 30 log^^ A 
where 

FMC = ferrite modulator constant 

A detailed description of the program can be found in 
reference 56. The calibration data and detailed evaluation 
of the 13.3 GHz scatterometer system is given in 
reference 57. 
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